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MONTHLY NOTICES 
OF THE 
ROYAL ASTRONOMICAL SOCIETY 


Vol. 117 No. 5 


MEETING OF 1957 MAY 10 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

Bibhu Prasad Dash, Department of Applied Geophysics, Imperial College, 
London, S.W.7 (proposed by J. M. Bruckshaw); 

William Russell Hindmarsh, University Observatory, Oxford (proposed by 
H. H. Plaskett); 

Sudhirendra Nath Saha, Geophysics Department, Imperial College, London, 
5S.W.7 (proposed by J. M. Bruckshaw); 

David Vyvyan Thomas, Royal Greenwich Observatory, Herstmonceux 
Castle, Hailsham, Sussex (proposed by P. A. Wayman); and 

Neville John Woolf, Department of Astronomy, University of Manchester 
(proposed by Z. Kopal). 


One hundred and fourteen presents were announc@d as having been received 


since the last meeting, including :— 














G. de Vaucouleurs, Discovery of the Universe (presented by the author); 

Royal Observatory, Cape of Good Hope, Cape Photographic Atlas of Southern 
Galaxies (presented by the Royal Observatory); 

R. A. Buckingham, Numerical Methods (presented by Sir Isaac Pitman & Co., 
Ltd.); and 

B. J. Bok and P. F. Bok, The Milky Way, 3rd edition (presented by the 
Harvard and Oxford University Press). 


MEETING OF 1957 OCTOBER 11 
Dr W. H. Steavenson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 

Douglas West Allan, Department of Geodesy and Geophysics, Madingley 
Road, Cambridge (proposed by E. C. Bullard); 

Bart Jan Bok, Mount Stromlo Observatory, Canberra, A.C.T., Australia 
(proposed by H. H. Plaskett); 

Sidney George Coomber, ‘ Belmont ’, Hen Parc Lane, Upper Killay, Swansea, 
Glam. (proposed by B. Featherstone); 

Digby Henry Christ, 744 Warwick Road, Solihull, Warwickshire (proposed by 
M. C. Johnson); 

Vincent Deasy, Ard Faill, Raheny, Dublin, Eire (proposed by H. A. Briick); 
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George Alfred Harding, Royal Greenwich Observatory, Herstmonceux 
Castle, Hailsham, Sussex (proposed by A. Hunter); 

Edgar Alpin Hill, 1 Bayswater Road, Sketty, Swansea, Glam. (proposed by 
B. Featherstone); 

William Aubrey Hill, Highways, Hill Grove, Caswell, Swansea, Glam. 
(proposed by B. Featherstone); 

David Reginald Lynn Jones, Bristol Aeroplane Company, 55 Salisbury Road, 
Maesteg, Glam. (proposed by G. M. Brown); 

Harold Edward Herbert Kirkby-Johnson, 1 Eastfield Road, Royston, Herts. 
(proposed by H. Thomson); 

Muiris Maclonnraic, 53 Mount Prospect Avenue, Clontarf, Dublin, Eire 
(proposed by H. A. Briick); 

John Dennis Pope, Royal Greenwich Observatory, Herstmonceux Castle, 
Hailsham, Sussex (proposed by R. v. d. R. Woolley); 

William Henry Thomas, 1 Trinity Place, Swansea, Glam. (proposed by B. 
Featherstone); and 

Clarence Benton Warrenburg, 5835 N 2nd Avenue, Phoenix, Arizona, U.S.A. 
(proposed by A. Knight). 


Two hundred and eighty-four presents were announced as having been 


received since the last Meeting, including :— 


S. K. Banerji, Earthquakes in the Himalayan Region (presented by the Indian 
Association for the Cultivation of Science); 

E. M. Burbidge and others, Synthesis of the Elements in Stars (presented by 
the authors); 

A. C. Clarke, The Making of a Moon (presented by Fredrich Muller, Ltd.); 

T. G. Cowling, Magneto Hydrodynamics (presented by the author); 

Y. Hagihara, Stability in Celestial Mechanics (presented by the author); 

H. M. Nautical Almanac Office, Sight Reduction in Tables for Air Navigation, 
Vol. 1 (presented by H.M. Nautical Almanac Office); 

International Astronomical Union, Symposium No. IV (presented by the 
International Astronomical Union); 

W. J. Luyten, A Catalogue of 9867 Stars in the Southern Hemisphere (presented 
by the author); 

M. K. Munitz, Space, Time and Creation (presented by the Free Press) 

M. K. Munitz, Theories of the Universe (presented by the Free Press); 

N. B. Slater, The Development and Meaning of Eddington’s Theory (presented 
by the Cambridge University Press); and 


H. C. van de Hulst, Light Scattering by Small Particles (presented by Chapman 
and Hall Ltd.). 
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THE MOON’S LIBRATION IN LONGITUDE 
Harold Jeffreys 
(Received 1957 June 5) 


Summary 


Analysis of Yakovkin’s results on a possible free libration of the Moon 
in longitude shows that the motion found agrees in speed and phase with the 
forced term whose argument is twice the difference of longitude of the perigee 
and node. If it is interpreted as due to this term 

Y =0°'0002049 + 0°0000009. 
If the motion is a free one 

y =0°0002098 + 0:0000022 
The results imply values of f close to 0-67 and therefore intermediate between 
the two groups of determinations made from the annual libration. 





Search for the Moon’s free libration in longitude had been unsuccessful 
until recently, and y, the dynamical ellipticity of the Moon’s equator, has been 
found from the annual libration. The results, when expressed in terms of the 
ratio f, cluster about o-5 and o°8, and no satisfactory explanation has been avail- 
able. The situation has been completely changed by two important papers 
of A. A. Yakovkin (1952) and K. Koziel (1949). Yakovkin appears to have 
detected the free libration, with an amplitude of about 50”. Koziel points out 
the possible importance of a term whose argument is twice the difference of the 
longitudes of the node and perigee ; for f= 0-662 this would have the same period 
as the free libration, and accordingly it might be large enough to affect the solution. 
‘The constants «, 8, y are defined in terms of the moments of inertia by 

C—B C—A A—B 
eee 


and f=a/8. ‘To a very close approximation 
at+y=f. 
Yakovkin assumes a trial speed of 0 -3318/day and determines a harmonic 


term from the Kazan observations, dividing them into four intervals according 
to observers. ‘The results are expressed in the form 


Libration = P’sin (0°-3318 (t—t)) +p°) . . ..... (1) 
where ft, =J.D.2412000. Uncertainties are standard errors. 
P p 
IQIO-15 5i+11 248+ 15 
1916-31 46+ 10 206+ 12 
1932-42 60+ 13 187+ 13 
1935-45 78+ 16 184+ 14 


A least squares solution is 
p =200°8 + 6-7 — (2:10 + 0°62)(t — 1930°5)/1 year | (2) 
P=55+6°0 rn 
x?=4'4 on 5 degrees of freedom, 
33" 
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The first question is whether the libration found is actually the free libration 
or a forced one due to the term mentioned by Koziel. We note that its speed is, 
in all, 0°-3261+0°-0017/day; and the speed of the forced term would be 
0°:3287/day. The phase of the observed libration at 1930°5 would be 

'200°°8 + 6°7 + 14160(0°-3318) = 13 x 360° +219°1+6%7 . of kee 
That of the forced term would be 212°:0. If the identification with the forced 
term is correct the comparison contributes 3-4 to x* on 2 degrees of freedom. 
Hence this identification is perfectly possible. Since a free libration might 
have had any phase whatever there is strong reason to suppose that it is correct. 

We make two solutions, according as the motion is interpreted as free or 
forced. If it is free it leads to an estimate of y in the usual way. If it is forced 
the amplitude will give an equation for y. Koziel’s form for a forced motion 
can be written 
2°963 y Hain(ht+-«) (4) 

(jn)? — 2-9198y . 
n=13°:176/day. If the motion is free, the denominator vanishes for h = (0°-3261 
+0°-0017)/day ; hence 





— 


y =0°0002098 + 00000022. . . . . . . « (5) 
For the forced motion Koziel gives (p. 78) 
= —0°00001165 radians= —2-403” . . . . . . (6) 
Substituting in (4), with A=0°-3287/day, gives 
y =0°0002049+0°0000009 . . . . . . . (7) 
Exact coincidence of speeds would give y=0-0002131. 
The inclination of the Moon’s axis to the ecliptic gives, according to de Sitter, 
B(1 + 0:0047f) = 00006286 + 0-000001 5(p.e.) 
whence, with the close approximation Bf=8—y and either above value of y, 


B =0°0006266 + 0:0000022(s.e.). 
and 


f=1-—y/B=0°665+0°0035 ...... =. (8) 


Ce ee 


The uncertainty of f still comes mainly from that of y ; but as there still seems 
to be some conflict between results about the value of 8 its uncertainty may be 
somewhat underestimated. Since the uncertainties of 8 and y are nearly in- 
dependent it seems desirable that they should be given separately and that esti- 
mates of y should cease to be expressed by way of f. 

The apparent accuracy of the two solutions for y far exceeds that of any based 
on the annual libration, but they seem to be the only possible interpretations of 
Yakovkin’s data. I think that the interpretation as a forced motion is the more 
likely. Both require a close agreement between the free and forced periods in 
question, so that this coincidence is equally applicable to both; but the agreement 
in phase requires a second coincidence if the motions is free but not if it is forced. 

On this interpretation the data still do not separate a free libration; its amount 
would be unlikely to exceed 10”. 

A possible objection may be that for a close coincidence of periods the first- 
order theory may not be valid and that higher terms should be taken into account 
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by something like E. W. Brown’s theory of resonance. However an amplitude 
of 1’ or so in a pendulum motion would produce a change of free period by about 
I part in 10’, and it seems that any departure from non-linearity must be negligible 
in comparison with the observational! uncertainty. 

Koziel, noting that the amplitude of the forced term would theoretically be 
infinite on a linear theory if f=0-662*, maintains that if a trial value of f exceeds 
this value, any solution found by the method of least squares would still exceed it ; 
and similarly for a trial value less than 0-662. I do not follow this part of his 
argument, since with a trial value of either 0-5 or o-8 this term would be negligible 
in all approximations by corrections linear in the changes of f. But if any use 
is to be made of the annual term in future it would be desirable to analyse the 
data in such a way as to give separate determinations of the annual term and one 
with a speed near 0°-33/day. Some of the series of observations that have been 
used are short enough for a wrong estimate of the amplitude of the latter term 
to produce an important bias in that of the annual term. 


* T get 0-660, presumably owing to a slightly different value of B. 
160 Huntington Road, 
Cambridge: 
1957 June 4 


References 


A. A. Yakovkin 1952, Trans. I.A.U., 8, 231 
K. Koziel, 1949, Acta Astronomica, a, 4, 61-193. 


Note added 1957 August 1.—Professor Yakovkin informs me that there is a 
numerical error in the calculation of P, which he has not yet succeeded in 
correcting. It should be about 100. ‘This would change the estimate of y, on 
the assumption that the motion is forced, to 0-0002081, and f would be 0-667. 
A fresh analysis for the annual term gives f=0-69+0-02. A summary has 
appeared in the Astronomical Circular of the U.S.S.R. Academy, No. 178, 1957 
March 12. 








THE SECULAR ACCELERATION OF THE MOON, AND THE 
LUNAR TIDAL COUPLE 


C. A. Murray 


(Communicated by the Astronomer Koyal) 
(Received 1957 July 24) 


Summary 


The discussion of the modern observations of the Sun, Mercury and 
Venus by Spencer Jones, when taken in conjunction with Jeffreys’s theory 
of tidal friction, is shown to lead to a total dissipation of energy in the oceanic 
tides which is three times that calculated by Jeffreys in his discussion of tidal 
data for shallow seas. 

The true (negative) acceleration of the Moon corresponding to this value 
of the dissipation has been included in the expression for the Moon’s tabular 
mean longitude in the Improved Lunar Ephemeris. An investigation of the 
value of this acceleration in ancient times using Hipparchus’ eclipse and 
equinox observations leads to a value twice as large as that given by the 
modern observations. This result is very satisfactorily confirmed by the 
observed magnitudes of partial lunar eclipses recorded in the Almagest. 

It is pointed out that if the acceleration does change, then ephemeris time 
derived from lunar observations will not be uniform in the Newtonian sense. 





1. Let the observed correction to the gravitational mean longitude of the 
Moon be expressed in the form 


85L=a+bT+(q+s)T?+B(T) (1) 
where q is the real (negative) acceleration due to the reaction of the lunar oceanic 
tidal couple, s is the apparent acceleration due to any secular retardation of the 
Earth’s rotation and B(7) is the fluctuation due to irregularities in the Earth’s 
rotation; a and 6 are corrections to observationally determined constants of 
integration. 

If N is the lunar oceanic retarding couple, M, m the masses of the Earth and 
Moon, and c the mean distance between them (all in c.g.s. units), and if we assume 
for simplicity that the Moon moves in a circular orbit in the plane of the Earth’s 
equator, then Kepler’s third law and the principle of conservation of angular 
momentum give the rate of change of the Moon’s orbital angular velocity as 


M 
—3N Tee radian/sec?. 





‘Thus, if q is expressed in seconds of arc and T in centuries, 


3 


q= - *Naa3 x 2°05 x 1074, 
Putting M=5-98 x 1077, M/m=81-5, c= 3-84 x 10", we have 
q= —2°88 x 10-8 N, (2) 


2. Jeffreys (1) has pointed out that, in the case of the solar tidal couple N’, 
the corresponding real acceleration of the Sun is observationally quite negligible. 
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Accordingly the observed correction to the gravitational mean longitude of the 
Sun may be expressed in the form 


dL’ =a' +b'T +(n'/n){sT? + B(T)} (3) 
where n, n’ are the mean motions of the Moon and Sun. _ In this case the apparent 
departure from a purely gravitational orbit arises solely from non-uniformity 
in the Earth’s rotation. 

3- In his paper on the rotation of the Earth, Spencer Jones (2) adopted 
definite numerical values for a, b and (q+5) in (1) and thus, from observations 
of the Moon’s longitude 6L, defined B(T) as an empirical function of the time. 
By eliminating B(7) from (3), and from similar equations for Mercury and 
Venus, he determined (n’/n)s from a discussion of observations extending from, 
roughly, 1680 to 1930. 

He actually adopted g+s= +5”:22 and derived (n'/n)s= +1":23 +004. 
Since n/n’ = 13°37, his discussion gives 

q= —11"-22 +0753. (4) 
Thus, from (2), 
N=(3'9+0°2) x 10% dyncm. 
Inserting this in Jeffreys’s equation for the dissipation (oc. cit. p.227), and taking 
N/N’' =3°4 we have 

dE/dt = (3-5 + 0-2) x 10! erg/sec. (5) 
This is nearly three times the rate of working calculated by Jeffreys from his 
discussion of tidal data; he remarks, however, that there is probably some dissipa- 
tion along open shores which he did not include in his calculations. 

4. The value of q derived in (4) is of course the correction which has been 
applied to the theoretical secular acceleration in constructing the Improved Lunar 
Ephemeris (3), as originally proposed by Clemence (4). It is an empirical correc- 
tion derived from observations extending over the last three centuries. There 
is no a priori reason why it should remain constant over long periods of time, 
and it is desirable to derive the value which best represents the available ancient 
observations. 

In (1), sT?+ B(T) is the apparent displacement of the Moon in its orbit due 
to non-uniformity in the Earth’s rotation. It can only be determined in practice 
from observations of the Sun or planets. (It clearly cannot be determined from 
observations of the Moon unless g is known a priori.) In order to derive a value 
of g it is thus necessary to have observations of the Moon and Sun (or planets) 
at as nearly the same epoch as possible. 

If 5L, SL’ are corrections to the tabular longitudes of the Moon and Sun at 
distant epochs 7',, T, respectively, the average sidereal secular accelerations v, v’ 
over the whole periods 7,, 7, are defined by 


v=T, SL, v=T, SL’ 
‘Thus from (1) and (3), neglecting a, 5, a’, 6’ 
g=v—(n/n')v' (6) 
provided that 7,*B(7,)— 7,?B(T;) is negligible. We assume that this will 
be the case if 7, and T, do not differ by more than a few years. 


5. The observations of Hipparchus in the second century B.C. afford perhaps 
the most reliable determination of q in ancient times. His series of equinox 
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observations has been discussed by Fotheringham (5) who derived 


v= +1795 +0"-27 (7) 
for the correction to Newcomb’s sidereal acceleration of the Sun for mean epoch 
—137. His observation of the solar eclipse of —128 Nov. 20 has also been 
thoroughly discussed by Fotheringham (6) who derived the following equation 
of condition for corrections to the sidereal accelerations of Brown’s and Newcomb’s 
tables : 


— rae +2”-10+0"°55. (8) 
Combining (7) and (8) we find from (6) 
q= —20°-6+ 3°-2. (9) 


6. As has been remarked by Jeffreys (Joc. cit. p. 224), it is also possible to 
obtain g directly from the magnitudes of lunar eclipses. Apart from the small 
effects due to differences of topocentric libration, the magnitude as seen from 
anywhere on the Earth depends only on the geocentric co-ordinates of the Moon 
and Sun. It is thus independent of U.T. 

If Ay, A’, are the ecliptic longitudes of the Moon and Sun, f, the latitude ot 
the Moon, and X, 2’, f their respective rates of change at an E.1’. fy, say, sufficiently 
close to opposition, the relative co-ordinates of the Moon and shadow at time f 
may be written 

A—A' + 180° =Ay— Ag + 180° + (A—Y’)(t — 4), 
B= By + B(t—t). 
The angular distance between the centres of the Moon and shadow is 
{(A—A’ + 180°)? + 87}? which has the minimum value 


Consider the change in this minimum value corresponding to a change g7” 
in the Moon’s mean longitude. We have 


Bo ape —Not 180')|, (10) 











a+BA=41 + apt By — wep teXot 180 )+ar*(55 x a) 
Now let ¢, be the actual E.'l’. of opposition; then (11) and (10) become mC 
add 11+ Po a+ 384) (12) 
oe bala Fb (13) 

where f, is now the latitude of the Moon at opposition and 
58, = at (3 - Pear): (14) 


For an eclipse to be partial, |8,| must be at least of the order 26’ whereas, in the 
eclipses discussed below, the maximum value of |8fp| is 6’-2. Accordingly 
|8o|> |88)| and we can write 


SA= + {1+ a eeat a (15) 


according as f, is positive or negative. 
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It can readily be verified from numerical formulae given by Cowell (7) that 
f B® \-2 / 0B Baa 
| ° (A-2X’?J ob (X-X)dL 
where g is the Moon’s mean anomaly and the upper or lower sign is to be taken 
according as the eclipse is at the ascending or descending node. Also, if a is the 
Moon’s semi-diameter, then 
ao  942"(1 + 0-064 cos g). (17) 
If 5G is a small change in the magnitude of an eclipse, we have 


) = ¥ 0°:0076(1 — 0075 cosg) (16) 


. I 
8G= - = 8A (18) 


or, combining (14), (15), (16) and (17), 

6G= +40 10° %(1—o 14c0sg)q7*. (19) 
‘The upper or lower sign in (1g) is to be taken according as the Moon is eclipsed 
after or before passing the node. 

7. The observed magnitudes of eleven partial lunar eclipses recorded by 
Ptolemy in the A/magest have been compared with tabular magnitudes by Cowell 
(8) who took as first approximations to the secular accelerations v,= + 4"-88, 
v= +4711, or from (6) ¢.= —50”. The observed minus tabular magnitudes 
are given in ‘lable I, together with the values of g and the epochs 7 measured in 
centuries from 1800. These have been taken from the tables on pp. 524-7 
of Cowell’s paper with the exception of the magnitude of eclipse No. 18, for 
which Fotheringham’s corrected value has been taken (9). 


Taste I 
Ref. no. Fi g +4 10°78 T? (1— 0°14 cosg) 8G R 
2 25°19 194 000289 O'21 0°13 
3 25°18 345 218 O07 ‘Ol 
4 - 24°20 162 264 03 "10 
5 -23°21 210 241 “50 "04 
6 ~23°00 184 242 02 05 
¥ | — 22°90 281 ~ 203 +- OI + +06 
14 —19°73 344 ‘i 135 T SF — a 
15 —19°'40 359 + 129 + +07 + *04 
16 — 16°75 72 — 108 — ‘02 + ‘or 
18 —16°65 245 117 03 "00 
19 —16:64 326 0:00098 “00 0°03 


A solution by least squares, giving unit weight to each observed magnitude, 
gives a correction of +27”+6" to q,. ‘The result from the lunar eclipse 
magnitudes is then 

q= —23"+6". (20) 
The final column, R, of Table I contains the residuals from this solution. ‘The 
predominance of positive residuals suggests that the observed magnitudes may 
be systematically too large by about 0-03; this is in addition to an allowance of 
0-02 for increment in the radius of the shadow due to the Earth’s atmosphere, 
which has already been included in 8G. The probable error of one magnitude 
is + 0-04 and the magnitudes were in fact only recorded to the nearest digit (0°08). 
In view of the indefinite nature of the phenomenon the presence of a small 
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systematic error is hardly surprising. In any case its effect on the final result 
(20) is inappreciable. 

8. The value of q derived from the ancient lunar eclipses (20) is in excellent 
agreement with that deduced from Hipparchus’ observations (g). Combining 
these two independent determinations we may adopt 


q= —21"+3" (21) 
as the value which best fits the ancient observations discussed in this paper. 
This differs from the modern value (4) by nearly three times the sum of their 
probable errors. There is thus strong evidence that the lunar tidal couple 
may have changed considerably during the last twenty centuries. This is not 
improbable. The dissipation is thought to occur in relatively few regions of 
shallow sea ; local changes in coastline and sea level in these regions will therefore 
have a large effect on the total dissipation. 

The latest solution for the secular accelerations involving ancient observations 
is that by Brouwer (10); he, however, imposed the condition that g should have 
its modern value (4) throughout the whole period covered by his discussion, 
the validity of which procedure would now appear to be open to doubt. 

Apart from changes in the mean longitude at epoch and the mean motion 
of the Moon, the revision of Brown’s mean longitude adopted in constructing 
the Improved Lunar Ephemeris consists of the removal of one empirical term 
(the “‘ great empirical term’’), and the substitution of another (q7?). This 
ephemeris is to be used as the standard of ephemeris time, all departures of the 
Moon’s observed position from the ephemeris being attributed to deviations 
of the rotation of the Earth from uniformity. Any change in q will therefore 
introduce a systematic departure of E.'T. from a truly uniform Newtonian time. 

I wish to thank Dr R. d’E. Atkinson for his helpful criticism of this paper. 


Royal Greenwich Observatory, 
Herstmonceux Castle, Sussex : 
1957 June 28. 
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THE EFFECT OF BLANKETING ON THE STRUCTURE OF THE 
SOLAR ATMOSPHERE 


Anton’ Przybylskt 
(Received 1957 April 2) 


Summary 


The effect of blanketing on the structure of the solar atmosphere has been 
investigated by computing a (physically impossible) model solar atmosphere 
with variable net flux of energy such as would result if the absorption in spectral 
lines were to cease while retaining the actual temperature distribution of the 
solar atmosphere. Pecker’s (Ann. d’Ap., 14, 152, 1951) values of the fraction 
7 absorbed in spectral lines have been taken as the basis for these investigations. 

The investigations give a temperature drop above the level of mean optical 
depth 7 0-088 and a rise of up to 225 degrees below this level. 





1. Introduction.—In a previous paper, which subsequently will be referred to 
as (1), a model solar atmosphere with practically constant net flux of energy was 
computed. ‘Thus the computed model satisfied the important condition of radia- 
tive equilibrium, but in spite of that the agreement between the observed limb 
darkening and the limb darkening computed theoretically from the model was poor. 
This disagreement between theory and observation is mainly due to the following 
three factors : 


(i) the model did not take into consideration the effect of the absorption lines 
(blanketing effect), since only the absorption in the continuum was taken into 
account ; 


(ii) there may be a departure from radiative equilibrium due to convection in 
the unstable zone ; and 


(iii) the adopted values of the monochromatic absorption coefficients may be 
erroneous. 

In the present paper it is intended to investigate the modification of the struc- 
ture of the solar atmosphere due to the blanketing effect. An investigation of the 
other two sources of disagreement between theory and observation will be the 
subject of a subsequent paper. 

2. Previous investigations of the blanketing effect.—\n his classical paper on the 
blanketing effect, Milne (2) assumed that the spectral lines were formed in an 
optically thin outer layer of the solar atmosphere (reversing layer), and that they 
absorbed a certain fraction of the radiation and returned it to the photosphere. 
This reflected radiation could not accumulate indefinitely but escaped in the form of 
increased radiation in the continuum between the lines because of the increase of 
temperature in the photosphere under the reversing layer. ‘The mathematical 
treatment of this problem gave an overall increase of the temperature in the solar 
atmosphere. 
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Chandrasekhar (3) assumed that the absorption lines were uniformly distri- 
buted in the whole spectral range and were formed throughout the whole atmo- 
sphere, avoiding in this way the artificial division of the solar atmosphere into a 
photosphere and a reversing layer. Assuming the ratio of the absorption coeffici- 
ent in lines to the absorption coefficient in the continuum to be constant and 
independent of the depth and using Eddington’s approximation, he solved the 
equation of transfer separately for the continuum and for the absorption lines under 
the assumption of a constant net flux of energy. His investigation showed that 
(i) true absorption lowered the temperature close to the boundary of the atmosphere 
and increased it in the deeper layers; and (ii) pure scattering did not change the 
boundary temperature but increased the temperature gradient in the stellar 
atmosphere. In addition, his treatment of Milne’s problem of a reversing layer 
showed an increase of the temperature gradient in the reversing layer without any 
alteration in the gradient beneath it. The temperature below the reversing layer 
was, however, raised. All these results were obtained starting from the Milne- 
Eddington model of the solar atmosphere. 

Hopf (4) supplemented the theory of Chandrasekhar by evaluating exactly 
the boundary temperature using the method applied in his tract (§) for the exact 
solution of the equation of transfer in a grey atmosphere. Miinch (6) solved the 
problems dealt with by Chandrasekhar using the methods developed by 
Chandrasekhar (7, 8, 9) for the approximate evaluation of the source function by 
means of Gaussian sums. In addition, he generalized the solution of 
Chandrasekhar for the case of pure absorption assuming lines to be unevenly 
distributed throughout the spectrum and applied the results of his investigations to 
the solar atmosphere. 

All the above attempts, with the exception of Miinch’s treatment, were of 
little practical value since the assumptions concerning the nature of absorption in 
lines and their distribution in the spectrum were too simplified and artificial. 
Pecker (10) and Labs (11) applied a method proposed by Strémgren and Unséld 
(12) for correcting a model stellar atmosphere to models which probably approach 
closely to the real solar atmosphere. This method gives the corrected tempera- 
ture at a given level from the assumption that the absorption in the continuum and 
spectral lines must be equal to the total emission 


I (Ke+,)I, dv= f(x, +«,) B, dv. (1) 


Laborious computations were necessary in order to find the absorption coefficient 
x, in spectral lines. It was, of course, not possible to deal with every line separately 
but good results were obtained where lines were grouped in classes according to 
their equivalent width, ionization potential, exitation potential and wave-length. 
In this way the absorption lines were classified into a number of line types and the 
work was greatly reduced. The method of improving model stellar atmospheres 
by means of the iterative equation (1) is, however, only slowly convergent and for 
this reason Labs (11), applying the formula only once, was not able to improve the 
temperature distribution below the mean optical depth 7 = 0-01. 

Pecker (10), assuming that formula (1) need be applied once only, and making 
other simplifying assumptions, derived the formula 


T( *) f d (+ —— T.\4 
7h | aad ss = oe OC ieee oF J) ae 
1 r (ATs F)| Jx, B, dv/{B, dv (7) . (2) 
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where AT is the required temperature correction, F the constant net flux of energy 
and AF the fraction of the flux absorbed in lines. Using this formula he calcu- 
lated corrections down to the optical depth t5999=0-9. However, since Labs, 
using the exact formula, was unable to obtain satisfactory corrections below 
7 =0°O1 it is extremely unlikely that Pecker’s corrections based on an approximate 
formula are in any way reliable. 

The disappointing results of Labs’ investigations prompted Bohm (13) to 
apply the “‘flux-iterative’’ method of Unséld (14) to the problem of blanketing. 
This method is based on the equation 

—7.AB=3AF(o)+ 3 | AF(@)d(r) 4-5 (AFC) (3) 
in which AB is the required correction of the source function at the mean optical 
level s and AF(r) is the difference between the actual net flux of energy and the 
required constant net flux.* Formula (3) has been derived by R. v. d. R. Woolley 
(15) for a grey atmosphere for which Eddington’s approximationis valid. It gives 
only an approximate correction AB but a better value can be found from it with the 
aid of the exact iterative formula 


1d 
AB=A(AB)— — (AF) (4) 
in which the A-operator is defined by the equation 
A, UfO}=3 | AOE (t—x)dt+3 | (OE (x—dad. 


‘The solution of equation (4) must satisfy the boundary condition 
AF(0) = 7®, (AB) 
where the ®-operator is defined by the equation 


, {f(t)}=2 | “ f(t) E(t—x) dt—2 | £0) E,(x—t)dt. 


The method of Unsdld gives good results since the iterative formula (4) is likely 
to be convergent for reasons given by Bohm in his paper. In the present paper, 
however, it is intended to apply another method to the blanketing problem and to 
compare the results with those of Bohm. 


3. Application of a variational method to the problem of the blanketing effect 


3.1. Introductory remarks.—The blanketing effect modifies the temperature 
distribution, which in consequence changes the continuum between the spectral 
lines. If, then, the absorption in lines were to cease while the modified tempera- 
ture distribution were maintained, the total flux of energy in the solar atmosphere 
would not be constant but would vary with the depth. The actual value of the 
flux F'(r) at the level t would be the observed constant flux F’,,.(o) increased by the 
amount of energy absorbed in spectral lines at thislevelz. It is, therefore, possible 
to allow for the effect of blanketing on the structure of the star’s atmosphere by 
computing a model of the solar atmosphere with variable total flux F(z) dependent 
on the optical depth and disregarding the spectral lines altogether, i.e. by taking 
only the continuous absorption coefficient into account when computing the net 
flux of energy. ‘This task of computing a model solar atmosphere with a variable 

* In the present paper the symbol F for the net flux is as used by Woolley and Stibbs in The 
Outer Layers of a Star (Oxford, 1953). In the Radiative Transfer of Chandrasekhar (Oxford, 1950) 


and the Basic Methods in Transfer Problems of Kourganoft (Oxford, 1952) 7. F is used for this 
quantity, 
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flux of energy is a priori more difficult than the computation of a model with constant 
flux. However, the application of a variational method described by the author 
in (1) makes this task very easy if a model solar atmosphere with a constant or 
nearly constant flux of energy is available. Such a model has been published in 
‘Table V of (1): the net flux of energy (without taking the spectral lines into account) 
varies by 4:3 per cent from its boundary value F®(o) = 62-68: 10° ergs cm™* sec! 
within the mean optical depth r= 11-528. In the present paper it is intended to 
improve this model which subsequently will be referred to as the initial model. 
The designation Fr) will be used for its net flux of energy in accordance with the 
notation used in(1). The values of F®(r) are shown in Table VI of (1). 

3.2. Theoretical remarks.—In the variational method of improving stellar 
models described by the author in (1), the required corrections of the source 
function AB(r) are represented as a linear combination of some elementary 
functions f,, fo,....- ye 


AB (7) =a, f,(7) + Gofo(7) + «+--+. +a, f,(7)- (5) 
The resulting change of the net flux is then 
AF (r)=a,7 ®,(f,)+a,7®,(f2)+ ..... +a,7®(f,). (6) 


Assume now that the change of the net flux A F (7) is equal to the difference between 
the required (variable) flux F(r) and the net flux F®(r) of the initial model. The 
system of equations (6) written for a sufficient number of optical levels + can be 
solved for the unknowns 4), dy,..... , 4,, which, inserted in equation (5), give the 
correction of the source function AB(7) and consequently the correction of the 
temperature AT(r). 

3.3. Application to the solar model.—An application of the above method in (1) 
was so successful that it could be applied with confidence to the present problem. 
‘The essential point of this method is that the same mean of the absorption 
coefficient is used in the computation of the mean optical depth for the initial and 
the corrected models. It can be expected that the temperature correction will 
have the simplest form if the true mean absorption coefficient 

yg 
“— % ’ x, F dv 
is used with the inclusion of the absorption in special lines. ‘This, however, cannot 
be computed until the final model isknown. An approximation of this mean can 
be computed for the continuous absorption coefhicent with the formula 


I the 
<= FRG, 7 Po Moy ) 


in which the monochromatic fluxes F(z) are known from intermediate results of 
the computations in (1). The results of the computations of the mean (con- 
tinuous) absorption coefficient « with the formula (7) are shown in Table I. For 
reasons explained in (1) the values of the logarithm of the electron pressure p, are 
given together with the values of x. 

According to Béhm (13) the absorption in spectral lines increases the mean 
absorption coefficient by about 10 per cent, and therefore the values of « from Table I 
multiplied by 1-1 will be used for the computation of the mean optical depth, which 
will be designated +* for both the initial and the corrected (final) model. The 
notation 7 will be used for the mean optical depth without taking the absorption in 
spectral lines into account. Since the optical depth of the initial model has been 
computed in terms of another mean absorption coefficient a re-computation of the 
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mean optical depth 7* has been made for the purpose of the present paper. ‘The 
results of this re-computation are shown in Table II where both the initial and the 
final model are shown side by side. 

In order to avoid misunderstandings it should be stressed that the correction 
of the mean absorption coefficient due to spectral lines is used only for the com- 
putation of the reference depths in equations (5) and (6). (In the following 
sections these equations will be used with the notation +* for the optical depth.) 
If a check is made to ascertain whether the corrected model has the variable net 
flux of energy F'(r), as mentioned in section 3.1, only the continuous absorption 
coefficient should be taken into account. No check of that kind, however, will 
be made in the present paper, since the results of (1) show that the variational 
method of correcting model stellar atmospheres applied here is a reliable one. 


TABLE | 


Values of «/p, in units of 10°** cm* dyne~ per one heavy particle as defined by Vitense (16) 








© logio Pp, Kip. 
1°04 0°149 9°989 
1°00 0°599 7°524 
ogo 1°O74 4°936 
o'So 1-647 3°164 
, 0°70 2°3792 2°113 
0°60 3°137 1°547 
0°50 38908 1°434 

















3.4. Fraction n of energy absorbed by spectral lines.—'lhe exact determination of 
the amount of energy absorbed by spectral lines is a.very difficult task, since it is 
almost impossible to draw the continuous background in some parts of the 
spectrum. For this reason the evaluation of the fraction absorbed varies con- 
siderably from author to author. ‘Thus for instance Mulders’ determination (17) 
gave 7 =0-083, Michard’s (18) 7» =0-'124, Wempe’s (19) »=0°0g1. Still less is 
known about the variation of this fraction 7 with the optical depth. Laborious 
computations have been undertaken by Pecker (10) and Labs (11) in order to 
determine this variation, but of course no great accuracy could be expected. Only 
the results of Pecker have been published and they will be taken as a basis for the 
present investigation. ‘The uncertainty connected with that kind of investigation 
is so great that an independent re-computation would hardly be justified. ‘This 
means, of course, that the modification of the structure of the solar atmosphere due 
to the blanketing effect can only be computed approximately. 

3.5. Numerical results.—The following eight elementary functions 


Ai(t*)=1, fAt*)=7*, f(7*)=7™, h(r*)=7", 
frr=e" fdrry=e™",  ffr*)=e™, file) = Br") 
have been used in order to find the corrections to the temperature distribution of 

the original model. ‘The ®-transforms of these functions are : 


O(f,)=2E;(7*), ® (f.)=$—2E,(r*), O(f,)= §7* + 4E;(7*), 

®(f,)=4°8 + 47*?— 12E,(r*), O(f5)=2{e-* (0°72963 —Inz*) — E, (7*) — E,(7*)}, 
® (f,) = 0-08 e*”* {9°59453 + E, (47*)} — 0°08 {5 (7*) +E, (7")}, 
®(f;) = 0-02 €-!”* {19°79932 + E, (97 *)} — 0-02 {10E, (r*) + E, (7*)}. 





| 
| 
( 
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The numerical values of the ®-transform of the function f,(7*) = E,(7*) can be 
found in Kourganoff’s handbook (2¢) up to 7* = 3°25. For greater values of +* 
they were computed numerically by the author. 

Equation (6) has been written for 21 levels of +* from o to 12°807. ‘These 
levels were not uniformly distributed but there were more levels for small values of 
+* where the effect of blanketing is more important. Only seven levels were 
taken into account fors* >3°0. It was intended to obtain a corrected model with a 
net flux of energy equal to the boundary value of the initial model F(o) = 62-678 
ergs cm™® sec~! corresponding to a solar constant 1°945 cal. cm~ min, 

The solution of the system of equations (6) by the least-squares method gave 
the following results: 


—- a, = +2°837, a, = — 0°348, a,= +0-0118, 
a,= + 13°888, a,= + 10°483, a,= — 11°320, a,= — 16-878. 


With these coefficients equation (5) gives the corrections of the source function AB 
andthusthe newtemperature. For7*=othe temperature is lowered by more than 
1000 degrees; it remains unchanged for +*=0-088, and then increases by an 
amount ranging up to 225 degrees. ‘The sharp drop of the temperature close to 
the boundary is, however, less pronounced if fewer elementary functions f,, are 
used in equations (5) and (6), even though little is changed for levels +* over 
0:2. ‘Thus if only the first seven functions are used the boundary temperature is 
lowered by only 50 degrees ; for the first six functions there is even a small increase 
of the boundary temperature. _ If only the first five functions are used, the solution 
gives an overall increase of the temperature by about 200 degrees. If, however, 
an additional function f,(7 *) = E; (7 *) is used the system of equations (6) with nine 
parameters a@,, is undetermined and this means that eight functions are sufficient 
to determine the temperature distribution of the corrected model solar atmosphere. 

It must be stressed that the method used in this paper does not permit the 
determination of the temperature distribution in the uppermost parts of the solar 
atmosphere, where a sharp drop of the temperature can be expected. ‘This is due 
to the fact that a large change of temperature in a sufficiently thin layer of the 
atmosphere produces only a negligible change of the net flux of energy. On the 
other hand the corrections in the deeper layers (7 * >o-2) are determined correctly. 

Since the observational data (21-25) indicate a low boundary temperature of 
the Sun (down to about 3700 degrees) the temperature distribution obtained with 
eight functions f,,, lowering the boundary temperature to 3708 degrees, has been 
assumed for the computation of the corrected model solar atmosphere. 

4. Model solar atmosphere with blanketing effect—The method of computing 
corrected model stellar atmospheres has been described by the author in (1) and 
the same method has been used here. The results are given in Table II together 
with the initial model. For both models the optical depth 7 is given in terms of the 
mean continuous absorption coefficient in the first column; a uniform increase of 
the absorption coefficient by 10 per cent is shown in the next column and allows 
approximately for the change of the mean absorption goefficient due to the 
absorption lines. ‘The mean continuous absorption coefficient for the computation 


of 7 has been taken from Table I; for the values of @ = a > 1°04 Rosseland’s 


mean absorption coefficient has been used. 
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Model solar atmosphere 


Blanketing effect on structure of solar atmosphere 


Corrected model 


A 


XL 
‘QI 
‘05 
13 
*2!I 
-20 
-30 
34 
‘37 
“40 
“42 
“47 
“50 
‘54 
56 
“59 
‘64 
‘68 
72 
75 
“SO 
84 
“$6 
“SO 
‘QI 
‘93 
“95 
‘97 
‘99 
ae) e) 
‘Ol 
03 
04 
05 
‘06 
‘06 
‘07 
‘07 
‘08 
‘08 
‘08 
‘08 


(p indicates total pressure due to hydrogen alone) 
Initial model 
Y, is T logioP logiop, T 

0°00 0°00 4511 XK 20 3708 
o'o! o'oll 43830 3°35 0°02 4076 3 
0°02 0°'022 4549 4°04 o'160 4305 4 
0°03 0°033 4867 4°14 0°25 4459 4° 
0°04 0°044 4885 4°20 0°32 4011 4 
0°05 0°055 4903 4°20 0°37 4723 4 
0°06 0:066 4921 4°30 O'41 4814 4 
0°07 0°077 4939 4°34 0°45 4591 + 
0°08 0088 4956 4°37 0°45 4958 4 
0°09 0099 4973 4°40 0°50 5018 4 
o'10 O*110 4990 4°42 0°53 5070 4 
oO'12 O°132 5025 4°46 0°55 5154 4 
O'l4 O°154 5058 4°50 0°62 5222 4 
o'16 0176 5091 4°53 0°66 5279 4 
j o'18 0'198 5123 4°56 0°69 5327 4 
: 0°20 0°220 5154 4°59 o'72 5368 4 
i O°25 0°275 5231 4°04 °°79 5455 + 
i 0°30 0°330 5305 4°68 o'85 5529 4 
i 0°35 0°385 5377 4°72 o-gl 5595 4 
i 0°40 0°440 5440 “95 0:96 5656 4° 
i 0°45 0°495 5511 4$°75 rol 5714 + 
oS 0°550° 5575 4°50 1°05 5771 + 
06 0660 5694 4°85 I'l 5881 4 
o'7 0°770 5807 4°88 12 5986 4 
o8 0880 5913 4°91 1°32 6085 4 
°-9 0°990 6013 4°93 1°40 6180 4 
1°O 1‘10 6110 4°95 1°48 6271 4 
1°2 1°32 62558 4°95 1°64 6442 4 
I*4 1°54 6452 5°00 1°75 6600 4 
1°6 1°76 6606 5°02 1‘9gI 6748 4 
3 1-8 1-98 6748 5°03 2°02 6886 5 
2°0 2°20 6880 5°04 a°t3 7016 5 
2°5 2°75 7177 5°06 2°35 7309 § 
3° 3°30 7437 5°07 2°54 7566 5 
3°5 3°85 7671 =5°08 = 270 7798 5 
4°0 4°40 7884 = 5°09 2°84 8008 5 
4°5 4°95 3083 5°09 2°96 8202, 5 
5°70 5°50 $270 5°10 3°06 8381 5 
; 6 6-60 8611 5‘10 3°25 708 5 
7 7°70 8923 511 3°41 goo2 5 
8 8-80 9207 5*11 3°54 9270 5 
9 9°90 9469 5§'11 3°66 9518 5 
10 11° 9714 5°12 3°76 97590 = 5 

As expected, the blanketing effect causesa sharp drop of the tem 
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5. Comparison with other models 

5.1. Boundary temperatures.—The method used in the present paper does not 
provide any means of determining the effect of the blanketing on the outermost 
layers of the Sun. Thus the derived boundary temperature of 3708 degrees may 
be considerably out. Pecker (10) and Labs (11) in their investigations of the 
blanketing effect derived a boundary temperature of about 4300 degrees, while 
Bohm (13) found a temperature of 3400 degrees. 
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Fic. 1.—The effect of blanketing on the structure of the solar atmosphere. 
initial model; —- — — — corrected model. 


5.2. Comparison with Béhm’s model.—The method used by Béhm in his 
investigations of the blanketing is, no doubt, correct and therefore a comparison of 
his final model with the present model is of special interest. Since, however, the 
mean absorption coefficient used by him is not identical with that used here, the 
optical depth 7, for A= 5000 A has been chosen as the basis for this comparison. 
As already stated in Section 5.1, his value of the boundary temperature is 300 
degrees lower than that of the present model. Since, however, the temperature 
gradient of his model is much steeper, the temperatures of both models are equal 
(3840 degrees) at ty=0°0007. In the interval from 7,=0°0007 to 7,=0°07 the 
temperature of Bohm’s model is up to 230 degrees higher while in the interval from 
7 = 0°07 to 7,=0°6 the difference is up to 130 degrees in the opposite direction. 
Below the level 74=0°6 Béhm’s model shows a higher temperature and the 
difference reaches 105 degrees for ty= 2-0. Béhm’s model is not given for higher 
values of 7». 

It is difficult to explain the differences between this model and that of Bohm. 
They may be due partly to the fact that Béhm’s model was computed for a solar 
constant 1-96cal cm~* min“, while this model is computed for a slightly lower 
value (1-945 cal cm~* min-'). They may also be due to the use of a lower value 
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for the fraction of energy absorbed in spectral lines y=0-109. In addition, the 
variation of 7 with the optical depth is certainly not the same for both models. 

Anticipating some results of further investigations on the model solar atmo- 
sphere, which will be published soon, it is possible to say that Bhm’s model seems 
to represent the uppermost layer of the solar atmosphere better, i.e. where ty) <0°07, 
while the present model may be better at greater depths. 

5.3. Comparison with empirical models.—Fig. 2 shows a graphical comparison 
of the present model with the empirical model of Barbier (26) derived from the 
limb darkening of the continuum, that of de Jager (27) derived from the limb 
darkening of hydrogen lines and that of Vitense (25) derived from the investigation 
of metal lines. Differences of the order of 100 degrees exist between the four 
models represented on the diagram between 7)=0°0 and t9=1°0. For ty>1°0 
these differences are of the same order of magnitude but the empirical models are 
no longer reliable here. 
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Fic. 2.—Comparison of the theoretical model with the empirical models. 
-_-—— Vitense’s model | —+—+-—+-— de Fager’s model VII 
—r+—+-—+-—- ~~ Barbier’s model — model under discussion 

(In the interval from 0-01 to 0°06 de Fager’s model is almost identical with that of Vitense.) 





6. Comparison with observational data.—A comparison of the limb darkening 
computed from the corrected model and the observed limb darkening shows 
considerable differences, even larger than for the initial model. This means 
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that taking blanketing effect into consideration does not improve the model in this 
respect. 

A further attempt to improve the model in order to remove the discrepancies 
between computation and observation will be described in a subsequent paper. 

7. Application of Unséld’s method.—The application of formula (3) with the 
subsequent use of the iterative formula (4) should give the same corrections of the 
temperature as the variational method employed in this paper. Indeed, for 
+*>o:2 the corrections of the source function computed with formula (3) do not 
differ by more than 8 per cent from those obtained in the present paper. It may 
be assumed that the subsequent application of the iterative formula (4) would 
remove these relatively small discrepancies, although no check has been made of 
this. For 7*<o-2 the differences are larger but there no determinations are 
reliable. 

It may be concluded that Unséld’s method and the variational method used 
in the present paper give practically the same results. 
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TEMPERATURES AND ELECTRON DENSITIES IN FLARES 
AS DERIVED FROM SPECTROSCOPIC DATA 


J. T. Jefferies 
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Summary 


It is shown that model flares with electron concentrations and temperatures 
in the ranges 5 x 10!!—10!cm~* and 104—1°5 x 104 deg. K respectively 
have emissions consistent with the observed hydrogen and helium line 
intensities in solar flare spectra. Some difficulties are shown to exist in 
interpreting the observed great line width of Ha and a possible solution is given. 





1. Introduction.—Over the past twenty years a good deal of observational 
data on solar flares has been obtained from studies both of their radiation at 
visible and radio frequencies and of their associated terrestrial effects such as 
SID’s. Surveys of various aspects of present knowledge are given, e.g. by 
Ellison (1, 2, 3) and Kiepenheuer (4). 

Data on spectra of flares are, however, rather meagre and show considerable 
variation from one observer to another and so, presumably, from flare to flare. 
Comparison of the spectra given by Allen (§), Ellison (6) and Richardson and 
Minkowski (7) shows that this variation is very marked in the helium lines 
D,(3°D-2°P) and A6678(3'D-2'P) ; the former appears, when at all, sometimes 
in emission and sometimes in absorption while the latter is seen only in emission 
from a few intense flares (3). Spectra covering the Balmer series limit, which 
will be of interest in the subsequent discussion, have been obtained by Richardson 
and Minowski (7) and by Suemoto (8). 

The interpretation of flare observations in terms of a single model consistent 
with more than one aspect of the whole data does not seem to have been under- 
taken. ‘The earliest attempt was by Giovanelli (9), who showed that, unless 
self-absorption of flare Lyman radiation were taken into account, temperatures 
of the order of 105 deg. K or more were required to account for the Ha intensity. 
In a more recent investigation Svestka (10) has interpreted the well known 
relationship between H- central intensity and line width in terms of a model 
with temperature ~1-4 x 10'deg. K and hydrogen abundance N,, of the order 
of 10!®cm~* ; Svestka’s analysis, however, assumed thermodynamic equilibrium 
in a flare and his deductions are consequently open to some doubt. Suemoto (8), 
confining his attention to an analysis of the Balmer line H1o, has concluded that 
if 7 =10'deg. K, the electron concentration is 1-1 x 10’ cm-, 

There is as yet no agreement as to the mechanism responsible for the wide 
Hx wings. Ellison and Hoyle (11) and Mustel and Severny (12) have suggested 
a Stark broadening mechanism. Goldberg, Dodson and Miiller (13) employed 
a damping type of absorption coefficient in their analysis although they rightly 
pointed out that the main problem in connection with the H« width is to account 
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for a sufficiently large number of 11-state H atoms in the line of sight. Their 
estimate for this number, ~10!°cm~*, has been derived on the assumption of 
natural broadening but they have made no attempt to link this figure with flare 
thickness, temperature and hydrogen concentration. 

Existing flare models thus seem to account for some observational features 
to the neglect of others. In attempting to gain fuller consistency with the main 
observed spectral features of flares, we are faced with considerable lack of know- 
ledge of the structure of a flare in depth. Probably the most important factors 
in deciding the emission characteristics are the maximum electron concentration 
and the temperature in this region, together with a thickness ; this may be the 
actual thickness for a flare which is more or less uniform in height or a scale 
height if conditions vary through the flare. 

We adopt as a flare model a layer, of uniform electron concentration V, and 
in local statistical equilibrium at an electron temperature 7, parallel to the colder 
solar photosphere which it overlies. ‘The geometric form chosen is unimportant 
in the present context, for the computed emergent intensity will be only slightly 
affected by changes in shape. Absorption in the solar atmosphere above the 
flare is neglected ; measured flare heights indicate that generally there will be 
little absorbing material above this height in the chromosphere. Although we 
can hope to gain only an estimate of conditions in flares by using such a simple 
model, elaborations hardly seem justified in view of the paucity of observational 
data. 

2. Emission intensities in Ha and D,.—If the model flare layer is optically 
very thick at a particular frequency, so that no photospheric radiation penetrates 
the flare, the emergent intensity /,(0, «) in a direction making an angle cos~!u 
with the outward normal can be found from the expression : 


L,(0,n)= | 8,(r)e-"" dr (1) 


where 93,(7) is a source function defined as the ratio of the monochromatic 
emission and extinction coefficients at the optical depth 7, i.e., 


53,(7)=E,/x, (2) 
For either coherent or noncoherent scattering in the medium, the emission 
coefficient E, at the line centre can be written in the approximate form,* 
Toko 


ai ae 5 + € (3) 
where A is a scattering parameter defined so that 1 —A is the fraction of absorbed 
quanta which are subsequently re-emitted in the same line, e, is a ‘‘ true” emission 
coefficient, i.e., it represents the emitted component which is independent of 
~ aan of radiation in the particular line, and ./, is the total intensity defined 

y 


, I. 1,(p) dw. (4) 


| On the assumption that «/« and A are constant throughout the flare, an expres- 
sion for Jy may be found using the Eddington approximation to the transfer 
equation 
d2J 
} de =AJ — 4n€/k, (5) 





* This follows, e.g., from the discussion given by Woolley and Stibbs (14) pp. 164 et seq. 
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where frequency subscripts have been omitted for convenience. Equation (5) 
has the solution, finite at large depths, 
J = 47€/KA + ae~ VO", (6) 
The integration constant « in (6) may be obtained as usual from the surface 
condition J =(2/3) dJ/dr and it is then easily shown from (3), (2) and (1) that 
1(0, 4) = (e/«A)[1 — (1 —A)/{(1 + 2 (A/3))(1 + HV (3A))}]- (7) 
‘lhe quantities €/«A and A may be evaluated quite readily using methods given 
by Giovanelli and Jefferies (15) provided we know the various rates for collisional 
and radiative transitions between the atomic states. Collisional rates for any 
temperature may be found from appropriate cross section data and are given 
for H and He in the author’s publications (16, 17).* Radiative transitions are 
more troublesome since the intensity varies throughout the medium and is in 
any case only found by solving the equilibrium equations and the transfer equa- 
tions. As pointed out in (1§), a simplification is possible if the atmosphere is 
optically either thick or thin in the various spectral regions of interest. In 
the following the flare layer is assumed to be thick in the Lx and Lf lines of H 
and in their counterparts in He, thin in the subsidiary continua, and thick for the 
principal (Lyman) continua unless otherwise specified. It can easily be shown 
from expressions for the optical thickness and ground state populations that, if 
the layer thickness is —10°cm and the helium abundance not improbably low, 
these conditions on the Lyman and Lyman type lines are met when N, ~10!® cm- 
and T’<2:5 x 10'deg. K ; anticipating the results of the present analysis it can 
be said that the derived conditions in flares conform to these limitations. 
TaBLe I 


Intensity of the radiation emitted normally by a thick layer at the centre of Hx. Surrounding 
continuum intensity at centre of disk is unity. Values marked with an asterisk are 
computed for a model optically thin in the Lyman continuum 


T(deg. K) 

N,(cm-*) 170 x 104 1°5 x 10° 2°5 x 10° 
10!! 0°24 o°35" 0°48* 
— 0°72 3 1°9 
10”° 1°8 2°8 4°0 
10!4 2°6 4°4 5°4 

Tasce II 


Intensity of the radiation emitted normally by a thick layer at the centre of D3. Surrounding 
continuum intensity at the centre of the disk is unity 


: T(deg. K) 

N,(cm-*) Io x 104 1°25 x 104 1°5 x 104 2°5 x 104 
ro” 0:28 0°34 0°39 0°50 
10! 0°99 I'l hy | 1°9 
10!3 2°4 28 3°4 4°3 
1ol4 36 4°0 5°4 6-3 


Knowing ¢€/«A and A, (7) may be evaluated and the emergent flare intensity 
expressed as a fraction of the appropriate continuum intensity at the centre of the 
Sun’s disk (Miinch (19)). Values obtained for the central intensities of Hx 
and D, in this way and so applying only to flares optically thick in these lines, in 
accordance with (1), are given as functions of N, and Tin Tables 1 and2. ‘These 


* Some collisional rates given in (16) are too large by a factor of two. Corrections are given 
by Jefferies and Giovanelli (18). 
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figures should be reliable to 50 per cent or better ; at electron concentrations 
~10!2cm- collisions generally predominate and errors in collision rates tend to 
cancel, due to corresponding errors in reverse transitions. 

3. Optical thicknesses in H and He spectra.—The results given in Tables I and 
II apply to layers with high opacities at the line centres. Before making any 
comparison with observation, it is necessary to find expressions for the optical 
thicknesses of the layers in terms of their physical thicknesses 2. 

At the centres of Ha, D, and 16678 the optical thicknesses are given, using 
a Doppler form absorption coefficient, by 

7(Hx) =4:9 x 10°! T-!?.N,(H)= 

7(D,) =6:3 x 10°! T-!2N,3(He)z (8) 

7(6678) =8-6 x 10°" T-!2.N,\(He)z 
where it has been assumed that the hydrogen and helium fine structure states 
of a particular principal quantum number are distributed in the ratios of their 
statistical weights. 

To evaluate the optical thicknesses from (8), values of the I1-state populations 
are required. For He these are given in (17) and for H they can be found in a 
similar manner. Equation (8) then gives the results shown in Tables III, IV 
and V computed on the assumption of an H to He population abundance ratio 
of 5. 

Tase III 


Values marked with an asterisk are computed for a model optically thin in Lc 
7(Ha) z (cm ) 


T(deg. K) 

N(cm-*) 1°o x10! 1°§ x10! 2°5 x 10 
10!! 1°2X 10°° 16x10°**" 66x10°°° 
z0** 4°5 x 1077 ¥°2 xX 107" EX te-** 
10! 73x10°°* 54x 107° 4°0X 107% 
10!! 8-0 x 107-5 gx10°° 42x10 

TABLE IV 
7(D;)/z (cm~') 
T(deg. K) 

NA(cm~*) 1'o x10! 1°25 x 104 ms X 10° 
101! raxse* 2 xX167° Ext 
10"* 3x67" 116x103 2°t X 107° 
so” 50x 10°° 19x 107° 2x10” 
10!4 3 xX 167° 1°9 X 107° 2°2 x 1075 

TABLE V 
T(A 6678)/z (cm~') 
T(deg. K) 

N (cm) 10 x 104 1:25 < 104 1°5 < 108 
1o!! aqxae" 26xso" 24 x10°° 
10!* 29x20" 36x 10° 6:0 x 10“ 
10!" 34xX10°" 54x10" Tex 126-7 
1o!4 3°4X 107% 5°7X10-* 7:0x10-% 


4. Comparison with observation. (a) Emitted intensities—The computed 
results given in the last two sections may now be compared with observational 
data. ‘These are meagre and rather conflicting as far as simultaneous observations 
of Ha and D; are concerned. ‘Thus Ellison (1) gives intensities for one flare only, 
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which in the above units are 3-0 for Hxand 1-2 for D;. Richardson and Minkowski 
(7) recorded D, only for limb flares ; there is evidence, however, that these 
observations did not include any very strong flares. Allen (5) made visual estimates 
of the total changes in the Fraunhofer lines during a flare but his results do not 
lend themselves readily to the present discussion of central intensities. 

Tables III and IV indicate that, particularly at the lower temperatures, the 
optical thickness in Ha will generally be considerably greater than in D,. It 
may be noted that, since an underlying sunspot is not usually seen through a 
flare viewed at the centre of Hx, the flare is normally optically thick at this wave- 
length. ‘There is no observational evidence for the D, optical thickness ; how- 
ever, it is probably never very large since the line is observable only occasionally. 
Thus while the Hz central intensity generally corresponds to that for an optically 
thick model this is not the case for Ds, and so, while the central intensity of Hx 
should be as given in Table I, for D, the intensity should lie between unity— 
corresponding to radiation from the underlying photosphere—-and the value 
given in Table II. 

The closely similar values of the saturation intensities given in Tables I and II 
indicate that, if the D, line appears at all, Hx and D, should appear at the line 
centre either both in absorption or both in emission. ‘The uncertainties in the 
computed data might allow for some departure from this expected behaviour, 
but it is interesting that Waldmeier (20) has found just this type of relationship 
between Ha and D,, the transition from absorption to emission in D, occurring 
with an H~« central intensity of about 1-o. 

It should also be pointed out that, with the noncoherent scattering mechanism 
probably applying for Hx, a double peaked profile will result—Jefferies (21)—as 
observed in most flares, and the maximum intensity across the line could exceed 
that at the centre by a considerable amount—a factor ~2 would not appear 
excessive. For D, this behaviour should not be so marked, since it will show a 
greater degree of coherence in scattering than Hx. Further, since the D, optical 
thickness is in general much smaller than Hx, the double peaks have less tendency 
to form. 

‘The line 46678, which is the singlet counterpart of the triplet Dg, is only seen 
in the brightest flares (1) ; presumably the flare layer is optically very thin in 
this line for all but the strongest flares, and possibly even for them. 

According to Richardson and Minkowski (7), and Suemoto (8), the flare 
emission at the Balmer series limit is small. If this is accepted as applying to 
most flares it may be used as a criterion for limiting their possible range of physical 
conditions. ‘Thus the emission per unit solid angle from an atmosphere thin 
in the Balmer continuum is, at the series limit, 


E,=2:1 x 10° #N 27-3223 ergs cm~ sec~(c/s)-? 
while the background radiation at the centre of the disk is, Munch (19), 
I,= 12 x 10 ergs cm sec~!(c/s)~" 
so that 
E,/I,=18x 10° PN JsT-", (9) 
If this ratio is to be -<5 x 10-® for Balmer continuum to escape observation, then 


N22T-3? <3 x 107" cgs units. (10) 
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We now attempt to determine, using (10) and the results given in ‘Tables I 
to V, the electron temperature and concentration for what might be considered a 
typical class 2-3 flare with low emission in 1 6678 and in the Balmer continuum, 
and an Ha emission at the line centre =1-0 times the neighbouring continuum 
at the centre of the disk. 

Fig. 1 shows the 7, N, space for the region 5 x 10°< 72 x 10* deg. K, 
10!! <N,<10"cm?, together with the zones prohibited to a model flare as 
giving theoretical results incompatible with observation. ‘I'wo values of the 
model thickness, z= 10° and 1o*cm, are used. The left hand curve, corre- 
sponding to the restriction on the Ha central intensity, is uncertain for 
T <10‘ deg. K as theoretical results are lacking ; this part of the curve is therefore 
shown by a broken line. The intensity of and opacity in D, do not materially 
limit the allowed T and N, values and so for clarity its contributions to the figure 
are omitted. 

The main virtue of the results given in Fig. 1 is that they show approximate 
limits to values of T and N, ; these values are subject to some uncertainty arising 
from errors in the basic data used in solving the equilibrium equations and also 
from possible error in the adopted helium abundance. For reasonable variations 
in these quantities, however, there will be only minor changes in the positions 
of the bounding curves. 














—— 3 =10’cm 
16- --—-- 2 =10% cm 
T (A6678) 
“Pr y TOO HIGH 
_ Ba.C EMISSION 
< Wr TOO HIGH 
°c Ho INTENSITY 
[ TOO LOW 
- 0 
ob 
L l l L j 
10"' 5 x10'' = 107? 5x10" 10° 5 x10"! 
Ne (cm~) 
Fic. 1. 





(6) The Ha width.—So far we have not attempted to explain the great width 
of the Ha line observed in solar flares, although clearly this most characteristic 
feature must be accounted for on any valid model. Observational data on the 
Ha line width are quite plentiful and have been investigated, e.g., by Goldberg 
et al. (13) and by Bruzek (22). These authors have shown that the Hx width 
depends strongly on the position of the flare on the disk, wider lines being nearer 
the limb, as might be expected. It appears from their analyses that we need 
to account for lines ~4A wide at the disk centre. The line widths referred to 
have been measured visually and so refer to the difference between the two wave- 
lengths for which the contrast between flare and background is just visible. 
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The estimated line width of a given flare thus depends on the contrast sensitivity 
of the observer, assumed here to be 2 per cent. 
In the H« wings, the flare intensity is given by an expression of the form, 


[P= 1a ~*) m le, ds (1 1) 


where /,”) is the monochromatic photospheric intensity and 7, is the optical 

thickness in Ha measured up from the photosphere and through the flare. We 

shall write r=7" +7° where 7“ is the opacity of the flaring region and 7” that of 

the material below the flare. Since Lyman radiation emitted by the flare will 

excite the underlying material, the contribution 7° may be quite considerable. 

The contrast P between the flare and its surroundings follows from (11) as 
L9-1® — fe,dz 





P,= [@ = [@ a (12) 
which is readily shown to be equivalent to 
(4N, exp(Xy,) —1 | 
r. =1,54 ee Hx) I > (13) 


where (N;/N,) is an average population ratio taken over the emitting and 
absorbing region; X,,,=Av/kT,, T,, being the radiation temperature (6150 deg. K), 
corresponding to the solar continuum near Ha, and 1 —/f,, represents the fractional 
depression of the Ha wings below the continuum ; at AA=2A, f,>0°8. 


TaABLe VI 


Population of the hydrogen u-state in cm=*. The value marked with an asterisk is computed 
for a model optically thin in Le 


T(deg. K) 
N,(em~) 1-0 x 104 1°5 x 104 
20°" 2x10" 4X 10% 
107 gx 10° 3x 10° 
to** 1°5 X10" 4x 10° 


For the flare models of interest here, the magnitude of the bracketed term 
in (13) is found to be only slightly dependent on N, and to have values of about 
2°5 and 5, respectively, for 7 = 104 and 1:5 x ro'deg. KK. Excitation by Lyman 
radiation will result in a similar value for this term in the region around a flare. 
To account for an Hz line of width 4A at the centre of the disk, we thus require, 
at AA=2A, 

7=8x 10-5, (T= 1-0 x 10* deg. K) 
T=4x 107%, (T=1'5 x 10' deg. K), (14) 

‘To check whether these are compatible with the models described by Fig. 1, 
consider first the broadening mechanism. Which of the two most likely causes, 
statistical Stark effect or radiation damping, is dominant depends on the ion 
concentration ; the former holds for N,210!%cm-%, the latter for lower N,. 
The optical thicknesses follow from standard expressions as 


7, =1°§ x 10-2 AA)-3? INNe dz (15) 





for statistical Stark and 
' ry =5°0x 10717(AA)-2| Ny dz 
for natural broadening, AA being in Angstroms. 


The component r" can be found using Table VI which gives, as a function of 
N, and T, values of N, obtained by the solution of the equilibrium equations 





a 
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for the 11-state population of hydrogen. Table VII shows the resultant 7” for 
two values of the model thickness and various electron temperatures and con- 
centrations. "Those models not compatible with Fig. 1 are indicated by an 
asterisk. 

The other component, 7’, of 7 can be found as follows. Neglecting the 
Lyman emission of the underlying material and treating it simply as a purely 
absorbing and scattering medium, it is readily shown that the total intensity J’ 
at the centre of Lx, say, below an optically thick flare is given by 

J' =(27€/KA) exp [ — /(3A)7'], 
where 7’ is optical depth in the underlying medium and it is supposed that the 
scattering parameter A has the same value as in the flare. ‘The population of the 
11-state in the lower region follows from the usual equilibrium equation, which 
in the present case is simply 
AN, =AgN3 

where 4,,N, is the rate of absorption of Lx radiation directed down from the 
flare. ‘The number of 11-state atoms in the line of sight is then given by 


| Ndz= Aw, dz/Azy. (16) 
Supposing the La line to be flat, of width gAv,, Av, being the Doppler 
frequency width, {A,,N, dz is then given by 
gAv ,(27€/KA) 


| 412M, ds = ino GN) (17) 


Tasce VII 


Values of the component t® (AA=2A) for various flare conditions. The values marked with 
an asterisk correspond to models incompatible with Fig. 1 


N, (cm~*) 
T (deg K) z=10'’cm z=108cm 
10! 108 10! 10!% 
1°0 x 104 5X 207° 4x 107% x s0°* 4X 107** 
1°5 » 104 4x 1075 1X 107* 41074 1x 10°** 


For flare temperatures of 104 and 1-5 x 10' deg. K, (17) has the approximate 
values 2 x 10!%¢ and 2 x 107g respectively, and these are only slightly dependent 
on the electron concentration of the flare. Taking g=5, corresponding to a 
flare of optical thickness 10% at the centre of Lx, we find, for T= 104 deg. K, 


| Ne dz=2 x 10!! cm? (18) 
and for 7’=1°5 x 10'deg. K, 


|Nedz=2 x 10% cm, 


The resultant values of 7 are given in Table VIII for natural and statistical 
Stark broadening. The electron concentration N’, in this table is that of the 
region below the flare and has a value given by the equilibrium equation 

2 ee, * 
N, ~ P, +P,’ (19) 
the P's being transition rates to and from the continuum. ‘The excitation rate 
P,, is readily evaluated, knowing the Le intensity, and is approximately 10°! N,” 
or 10°°N,” for flare temperatures of 104 or 1°5 x 10! deg. K, respectively NY 
being the flare electron concentration. The spontaneous rate P.,+ P ‘is 
s a el c2 
approximately 4x 10-"N,.’._ At 500km above the photosphere, N,~10!4cm=3 
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(23) while, at the photosphere, Ny~10!*cm-*. If the Lc radiation penetrates 
to these levels, it follows from (19) that the ionization in the underlying region 
will be almost complete for a flare in which 10¢<7<1-5 x 10'deg. K and 
102 <N, <10% cm’, 

From ‘Tables VII and VIII and the above discussion, it appears that in 
general for flare temperatures ~1-5 x 10'deg. KK, we should expect 7+ >7’, 
while at lower temperatures r*>7°. In either case, it seems possible to find a 
flare model which satisfies the requirement (14) and is compatible with Fig. 1. 
However, consideration of the observed relationship between Ha width and 
central intensity and the variation of this width with time indicates that the 7° 
component predominates at least in the initial stages of large flares, for the following 
reasons. At temperatures >10'deg. K, hydrogen is largely ionized and, as 
shown in ‘Table VI, a decrease in temperature increases N, in the flare and so 
increases rt”. While the bracketed term in (13) decreases with decreasing 
temperature, the net result is that, if 7” were the principal component of 7, an 
increase—or at best a very slow decrease—of line width would accompany a 
reduction in temperature and so in central intensity. This is in disagreement 
with observation, and it is especially difficult to see how the observed rapid 
initial decrease in line width from the maximum value could be accounted for by 
Ha absorption and emission processes in the flare itself. 


Tasce VIII 
Values of the component 7°(AA=2A) for two flare temperatures 


T (deg. K) Tyo Tc 
170 x 10 2x 107° 6x10 "N, 
1°5 x 104 axs2or* 6x 107!"N,’ 


The extension by Lyman radiation of the flare excitation conditions into the 
surrounding material can, however, explain these emission characteristics in 
the Ha wings. Asa flare cools from an electron temperature of 1-5 x 104 deg. K, 
the Lyman radiation intensities decrease rapidly and so, according to these 
arguments, does the Ha line width. As the cooling progresses, excitation in the 
surroundings becomes too small to be of significance, and the H« width is then 
controlled by excitation in the flare itself. It appears, then, that at the maximum 
phase of a large flare, the Ha line is Stark-broadened and, as the cooling progresses, 
natural broadening takes over since, from Fig. 1, an electron concentration 
<10%¢cm- is suggested. In this later phase, where the line width is controlled 
by Ha emission and absorption processes within the flaring region itself, one 
would expect only slow changes of line width with changing Ha intensity. 

(c) The Lyman spectrum and ionospheric disturbances.—The ionospheric 
effects of solar flares have often been ascribed to enhancement of emission in the 
Lyman series and continuum and it is thus of interest to compare the computed 
emission with estimates of the normal solar background in these frequencies. 

The total intensity J , of the emergent radiation from a scattering and emitting 
medium of high opacity is given by the approximate relation : 

__2v(X3)_— ne 
+. say a (20) 
For La and Lc, (20) may be readily evaluated by the methods indicated in (15). 
Resultant values of J, at the centre of L« and at the head of the Lyman continuum 
are given in Tables IX and X. 
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From an assessment of results of rocket observations, Watanabe et al. (24) 
have concluded that the La flux at the earth is ~2 x 10-1 ergscm™ sec"'. Byram 
et al. (25) have given more recent values which indicate a considerable temporal 
variation of this flux, figures as high as 9 ergs cm™ sec~' having been recorded 
in recent rocket flights. Taking the earlier figures as a lower limit and adopting 
a line width of o-2A (26), the equivalent black body temperature of the sun at Lx 
is found to be 6-6 x 10% deg. K and J~2 x 10°* c.g.s. units. For Le no observa- 
tions have been made, but for an effective black body temperature of 
6-0 x 10° deg. K, J~1 x 107"! c.g.s. units. 

Unpublished results of Gardner (27) on ionospheric observations at times 
of large solar flares show that the electron density at about 70 km can increase by 
a factor ~20. Since the equilibrium value of the electron density varies as the 
square root of the intensity of the ionizing radiations, an increase of ~4 x 10? in 
the incident intensity is required if the same radiation is involved in causing 
normal and disturbed ionization. Since a large flare has area ~10-* of the 
visible hemisphere, the intensity of the incident radiation must then increase 
by a factor ~4 x 10° at the flare. Comparison of Table IX with the background 
La flux shows that such a large increase is not obtainable with models compatible 
with Fig. 1 unless the line width increases to an extreme value. For Lc, results 
cannot be definite since the background is unknown ; it appears from Table X 
that, if this background corresponds to that from a black body at 6 x 10° deg. K, 
the Lc flux from a flare may be some 10° times the normal background. 


TaBLe IX 
Emitted intensity in c.g.s. units at the centre of Lx for optically thick layers 
T(deg. K) 
N,(cm-*) 10x 104 1°5 x 104 
10!” I'2* 1077 5°9 X 1078 
10} 3°4 x 1077 I‘9X 1075 
TABLE X 
Emitted intensity in c.g.s. units at the head of the Lyman continuum for optically thick layers 
T(deg. K) 
N,(cm-*) 170X104 I°5 X 104 
10! 3*gx10°° 9°7 X 1077 
10” 8-7 x10"° 7°8 x 107° 


It seems from the above that we cannot ascribe both normal and enhanced 
D layer ionization to the absorption of L« nor by inference to any other of the 
Lyman lines. ‘There is insufficient observational evidence to make any similar 
conclusion for Lc radiation. 

A possible alternative source of ionizing radiation might be the helium 
resonance lines and continua which, because of their higher excitation potentials, 
will be enhanced to an even greater degree than the corresponding hydrogen 
emission. There is, however, as shown by C. Warwick (28), an additional 
condition which the ionizing radiation must satisfy. Her analysis of the associa- 
tion between heights of limb flares and the occurrence of SID’s suggested that 
the ionizing radiation responsible for the production of an SID is absorbed in 
tangential passage through the solar atmosphere unless formed at a height above 
about 15 cookm. From an investigation of the opacity of the chromosphere in 
various possible ionizing radiations—X-rays, Le, LB, and La—Warwick con- 
cluded from the above criterion that Lx is the most probable cause of the SID’s. 
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Since, however, the first helium resonance line at A 584 has an absorption coefficient 
almost equal to that of Lx and since, also, the number of ground state He atoms 
at 15 oo0km above the solar surface is at least as large as for hydrogen, the 
tangential opacity in A 584 is at least as large as that of Lx. Warwick’s observa- 
tions, therefore, could be interpreted equally well in terms of the helium 
resonance line as the SID producing agent. 

There are, however, difficulties associated with this interpretation. — Firstly, 
the actual flux in the helium resonance lines is low even for a flare at 15 000 deg. K, 
although the flux increases rapidly with increasing temperature and may be 
sufficient to produce SID’s with a large flare. Furthermore, published data 
on absorption cross-sections and atmospheric concentrations suggest that these 
resonance lines will be absorbed at ~180 km above the Earth’s surface. The 
usual statement that the radiation cannot penetrate to the 70 km region is, however, 
not necessarily correct; the radiation is not all lost on absorption since, as the 
atom density is low, a good deal will be re-radiated by spontaneous recombination 
so that the radiation is effectively scattered and may penetrate to much lower 
levels than that at which it was first absorbed. However, excess ionization 
would be produced from the 180km level down, and the fact that this is not 
observed argues against attributing SID’s to radiation, such as the helium 
resonance lines, which is absorbed first at high levels in the atmosphere. 

5. Conclusions.—Generally well verified observational results on the hydrogen 
and helium emissions of class 2-3 flares have been shown to be consistent with a 
uniform and static flare model whose electron temperatures and concentrations 
lie respectively in the ranges 1-0 x 10'to 1-5 x 10' deg. Kand5 x 10! to~10% cm’, 

It has been shown that excitation in lower regions due to Lyman radiation 
from the flare is adequate to explain the Hx line width for a flare temperature 
~1+5 x 10' deg. K, the line being Stark-broadened at these levels. On this basis 
the observed relationship between Ha central intensity and the line width may 
be accounted for, and the observed rapid decrease of line width away from the 
maximum receives a ready explanation. It has also been shown that the La 
flux flare is insufficient to account for the increased ionization observed simul- 
taneously in the D layer, and that there are difficulties in attributing it to the 
helium resonance line A584. 

For a complete analysis along the lines indicated in this work far more data 
are required on flare spectra. In particular, simultaneous measurements of 
the helium and hydrogen spectra for both limb and disk flares are urgently needed 
and observations in the region of the Balmer continuum could furnish very useful 
data for limiting the possible range of physical conditions in flares. 
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Notes added in Proof 

1. The referee has drawn the author’s attention to two interesting recent 
publications by Svestka* on the physical conditions in flares ; the results of these 
are generally in agreement with those found above. It is interesting to note 
that he also has been forced to the conclusion that Lyman « radiation into the layers 
below the flaring region must be considered to explain the Ha width. He does 
not, however, consider the influence of Lyman continuum in ionizing this material 
and so producing Stark broadening. Further if the Hx opacity arises largely 
from a region below the flare, then Svestka’s analysis based on the assumption 
that the 7 = 10 state atoms—whose total number is derived from the Ha width—are 
distributed uniformly through the whole flare cannot be valid. Although, in 
general, Svestka recognizes the importance of departures from thermodynamic 
equilibrium, he implies in places an equality between kinetic and excitation 
temperatures. ‘This is not acceptable and the run of temperature with depth 
derived this way is invalid since the excitation temperature must have an intrinsic 
increase with depth due to the increase in radiation intensity into the flare. 

2. A recently completed investigation by the author on the influence of non- 
coherent scattering on the shape of emission lines indicates that the emergent 
Lyman «@ line from a flare will be strongly self reversed and that, while the central 
intensity given here should be correct, in the wings the intensity may increase by a 
factor ~10*. Even so it seems difficult to account for the observed effects in the 
D region in terms of excitation by Lyman «. 


*Z. Svestka, Czech. Acad. Sci. Astr. Inst., Pub. No. 32, 1957; and B.A.1. Czech., I, 130, 1956. 
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POSITIONS OF COMET AREND-ROLAND (1956h) 


R. H. Garstang 
(Communicated by the Director of the University of London Observatory) 


(Received 1957 June 22) 


: The following positions of Comet Arend—Roland (1956 h) have been obtained 
: with the Radcliffe 24-inch photographic refractor at the University of London 
Observatory. 
‘ 1957 Ue i 1950°0 Parallax Factor 
3 x ; 8 X é5 
4 
t hm ss gy s 
: 1 Apr. 24°87249 2 33 32°26 43 12 12°4 0°394 +7°31 
$ 2 Apr. 27°89564 3 11 26°74 51 9 28°5 0°443 77°04 
2 3 Apr. 27°91900 3.11 45°11 +51 12 30°2 +-0°387 -+-7°47 
i 4 Apr. 28:87160 3 24 16°63 +53 9 35°1 -0°532 + 6:24 
j 5 Apr. 29°89831 3 37 56°04 54 59 39°1 0°513 + 6°59 
‘ 6 Apr. 30°91324 3 51 28°65 +56 33 25°2 0°513 +6°67 
: 7 May 13°89651 6 16 26°12 +63 37 22°3 -o'809 + 3°85 
i 8 May 13°90794 6 16 31°79 63 37 23°5 -0°797 + 4°22 
9 May 22°92768 7 16 40°00 +63 16 41-0 0°785 +4°34 
10 May 22°94846 7 16 46°67 63 16 34°4 +0°756 +-4°95 


The comparison star places were taken from the Zweiter Katalog der Astro- 
nomischen Gesellschaft (AGK 2), proper motions not being included. 


Comparison Stars 


I 43 296, t+ 43 299, 42 271 
2,3 51° 303, 51 308, r50 351 
4 53° 316, FS3 320, 52° 349 
5 55° 352, +55 358, +54° 372 
6 +57 436, +56° 437, +50° 434 
9% 64° 346, +63° 410, +63° 409 
9, 10 63° 445, 63° 450, +63" 448 


I am indebted to Mr T. Kiang for making his plate of April 24 available 
to me for measurement, and to Lieut Cdr L. M. Dougherty for voluntary 
assistance on several evenings. 


University of London Observatory, 
Mill Hill Park, 
London, N.W.7: 
1957 June 17. 
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IMPERFECTIONS OF ELASTICITY IN 
THE SMALL BODIES OF THE SOLAR SYSTEM 


Harold feffreys 
(Received 1957 June 4) 


Summary 


‘The rough estimate of the damping of the 14-monthly variation of latitude 
gives two types of solution on the hypothesis of elastic afterworking, in which 
the times needed to attain the strain for long-continued stress are of the 
orders of a few weeks and nearly a year. ‘These are applied to tidal friction 
in the Moon and Mercury. The long scale leads to quantitative difficulties. 
For the Moon the straightforward application, in which it is supposed that 
the imperfection of elasticity is elastic afterworking and applies to the whole 
body, gives satisfactory results, but if it applies only to the inner part it is 
apparently necessary that the rotation and revolution periods were brought 
into close coincidence when the Moon was much nearer the Earth than it is 
now. For Mercury the effect is too small to have produced the required 
result in the time available, and a greater degree of imperfection of elasticity 
seems to be required, at least for periods of the order of a day or at some 
earlier stage of its history. The shorter time scale for the relaxation leads to 
no serious difficulties. 

Comparison with satellites of the outer planets shows some peculiar 
features, but these may be attributed to differences of composition. 





1. It is well known that many bodies in the Solar System keep a constant face 
towards their primaries, and that tidal friction provides a qualitative explanation ; 
further, no other explanation has been suggested. The work of Darwin was 
mostly concerned with tidal friction in the body of the Earth, treated as a highly 
viscous liquid; but there was no independent evidence of the value of the 
viscosity. Later work has shown that bodily imperfections of elasticity in the 
Earth cannot account for the apparent secular accelerations of the Moon and Sun, 
which must be attributed to turbulence in seas with fast currents. Most of the 
bodies concerned have no such seas. Positive evidence for any imperfection of 
elasticity of the Earth under small stresses has been hard to find, but some has 
at last appeared from the damping of the 14-monthly variation of latitude. It is 
therefore desirable to examine whether this degree of imperfection of elasticity 
is adequate to account for the rotations of the Moon and Mercury and possibly 
other satellites. In all cases, besides the steady rotation and revolution, there 
are theoretically forced and free vibrations. Some of the forced ones are per- 
ceptible for the Moon ; the free ones are too small to have been detected. For 
other bodies even rather large librations could not be detected. A further 
problem for the Moon is to explain how the free librations, probably originally 
large, could have been damped enough to be imperceptible. 

In an early paper (Jeffreys 1915), I discussed changes of figure of the Earth 
and Moon under the influence of changes of rotation due to tidal friction, and 
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the possibility that the Moon’s excess ellipticities are the result of its having 
solidified when much nearer the Earth than it is now. I found that even if the 
Moon has behaved as if perfectly elastic the smallness of the free libration in 
longitude is consistent with this hypothesis, provided that the periods of rotation 
and orbital motion had been brought into coincidence before solidification. 
Even if there was considerable libration at solidification it has had ample time 
to die down if the Moon has the viscosity of pitch, which was taken as of order 
10'°c.g.s., but there were difficulties about the persistence of the ellipticities 
on the hypothesis of finite viscosity, and in any case this viscosity is far too small. 
2. Libration.—In the first place I consider the elastic case, since the equations 
for this case can be easily adapted to various types of linear imperfection of 
elasticity. ‘The analysis is similar to that of Jeffreys and Vicente (19572, b, 
referred to as I and II) for the Earth’s nutations. The displacements were 
expressed by small rotations as for a rigid body, with an elastic set superposed 
that makes no contribution to the angular momentum. ‘The elastic displacements 
correspond to co-ordinates that would have short periods in a free vibration, 
and can be eliminated at an early stage. Additional displacements had to be 
introduced to represent the motion of a liquid core, but these need not be con- 
sidered here. ‘The main other differences are as follows. (1) A full treatment 
of the elastic displacements of a planet must take account of initial stress. Love 
(1911) introduced the hypothesis that the initial stress is hydrostatic pressure, 
and we adopted this for the Earth, for which it is a good approximation. The 
librations of the Moon, however, depend on the ellipticities, which far exceed 
the hydrostatic values, and it is desirable to examine directly whether the associated 
initial stress-differences have any importance. (2) The usual theory treats the 
Moon as a rigid body, and the initial stresses do not matter. But elasticity 
must be considered before we can consider damping by imperfections of elasticity. 
There are theoretically three free periods of libration, none of which has been 
established by observation. ‘Their ratios are of the order of 60, and for some 
types of imperfection of elasticity they might be very differently damped. ‘The 
fact that they seem to have disappeared, combined with the persistence of the 
ellipticities themselves, should give useful inequalities for the imperfections in 
question. (3) For the Earth a term representing the changes of obliquity has 
to be added to the gravitational potential in determining the elastic displacements. 
For the Moon there is another term due to the potential produced by the Earth. 
(4) For the Earth the angular velocity about the axis of greatest moment can be 
treated as constant, but for the Moon extra displacements in longitude are an 
essential part of the problem, and a new co-ordinate is needed to represent them. 
2.1. Let the initial stress be —p 3;,+49),, with qg,=0. ‘The undisturbed 
geopotential is ¥. The equations of equilibrium are 
ov Po 4 OFik =, (1) 


POdx, ax, | ax, 





The surface stress in the steady state is 
—lipot lidins (2) 


and is zero since the surface is free. ‘Taking the work done by the stresses q,,. 


to the first order, we must replace W,, of I by 
* 


3 


wv 
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W,= || pol ta + May * — Ha a) 


- {Ufot +s, ‘am (a) 
i [fam Ga. (3) 


The last term is 
— | { dgunedS + {| mx Bar (4) 


The surface integrals cancel to the first mms by (2) and the volume integral in W, 
arising from ¥ is, to the first order, 


y 
[ff pow = dr. 


The contributions from pp, gj, and Y cancel by the equations of equilibrium. 
Thus, as we should expect, initial stress introduces no new first-order terms. 
The second order terms are of the form 4j7.€n€ym3 but the elasticity terms are of 
the same form with the elastic moduli replacing g;,. Hence the second order 
terms arising from initial stress-differences would only be of the order of 1/1,000 
of the elastic terms even if the material was on the verge of fracture. ‘They are 
therefore negligible and the previous form of W,+ W, holds. 

2.2. As for the Earth we take the axis of x, in ‘the mean direction of the minor 
axis and the other axes rotating about it with constant angular velocity n. Rota- 
tions /, m of x, towards x, and x, are introduced ; we take also an extra displacement 
x about x;. A particle at (x,, x, x3) goes to (€,, &, £3), where, by 2.2 of I, 

f= xy(1— HP by?) tu) + Uy +) — (x + $lm)(x9 +05), 
Eg= x9(1 — 4m? — 3?) + uy + m(x3 + us) —(x— $lm)(x, +4), 
Ey = %y(1 — 3? — bmn*) + us — (1+ mx)(x, +0) — (m —Ly)(%2 + U9). (5) 

2.3. If the standard position of the Earth is on the axis of x, the extra terms 
in ¥ are 

2 


1p 2( 2 4. x2 ad a ae 
jn (m+ 93)+ Sy %3— x3), (6) 


where « is the ratio of the masses. Their contribution to W is 


$n*[1?2(C — A) + mC a B)- | J pol; = {x3(lx, + mx,)}dr7] 


— rp SMC-A)+3x(B-A)~6 {|| pou (llrs—x)}de] (7) 


Since we shall consider only free vibrations we need not consider displacements 
of the Earth from the axis of x. 
Apart from a perfect differential, 


n | | | PolErés—€,£,)dr = 4n(B + A —C)(Im—im) 


re . : = 
+n | J J Pott; oa, ( — Lens + mx 1X3 + 2X%,X,)d7, (8) 


||| po Gi+ E+E )dr= (BP + Ant + CZ), (9) 
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We have used the conditions that u; are chosen to give no angular momentum and 
that the associated speeds of free vibration are high. Their contributions to W 
are equivalent to introducing an extra ee potential 

U) =n?x,(lx, + mx,) + G a ae (tes XX 1X2) 


+n(— lxgsxs + MX Xs + 2XX1%Xe). (10) 
and the Lagrangian function is 


= (BP + Ari + Ci) +In(B + A— C)(lin — Im) — }n®{(C — A) + mC — B)} 


~ re GC-A)+3x(B-A)} 


+| {| eau (SS 1 2 \ar—| {I (AA8 e+ ann) Se dr 


+ [ feowi( su = HF LS sat) dr (11) 


Ox,OX;, ~ OX; OX; ~ Ox, Ox 








u, are quasi-statical co-ordinates, and if a statical solution is used for them all 
quadratic terms in them disappear and the linear term (namely that in U/) is 
halved. Then the terms in wu’, can be replaced by 


Ayr 


pia al - eet ap a) . 
S| | | pou; adr = 3 | | polje;U dS — } || | U; 5, (Potti)dr- (12) 


It is enough to treat the Moon as homogeneous and incompressible ; then the 
volume integral vanishes. For U| =k,K>, with K, any solid harmonic of degree 2, 
the normal displacement is 


lu;=€K,/a, (13) 
where 
_ _ 5Rapa® 
~ Tou -+ 2gpa (14) 
In the case 
kyK, = jn*r*(} — cos), (15) 
corresponding to the effect of a steady rotation, we get 
— 
“~ 2(19n+2gpa) ” ie 


say ; v is the ellipticity, here equal to the dynamical ellipticity, produced by the 
rotation. Then the terms in L arising from the deformation reduce to 


tvC | B+ nit + git (2+ +) n(Im — Im) 











2 Be ny, RE 
+n {(1+ +.) P +m? + +a (17) 
The terms in /, m, x are equivalent to increasing A, B, C by 
2vC ) 
5A =4vC+ oo 
vC 
dB=hvC— — , (18) 
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The changes of the differences of the moments of inertia due to the potential ¥ 
are just those given by these expressions with the sign reversed. If we write 
A+68A=A, and so on, 


v ae a. ae vC \.,) 
L= + {(By+ = = + an Pe +(4y+ Se. — ~~) +(¢ or 3" so Peres Jee} 
8 Ae 
+ nf By + Ay—Co+Co(- + 4) (in im). 


— £ n{Cy— Ag) + m*(Cy— By) }— 3m ry Co Ay) + x7( By — Ag) 5 (19) 


2(1+K 
We write 
Co— Bo= Age 3 Co—Ap= BoB ; Bo— Ao = Coy: (20) 


The free speeds, where squares of the ellipticities and their products by v are 
neglected, are 


(1+ erat ae): m(a+; ep Gk) eo 


The two slower speeds are to this accuracy the same as for a rigid body with 
moments of inertia Ay, By, Cy. The speed near n is slightly increased. As an 
indication of the probable amount, and for further comparisons, we notice that 
B is about 0-0006, « about 00004. If the Moon had its hydrostatic oe gaaras 
B would be about 0:00004, « about 0-000009 ; and »* is about 1-5 x 10°°. ‘The 
term in v for the speed near n is therefore of the order of 10~* of the 8 term. The 
uncertainty of 8 as found from observations of the inclination of the axis to the 
ecliptic is about 1 part in 400. The point may possibly become important if the 
accuracy should be increased. A speed near n represents a slow motion in space ; 
and £ is determined from a 19-year period, so that a small change in the resonating 
frequency might matter. The point is analogous to the difference between the 
effects of elasticity of the Earth on precession and the variation of latitude. 

3. Correction for imperfect elasticity.—T he simplest type of elastic afterworking 
is the linear form specified by replacing the rigidity » by the operator (Jeffreys 
1952, p. 10) 








1+1/pr’ 
T+ i/pr’ “ 


where for a time factor expiot, p will be replaced by ic. For astress P applied 
at t=o and kept constant, there is an initial elastic strain P/u), which increases 
asymptotically to (P/u»)(7'/7) if the stress is maintained for a time longer than 7’. 
Since for the Moon gravity has a small effect compared with rigidity we can take v 
to vary as 1/u. Since we have the solution for constant v, say vp, we can adapt 
by putting 


v=Vy+V, (2) 
7 —T 
"= Ye 4a)’ (3) 


adding terms in v, to the Lagrangian, and solving the modified period equation 
by successive approximation. As the data for 7 and +’ are very rough in any 


* (u/p)*/* is the velocity of transverse elastic waves; in the greater part of the Earth’s shell it 
ranges from about 4°4 to 7 km/sec. I generally take 4-5 km/sec. 
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case we may as well drop the small number «. The time factors for the free 
modes are respectively 





Ca 3 ay _ , M¥o(7 —7)7’ 

tin(1+ > B) 2 (nr? +1)’ (4) 

ia -_— re. vor —T)r (5 

+ in (4B) 4n*ap ( + =) 7(4n®aBr? +1)" (5) 
2 , =a , 

2 in (gy)! — 2? (6) 


2 1(3yn?7’"? +1) 
‘The case 7’ = co is that of elasticoviscosity. 


The only phenomenon in the Earth’s shell that seems to give positive evidence 
for imperfections of elasticity at small stresses is the damping of the free variation 
of latitude. There are two types of solution according as or, o7’ are large or 
small for this speed. If 7, 7’ are large (say > 60 days) the solution should satisfy 
(Jeffreys 1952, p. 247, (21). 10-* should be 1o~*.) 

Ee I I 


6 7 7’ 2x10%sec. 6years 


(7) 


This may be wrong by a factor of 2. If the term 1 in the denominators of the 
damping factors is dropped (corresponding to infinite 7’), these factors become 


simply 
2% 2 1 \ % 3 Yo 
-9-Geae (8) 


With the numerical values already adopted the times of relaxation are, roughly, 





2x 10° years ; 10 years ; 3 x 10° years. (9) 
If we take r’=1-17, we get 7’ =0°6 year. The free periods are roughly 1 month, 
100 years, and 3 years. Then the first relaxation time is hardly altered ; the 
others are multiplied by about 600 and by 3. 
The alternative solution for o7, 7’ small gives 
7’ —7r=1°8 days (10) 
In this case the relaxation times are 10°, 4 x 10’, and 4.x 10° years. The ratio r’/7 
does not affect the results appreciably. 

On the facg of it the damping in any case is enough to account for the smallness 
of the free libMtions if the time since the Moon became solid is comparable with 
the age of the Earth. 

3.1. Subsidence of a surface inequality.—There are however two possible 
difficulties. If the Moon is imperfectly elastic, the ellipticities themselves must 
be tending to subside and leave the Moon in a hydrostatic state. ‘The differential 
equation satisfied in a free body for a surface harmonic of degree n is 


+ ye bem, 
(2n* + 4n + 3)m 


where , must be replaced by an operator as in (1). For a second harmonic on 
the Moon this leads to a damping factor 


I O15 
A=[- ~ I‘O15 I 
(3 + 225) /1015 (:1) 


With 7’ = », r=6 years, this would give a time of relaxation of about 400 years ; 
the inclination of the Moon’s axis to the ecliptic would have changed considerably 
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even during the interval of accurate observation. Elastic afterworking, with 
the same value of 6, would make the time even shorter. It looks as if the excess 
ellipticities could not have persisted long enough for the librations to become 
appreciably damped. 

The solution appears to be that the stress-strain relation that leads to (1) 
assumes that the body is not self-stressed. If there are stress-differences in a 
state of equilibrium, there is no reason why there should not be a linear relation 
between stress and strain when both are measured from this equilibrium state. 
Then the times just found would represent the time-scale of the return to this state. 
There is even no conclusive reason against using the elasticoviscous relation in 
this way. Onthe other hand if we made the usual assumption that change of shape 
and stress-difference are always measured from the hydrostatic state, and adopted 
the elasticoviscous law with the viscosity derived from the Fenno-Scandian 
uplift, the time of relaxation would only be increased to something of the order of 
40000 years, and it would still be impossible to suppose that the ellipticities have 
existed for anything like the age of the Earth. 

3.2. The second difficulty is as follows. The evidence for the Earth points 
to a strength comparable with that of surface rocks to a depth of order 600km. 
The viscosity inferred from the variation of latitude may refer only to the lower 
part of the mantle. If the difference is due to temperature, we should expect 
the transition to be at about the same depth in the Moon; but this depth is 
nearly half the radius. Couples due to tidal friction in a uniform body are 
proportional to a high power of the radius, and it is possible on this ground that 
the times should all be multiplied by, say, 2®. ‘This would bring the longer ones 
to the neighbourhood of 10° years. ‘The initial amplitudes of the librations are 
of course unknown, but it is reasonable to suppose that those in latitude were once 
of the order of 5° and that in longitude near go° ; and a factor of e~*® would not 
go far towards reducing these values to their present values, which can hardly 
be more than 1’ or so. This is not fatal, since the factor 2° is dubious and the 
viscosity may be wrong by a factor of 2 ; if the damping has been by e~® it would 
probably be sufficient. 

However, the damping factor contains a factor n*. If, as I originally supposed, 
the Moon solidified when much nearer the Earth than it is now, nm? was much 
larger. Also the cooling might not have reached so great a depth. On both 
accounts the damping may have formerly been much more rapid than it is now, 
and the difficulty disappears. 

4. The theory is not difficult to adapt to a freely rotating satellite. Only 
the effects on / and m arise, since the amplitude and phase of y are replaced in 
free rotation by the third Euler angle and its rate of increase ; a forced motion is 
superposed. ‘The results remain of the same order of magnitude so long as the 
mean motion and the rate of rotation are comparable. The secular effects con- 
sidered by Darwin are of course of the second order in the tidal potential. 

5. Tidal friction and rotations of satellites—In ‘The Earth’, 8-02 (13), 
I show that for a homogeneous body the tidal frictional couple is 


a’ 
= — <0> pH sin 2«, (1) 


where p and a are the density and radius of the disturbed body, m, c the mass 
and mean distance of the disturbing body, H and 2 the maximum amplitude 
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and lag of the semidiurnal tide. For a satellite with rotation speed w and mean 
motion n 
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Spa®—3,fma 


2 
yi , a 
ion +anee” at sin? 6 exp{2i(y+wt—nt)} (2) 


H sin? 6 exp{2i(x + wt —nt)}= 





in which we must write 
1+1/2t(w—n)r’ | I 
BVO T+ 1/2i(w—n)r =H ~ 2i(w— nyo ) (3) 
provided (w—n)r, (w—n)r’ are large. For the Earth the gravity term is rather 
small compared with the rigidity term where they are combined, and for smaller 


bodies it is much smaller. Then effectively H has the value for a perfectly 
elastic body and for large 7, 7’ 





I 








2€= 2(w — nyo . (4) 
We write as before 
fm/3 =n? /(1 +x)=n’, (5) 
and in all the angular acceleration is 
; N go7 / n’a \2 p I 
eee = (7) 19fy + 2gpa 2(w—n)b- (6) 


As trial values for the Moon we take a rotation period of 6 days and a revolution 
period of 12 days. We then find 

w=2 x 10°!*/sec?, (7) 
and a change of the order of the initial value would need something of the order 
of 2x 10° years. At the present distance this time would need to be multiplied 
by about 20, and if the imperfection of elasticity is confined to the inner half of 
the radius a further factor of the order of 100 might be needed. It is therefore 
very unlikely that the Moon’s present rotation could have been established 
unless it was formerly much nearer the Earth than it is now. 

We may notice that the differential equation for w is of the form 

w= —A/(w—n), (8) 
which has the solution 
(w —n)? = (wy —n)*? — 2At. (9) 
With the above values w would become equal tom in about 1o*years. ‘The approxi- 
mations would however fail when w—n is small enough for € to become a large 
angle, and the approach would become asymptotic ; and at a later stage still the 
libration case would be reached. 

For Mercury, if we take the same initial rotation period, we get a time near 10° 
years. With a larger initial rate of rotation, as is probable, this would be much 
prolonged. It appears that there is a serious difficulty in accounting for the 
present rotation of Mercury even if the imperfection of elasticity adopted for the 
Earth is supposed to hold right up to the surface. 

In the case (w—n)7, (w—n)r’ small, the lag 

2e=tan-1{2(w —n)(r’ —7)}=1 (10) 
with the values from 3 (10); whereas (4) makes it 1/2000. Thus in this case the 
value of w just given would be multiplied by about 2000. ‘Thus the smaller values 
of 7, 7’ afford adequate tidal friction to account for the rotations of the Moon and 
Mercury without further adjustment. 
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6. Other satellites. —Far the most important factor in w for satellites is (na)?, 
apart from hitherto unknown variationsin#@. Itis therefore of interest to tabulate 
rough values ; I give (diameter/period*)~+(km/day*). For most of the dia- 
meters of Saturn’s satellites I use results communicated to me in a letter from 
Professor Kuiper; for the rest I take those given by Russell, Dugan and Stewart. 


Moon 5 Mimas 560 Ariel 8o 
Phobos 160 Enceladus 270 Umbriel 22 
Deimos 6 Tethys 230 Titania 13 
JV 600 Dione 110 Oberon 5 
JI 1200 Rhea 8o 

Jil 300 Titan 18 ‘Triton 130 
jill 100 Hyperion o9 

JIV 18 lapetus o-2 

Jvi 0°002 Phoebe o'0o1 


The great majority of satellites’ rotations should be more sensitive to tidal 
friction than that of the Moon. It is therefore not surprising that all whose 
rotations are known keep a constant face to their primaries except those of Uranus 
(with a possibility of a very small difference for JIII towhich Dr W. H. Steavenson 
‘has called my attention). The best-authenticated case of all, however, is lapetus, 
for which the effect should be much less than for the Moon. 

Dr Steavenson informs me that Titania and Oberon have axes not normal 
to the orbital planes and possibly in them. ‘There are cases where this state is 
theoretically stable, as Stratton (1906) has shown, but the physical properties 
required are very different from those that appear to belong to the Earth. There 
would be nothing surprising in this, since the small densities of some of the 
satellites of the outer planets would already suggest very different constitutions, 
possibly largely ice at very low temperatures. But it seems worth while to state 
the result that the rotations of Iapetus and the ‘satellites of Uranus provide 
additional evidence. 

7. The law of elasticoviscosity, proposed by Maxwell and occasionally used 
by Darwin, was originally based on a model in which, under constant stress, 
the molecules break down at a rate proportional to the stress and permit flow 
at a steady rate superposed on the initial elastic yield. In elastic afterworking a 
steady stress beginning at a certain time produces an initial elastic yield, but the 
strain increases asymptotically to a new constant value. This behaviour is 
represented in the simplest possible way by the law 3(1). Qualitatively it has 
experimental support for some rocks, but in laboratory conditions the stresses 
have to be rather large and the response is not linear. In the present problems 
all the stress-differences that concern us are so small that we are entitled to assume 
a linear law. It has been shown that if a perfectly elastic material encloses 
pockets of an elasticoviscous one, the response as a whole can follow the law 
suggested. . 

A perfectly elastic body can be in equilibrium under internal stress without 
external force ; a rubber ring gripping a rod is a simple example. It would 
be quite possible to have a perfectly elastic solid with such internal stress 
initially, and containing elasticoviscous pockets. In such a body the law for 
elastic afterworking would hold, the relation being now between departures of 
stress and strain from their initial values. The elasticoviscous connection between 


them could then be regarded as the formal limit when 7’ is made infinite. I know 
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of no suitable model for this case, but there seems to be no harm in considering 
the possibility. 

8. The results for the Moon and Mercury with the larger values of 7, 7 
are none too satisfactory. Application of the imperfection of elasticity derived 
from the variation of latitude to the Moon in the most straightforward way 
indicates that the rotation could easily have been brought into coincidence with 
the revolution and the free librations damped out. But, if, as is probable, it is 
effective for only the inner half of the Moon’s radius, the damping would be 
much slower and it is not obvious whether the time available would be enough. 
If the Moon was formerly very much nearer the Earth the difficulty disappears. 

For Mercury similar analysis suggests that if the original rotation was com- 
parable with that of the Earth and Mars, and 7, 7’ are long, the time needed to 
reduce its rotation to the present value would be much too long. We cannot 
avoid the difficulty by supposing Mercury to have formerly been much nearer 
the Sun. If it was once a satellite of Venus there is no difficulty in accounting 
for the slow rotation, but the time needed for it to have receded so far that the 
satellite orbit became unstable would probably be over 10!! years. It might 
be necessary to suppose that imperfection of elasticity has at some time been 
greater in Mercury than in the Earth. This would indicate internal temperatures 
higher than at present. On the other hand, if the smaller values of 7, 7’ are 
correct, all these difficulties disappear. 

The experiments of D. W. Phillips suggested r’/7 = 1-6 at rather large stresses, 
the new steady state being reached in some weeks. For small stresses 1-1 is 
likely to be an overestimate, but as it leads to rt about 18 days for the smaller 
solution it is possibly not wildly wrong. Further, it is likely to be applicable 
to cold materials, and there may be no reason to suppose that the dissipation is 
confined to the central regions. The indications, especially for Mercury, are 
that the solution 5(10) leads to fewer difficulties than 5(4) does. 

The satellites of the outer planets are probably largely ice and the rigidity 
is probably much lower than for rocks ; with equal times of relaxation (7, 7’) 
this alone would greatly increase the effects of tidal friction, and the imperfection 
of elasticity may also be greater. If so, the apparent smallness of the theoretical 
effect on Iapetus may well be compensated. 


, 
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FINITE ATMOSPHERES WITH ISOTROPIC SCATTERING 


III. CorRIGENDUM AND ADDENDUM 
I. W. Busbridge 
(Received 1957 July 24) 


Summary 


The conclusions of the Appendix of Paper I are wrong. A rigorous treat- 
ment has established the correctness of all the main formulae and it has led to 
some new and important relations. ‘These are summarized, and they are 
used to re-establish any results depending on the Appendix. 





1. Corrigendum.—Hopf has proved in (3), Section 21, that there is no non- 
null, bounded solution of the homogeneous equation 


Ita) = Jeo) “J (E( |t—z|)dt (o<w<1), (1.1) 


i.e., in the notation of Paper I (see (1)), of the equation 

J(r)=aA,{J(t)}. 
It is not difficult to go further and to prove that, if J(7) is O(/n7~) as r-> +0 and 
O(In[r,—7]}"1) ast >7,—0,thenJ(r)=0. These results contradict the conclusions 
of the Appendix to Paper I. 

In the conservative case the error is due to the non-uniqueness of the solutions 
of the X- and Y-equations, so that the argument cannot be reversed. In the 
non-conservative case it is due to the fact that X(u) and Y() have no poles. 
(See Section 2, below). 

In the main part of Paper I, all the final results are correct. ‘This is even true 
when, as in Sections 8 and 10 (conservative case), they are apparently based on 
the Appendix. In the non-conservative case the solution of (1.1) has been taken 
to be zero and, provided that X(u) and Y() are defined by I(4.21), the analysis 
of these sections can be justified rigorously. In the conservative case the analysis 
has to be altered. After summarizing some new and important formulae (Sec- 
tion 2), the revised analysis is given in Sections 3 and 4. 

2. Addendum.—If there is no non-null solution of the equation (1.1), the 
only solution of the auxiliary equation 


I(r, 0)= aK ,{J(t,0)} + exp(—or) (2.1) 
is the N-solution (obtained by iteration), viz. 


J (1,0) = yori "{exp(—ot)}. (2.2) 


Since 


A,{t}=1-4$E (7) —4E,(7,—7),| 
therefore 


A,{t} <1 — E,(§7)=p <1. (2.3), 
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It follows that, if f(r) is continuous for o<r<7,, so that |f(r)| <K (say) in 
(0,7,), and if o<w<1, thent 


|w"A,"{f(t)}} | <KA,"{1}<Kp”. (2.4) 
Hence > w"A,"{f(t)}is(absolutely) convergent foro <@ <1,0<r<7,. Itfollows 
n=0 


that the series (2.2) converges for o<mw <1, 0<17 <7, and any finite o. 
Let X(u) and Y() be defined by I(4.21), viz. 


X(u)=J(0,u), ¥(u)=J(r 0). (2.5) 
Then X(u) and Y(,) are positive and continuous for 4. 0 and it is easy to show 
that 
X(u)>1, Y(p)->0o as p> +0 (2.6) 
X(4)>J(0,0), Y(u) J (0, 0) as |u| oo. (2.7) 
The equality of these limits as || 0 is most easily seen from the relation 
I (4.26), viz. 
Y(u) = exp(—7)/H)X(—p), (2.8) 
together with the fact that /(7, oc) is a continuous function of o at o=o. 
From (2.2),we have 


oo a)” 2 


nad ~ 2 (—oct)™ +s a a 
Ira)= Sarin) YAU} CAS wr n(enp 
° 0 





n=0 u=o ™ m=o0 ™ n= 
The inversions of the orders of summation and integration can all be justified 
by absolute convergence on using (2.4). Putting s=o0 and 7,, and o=1/p, we 
have the expansions (valid for | |> 0) 


, < ( wet 1)"Cy —m , ~ _- 1)"Dy, —m 
3-2 ape. Te 2 ee (2.9) 
where 
Cc. -| > wo" A ne] : (2.10) 
n=U r=0 
D,,=| > aw" A” | ; (2.11) 
n=0 TT 


Relations between the coefficients C,,,, D,,, the moments of the X- and Y-functions 


(see 1(7.14)) and 7, are found by substituting into I (4.22) and I (4.26), expanding 
in powers of »~! and equating coefficients. It is found, in particular, that 
Cy = Dy =J (0, 0) = 1/(1 — Xp + Yo): (2.12) 

The expansions (2.9), which converge for || >0, show that X(u) and Y(p) 
are analytic functions of the complex variable which are regular for |n|>o 
{infinity included). 

Conservative case (w=1). In this case it can be shown that the moments of 
X(y) and Y(,) satisfy the relations 


Xp+Yo=l, (2.13) 
Xy—-¥1=7T1¥0» (2.14) 
YolX2 + Yo 374(%1 +91)) =¢- (2.15) 


t Since the kernel of A is positive, if g(r) </(r) for o& r< 74, then A+{e(t)} <A, f{h(o}. The 
inequalities in (2.4) follow by repeated application of this result. 
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Since J(o, 0) is finite and not zero, it follows from (2.12) and (2.13) that x91, 
yoo. Thus the functions defined by (2.5) are not the standard solutions. 

By the method indicated above, it can be shown, after some lengthy calcula- 
tions, that 


Cy = Do= 1/20 (2.16) 
Cy=4(1Co—y), Di=HuCot+y), (2.17) 

where 
y= 47,370 + 6(x3 — Vs) + 371(%2 — Ye). (2.18) 


3. Section 8 of Paper 1.—With the X- and Y-functions defined by (2.5), 
the analysis of Section 7 (0o<@<1) of Paper I holds also when w=1 up as far 
as 1(7.18), viz. as far as the equations 

J(0)(1 — X9) + I°(71)¥0 = 40> 
J(0)yq + I%(71)(1 — Xo) = ao. 
These are identical, since x9 += 1, and we have only one equation 


J%0) + I(r) = ag/yo (3.1) 

To get a second equation, we can use 1(8.18) with 
(0, + u)=jo(u2) =I%0)X (ue) -J%r,)¥(u), (3.2) 
[ry — w)=fr'(H) = I(r) X(n) —I%(0) ¥(u), (3-3) 


and fo()=fi4)=0, By(r)= ay. Then we get 
[J°(71) —J°(0)][%2 + Yo + $74(%, +y;)] =0.- 
By (2.15) the second bracket is not zero, and hence 


IO) =I%(71) = ag/2V0.- (3-4) 
Alternatively, we can use the fact that 
sYor(s)=| ; J*(t)exp(—st)dt — (n=0,1,...) (3.5) 


is an integral function and therefore regular at s=o. [For.J"(z7) is the N-solution 
of 1(7.2) with w=1 and therefore 


ee at" a,t*"! 
J™(7r)= FA” ——. 
(7) 2A. { as (n—1)! Hues tay}. 





By (2.4) this converges for 0 <r <r,, and |J"(7) | is bounded in (0, 7,) by 
(do |7i"/n! +... + Ja, [M(t—p). 
Since, by the smoothing property of A, J"(r) is continuous for 0<7<7j, it 
follows from a well-known theorem in analysis that s~1j,"(s) is an integral function.] 
By (2.9), (2.12) and (3.2), 
5 oo 1 
Fels) =CofI"O)— IMM 4+ ¥ (= 1)"[Cyed (0) — yl) — 
m=1 . 

pong can contain no negative power of s. Hence J(o)=J%(r,) and thus we 
get (3.4). 


_ When n=1, as in section 7 of Paper I, we get one independent equation, 
viz. 


J*(0)(1 — x9) +I*(73)¥9 —J%(0) x, +’ (71) = 41. 
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By (2.13), (3-4) and (2.14), this can be written 
YolI"(0) + I"(74)] = dag7, + a1. (3.6) 
‘To obtain a second equation we can either use I (8.18) as above, or we can expand 
S~"Jo'(s) =s*X(s)[I%(0)s“* +I M0)] 
— SAY (S*)[I%(ry)s“* +I (74)], (3-7) 
which comes from 1(7.10) with n=1, in ascending powers of s and equate to 
zero the coefficients of the negative powers. On using (3.4) we get 


JMO) —I'74) = a(C, — Dy) = — ay, (3.8) 
where y is given by (2.18). From (3.6) and (3.8),we have 


J"(0) = {ao(371 — Voy) + 41 }/2¥0, } (3.9) 
JI*(71) = {ao(371 + Yoy) + 41}/2Y0- 
Thus, when w = 1, the formulae I(7.21) and 1 (7.22) become 
I(o, +) = aa {out + a, + ao($7,—Yoy)} 
10) {Qype + a, + ao(47, + Voy) }; (3.10) 
X 
I(7,-»)= aw = Agp + a, + 4o(47, + Yoy)} 
09 
= Qype + A, + Ao(37, — Yoy)}- (3-11) 


The number y, can be adi from (2.18) or from the relation 
y= dD, a C, = [ SA] — >A." | , ° (3.12) 
n TE Ty a= 7r=0 
‘To change from the X- and Y-functions used here to the standard functions 
X*(u), Y%() with moments x,%, y,*, such that 


Xy =I, Yo =0, (3-13) 

we have to write 
X(p) = X*(H) — YoulX*(H) + Y%(u)]/(mr° +91°), (3-14) 
Y(u) = Yu) + you [X%(H) + ¥*(H)]/(1* +91") (3.15) 


“in (3.10) and (3.11). After a lengthy reduction the equations 1/8.11) and 1(8.12) 


are obtained, the constants /1(0), J'(7,) being given in terms of the moments 
x,*, y,® by 1(8.19)-(8.21). On comparing the results of Paper I, Section 8, 
with (3.10) and (3.11) it will be seen that it is much simpler to work with the 
X- and Y-functions defined by (2.5) than with the standard functions. Un- 
fortunately it is the latter which have been tabulated. 

4. Section 10 of Paper I.—When w=1, the analysis of Section g holds as 
far as I(9.12), which becomes 


(0) + O7(7,)=1. (4.1) 

By expanding s~'@,*(s) in ascending powers of s and equating the coefficient of 
s-' to zero, we get 

*(0) — (71) = 1 — 2Yp- (4.2) 


Hence 


O%(0)=1—yo=Xy, O(71)=Yo, (4.3) 
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and these are 1(9.14). Thus the equations I(9.15)-(9.18) hold also when w= 1, 
if X(u) and Y() are defined by (2.5). On changing to the standard functions, we 
get the results of Section ro. 

Concluding note.—It has been possible to give only a brief indication of the 
methods by means of which results can be proved rigorously. ‘The treatment 
is based upon N-solutions ; there are few difficulties and the analysis is sur- 
prisingly simple and elegant. The need for complete rigour is shown by my 
own mistake. 

St. Hugh’s College, 

Oxford : 

1957 July 9. 
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‘THE EFFECT OF REFLECTION ON THE DETERMINATION 
OF MASSES OF CLOSE BINARY SYSTEMS 


Alan H. Batten 





(Communicated by Z. Kopal) 


(Received 1957 July 4) 


Summary 


The “ reflection effect ’’ in close binary systems is briefly described, with 
particular reference to the errors it may introduce into the determination of 
masses and mass-ratios of spectroscopic binary systems. 

Recent work on the photometric effects of reflection in systems in which 
both stars are sensibly spherical and revolve in circular orbits is surveyed, and 
its results extended and adapted to the purpose of estimating the correction 
which should be applied to the velocity curves of such systems at any phase. 
To this end, it is necessary to consider the contribution to the observed 
radial velocity made by the rotation of the reflecting star. This is considered 
for a star which rotates about an arbitrarily inclined axis and with an 
arbitrary angular velocity, the only restriction being that the star must rotate 
as a rigid body. 

From these results a phase law for the effect of reflection on the radial 
velocity of a binary component is derived. From this law, the correction to 
be made to the observed masses is derived in terms of the luminosities, radii 
and separation of the two components. In deriving this correction the 
effect of any “‘ spurious eccentricity ’’ introduced by distortion of the velocity 
curve by the phase law is ignored. Finally the magnitude of the correction 
is roughly estimated for the case in which each star rotates with the Keplerian 
angular velocity about an axis perpendicular to the orbital plane. For most 
of the systems for which the theory of this investigation is valid the total 
correction to the total mass of the system would seem to be ~10 per cent of 
the mass of the heavier component. By extrapolation to closer systems, how- 
ever, it would seem that much larger corrections may sometimes be expected. 





1. Introduction.—The so-called “ reflection effect’’ in close binary systems is 
by now well known, one of the earliest comprehensive treatments of it being that 
of Eddington (1) in 1926. The term “‘ reflection”’ is to some extent a misnomer, 
the processes involved being absorption and re-emission in an unknown combina- 
tion with scattering. Nonetheless, “reflection” is a descriptive term, indicating 
that some of the light of each component is intercepted by the other, and re-emitted 
in the direction of the observer. ‘Thus, the hemispheres of the two stars which 
face each other are brighter than the rest of their respective surfaces, and differing 
amounts of the bright hemispheres are revealed with changing phase producing 
photometric effects with which all observers of these systems are familiar. ‘The 
law controlling the “‘ reflection” is not accurately known, but Kopal has shown in 
(2) that Lambert’s Law is at least a reasonable approximation, and on the basis 
of this law he has recently made a thorough investigation (3) of the photometric 
effect of reflection which takes into account the finite radii of both stars and is 
accurate to the extent that both stars can be regarded as spheres; that is to say, 
to terms of the fourth power of their ratio of their radii to their separation. 

Besides the effect of reflection on light curves, however, there is also an effect 
on velocity curves, ‘This arises because the distribution of light over the apparent 
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disk of the star no longer possesses radial symmetry and so the velocity measured 
by the spectrograph is not that of the centre of gravity of the star. The light 
centre is displaced from the geometrical centre of the disk into the bright hemi- 
sphere, and once the radial symmetry of the light distribution is destroyed, the 
axial rotation of the star contributes to the observed radial velocity. ‘The axial 
rotation of the star may increase or diminish the observed velocity according as 
the star rotates in a direct or retrograde fashion (though the latter has never been 
observed), and this correction will depend on the angular velocity and inclination 
of the axis of rotation. It is the purpose of this investigation to discover the 
magnitude of the correction for this effect and so to obtain some idea of the 
error in the estimation of the mass of the system. 

This object has already attracted a number of other investigators. Some time 
ago Kopal (4) examined this effect, considering the illuminating star as a point 
source. In 1954-5, Kitamura (5,6) published two papers in which the finite 
size of the illuminating star was taken into account. However, he derived a 
correction of only one or two per cent to the masses of u Her and yu." Sco, where 
Kopal derived corrections of up to six per cent. Recently, Ovenden (7) has 
suggested that, in some systems, corrections of the same order as the observed mass 
itself should be applied. It seems, therefore, highly desirable that a thorough 
investigation of the range of magnitude of this effect should be made, and, as 
a beginning to such an investigation, we have considered the effect under the 
following restrictions: both stars are to be regarded as spherical, and to be 
revolving in circular orbits; no restriction will be placed on the inclination of the 
axis of rotation of either star, or on their speeds of rotation, except that we are 
precluded by the first condition from considering very fast-rotating stars, and it 
will be assumed that the stars rotate as rigid bodies. ‘The law of “ reflection” 
will be assumed to be Lambert’s law. With these assumptions we shall be able 
to use many of the results obtained by Kopal in (3)—which will be quoted in the 
next section—to obtain a phase law which could, in principle, be used to rectify 
velocity curves before orbital solutions are made, just as light curves are rectified. 
If the spectroscopic binary is also an eclipsing binary, the phase law will not apply 
during eclipses—this situation presents a somewhat different problem, which 
has recently received an elegant treatment at the hands of Hosokawa (8). We 
shall not be able to apply our results to such close systems as Ovenden considered, 
as the assumption that the stars remain spherical would not be valid. We may, 
however, obtain some indication of the order of magnitude of the correction to 
be applied, even for these systems. 

2. Summary of previous results needed in this investigation.—In all that follows, 
unless stated otherwise, the suffix 1 refers to the illuminating star, and 2 to the 
reflecting star, L, m, r stand for the luminosities, masses and radii, respectively, 
of the two stars, R is their separation, 7 denotes the inclination of the orbital plane 
to the celestial sphere, and y is the phase angle measured from superior conjunction 
of the reflecting star. The quantity ¢ will also be useful, where 

cos€=cos¥sin?. 

The basic equation giving the correction Av(e) to be added to the observed 
radial velocity has already been given by Kopal in (4). It is 


| Vn (¢,) ¥ (¢,)do 
Av(e) = — 2(€) (1) 





L,+L,(e) 
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where X(e) is the total visible reflecting surface, and do is the element of this surface. 
The difference in radial velocity between the point (¢, ) and the geometrical 
centre of the disk is given by V2(¢, 7), and .4(¢, 7) is the intensity of the reflected 
light at the point (¢, 7). The coordinate system (¢, 7) will be more clearly 
defined presently. The total amount of light reflected by the reflecting star at 
the modified phase angle « is L;(e). The minus sign adjusts our definition of « 
to the convention that velocities of approach are considered negative. Kopal 
also showed in (4) that Av(e) is independent of the effective wave-length of observa- 
tion, if Z refers to the bolometric luminosities of the two stars. 

It is clear that aknowledge of the quantities V (¢, »), 4(¢, n) and Lp(e) and 


of the nature of the integral operator a ()do will provide us with the solution 


of the problem, and all these, except the Fal may be found in Kopal’s investigation 
(3). His results will be quoted without proof here. 

First we must define the coordinate system to which the coordinates (¢, 7) 
refer. ‘They are spherical polar coordinates defined in the usual way by 


F=Tih, Y= F=THs 
where fy =Cos ¢ sin n, g=sin ¢ sin 7, fg = COS (2) 
and where the origin is the centre of the reflecting star, while the x-axis always 
points to the centre of the illuminating star. ‘The z-axis is perpendicular to 
both the x-axis and the line of sight, and the y-axis is defined so that its positive 
direction points directly away from the observer for «=go’. 

In this coordinate system, Kopal has shown that for any point (¢, ») from 
which the whole of the illuminating star is visible (the “ full-light ” zone) .4(¢, 7) 
is given by os % r } 

Ibn) = pd Palin) + 7B Palen) + Fa Palen) » » «} (3) 


correctly to terms of the order (r,/R)*, where the P,,’s are Legendre polynomials. 
For any point (¢,) at which only a part of the illuminating star is visible 
(the ‘‘ penumbral ” zone) we have 


F(G,9)= =D 4450) + SA HGby2) (4) 
where wu is the coefficient of net darkening, “in 
R 
IN (bs “ai, SC mn y (5) 
R 
” I(bsn)= "et & C2 (HY (6) 
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_ 3 _ 1(%\ , 3(%\ - (2) 
ce= 3-32) +42) 16\7, 


and 


2 
Dp. > _3(%\ | 
Cr=§ 3(2) } (8) 
1/7 
@- +2 
I 
cP=- 


all correct to our order of accuracy. 

The problem of determining the total amount of reflected light thus reduces to 
that of integrating powers of :, over the total amount of reflecting surface visible 
at any time. Again, Kopal has shown in (3) that, except during eclipses, this can 
be done, in the full light zone, by applying the integral operator K (sin ,) to the 
powers of yw, and in the penumbral zone by applying the operator 
K (sin ng) — K (sin ,) where 


cos sing, n—sin-*(sin ™ secd) 
K (sin )u8}=14 [08 (e— 4) cos” 4 [sin +*y dn ap (9) 
e—n/2 sin? (sin sec ¢) 


and 


. Tot? ° To—T 
SIN 7), = 7 , SIN 72 = *R ;, 








n, and 7, being, respectively, the angle made by the inner and outer sets of common 
tangents with the line joining the centres of the two stars. 

Finally, by applying this result to the expressions in equations (3)-(6), and 
duly combining the completed integrals we obtain 


L,Ae)=L, {5 (2) eaSeerise 





3 
+ 3(%) (3 cos*e + 2. cose — 1) 


‘ (3) sin € cos € " (Gr) a 5+(7—5)u, oe ] sine} . (10) 
7 R 57 3-uy, 7 

This completes the survey of previous results necessary for the present 
investigation. 

3. The distribution of radial velocity over the surface of the reflecting star.—We 
wish now to investigate the nature of the function Vp (¢, 7) which is the difference 
at any point (¢, 7) between the observed radial velocity at that point and that of the 
orbital motion of the centre of gravity of the star. As has already been made clear, 
this difference arises from the axial rotation of the star. We shall suppose the re- 
flecting star to have an angular velocity of rotation nw, where w is the orbital 
angular velocity, and mg is any positive number. Since we permit the axis of 
rotation to be arbitrarily inclined to the orbital plane, there is no need to consider 
the possibility of retrograde rotation (m_<o) as retrograde rotation about one pole 
can always be represented as direct about the other. 
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It is impossible to specify directly the inclination of the axis of rotation with 
respect to the coordinate system of the previous section since all three axes change 
their direction with changing «. Consequently, we introduce a new set of axes 
x'y’z’ having the same origin as the previous set, and whose x’-axis coincides 
always with the x-axis, but whose 2’-axis is the normal to the orbital plane through 
the centre of the reflecting star; the y’-axis is now defined uniquely. The angle ¢ 
between the z- and 3’-axes (or between the y- and y’-axes) changes uniformly with 
phase, as the z-axis traces a great circle on the surface of the reflecting star, which is 
inclined at the angle 7 to the intersection of the orbital plane with the star’s surface. 
We define the positive directions of both the z- and 2’-axes so that £ = g0° —i when 
s=g0°; we shal] then obtain, by the standard formulae of spherical trigonometry 








: sec tanytanz 
ami erect == V(1 See 8) 
‘The transformation between any two vectors x and x’ is now given by 
x= Ax’ (12) 
where A I fe) O° 
° cosf —sing 
° sin f cos £ 


In the x’y'z’ coordinate system we shall construct a system of spherical polar 
coordinates (¢’, n’) precisely analogous to the (¢, ) coordinates of section 2 and 
we will now define the position of the pole of rotation by the equations 

¢=9, =H, =o 
so that at any phase ~% the coordinates of the pole in the (¢’, 7’) system are 
[(® —), H) and the components of angular velocity of rotation about the x’, y’ and 
2’ axes are, respectively, 
w, =n, w cos (® — 4) sin H 
w,, = N,w sin(® — 4) sin H (13) 
w,=n,w cos H 
from which, by means of equation (13) it follows that the components about the 
x, y and z axes are, similarly 
w, =n, w cos (®—)sinH 
w, = Nz w sin (® — 4) sin H cos { — nw cos H sin ¢ (14) 
w, =n, w sin(® —)sinH sin{ +n, w cos Hcos ¢. 
If we denote the linear velocities arising from these three components as 
V,, V,, V, the radial velocity V; of each point on the surface is given by 
Vp= —V,cose+V,sine (15) 
where V ,, V, are given by the well known equations 
V ,=2w, — yo. 
V, = XW, — Wy. 
Equations (14) and (15) can then be combined to give us 
V ,(¢, 7) = nergw [cos ¢ sin n{sin ¢ sin (® — ) sin H sin £ + cos H cos {sine} 
+sin ¢ sin n{cos « sin (® — y)sin H sin + cos cos H cos ¢} 
+ cos n{ — cose sin (® — 4) sin H cos {+ cose cos H cos ¢ 


— cos (® — ¥) sin H sin e}]. 
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We can write equation (16), 
V ls 0) = Fp + Foye + Fops (17) 
where 
F, = ngrqwisin € sin(® — 4) sin Hsin f+ cos Hcos {sine} 
F, = ngrw{cos € sin (® — #) sin H sin f + cose cos H cos ¢} 
F = ngrgw{ — cos sin(® — 4) sin H cos £ + cose cos H cos  — cos (® — p) 
sin H sine}. (18) 
All the functions of equation (1) are now known, and the solution of the problem 
now consists only of applying the K-operators to the product V,(¢, ») 4 (¢, »). 
This in turn reduces to applying the operators to terms of the form u7(1 <p <5), 
Hus and p?us(0<p <4) where p is an integer. ‘This is done in the next section. 
First, however, let us define 
1, = cos fs sin? = Cos € 
my = sin yp sini 
My = COSt. 
Then, introducing these, together with equations (13) into equations (18) we obtain 
F,=1(w, mo + w, mo) 
F, =1,(w, mp + w, mo) cote 
F,=12{(w, — w,) my) — w, sine tan} cote. (19) 
4. The K-integrals of p2, p?p, and piys.—The integrals K(sin,){u?} and 
{K(sin n,) — K(sin 7,)}{u?} have been determined for some of the smaller values of 
p, and the results have been used in the summarized data of section 2. Kopal 
showed in (3) that the integrals K(sin){u?} could be represented entirely by 


. four functions of « and 7, (or 2 for K(sin 72){u?}) for o<p <3 and the present 
writer has verified that this remains true forp=4 and 5. ‘The four functions are 


®, = cos“! (sin y, cosec e) 
7 , I. 3 Das “ 
= — —Sin y, cosec e — = sin* y, cosec®« — — sin® y, cosec’«... 
2 6 40 
®, = cos“! ( — sec n, cose) 


I , k 
=m-—€e+ 3 cote sin? y, + g Sin* n, (cot* € + 3) cot € 


I 


+ pam ncate(15 + 10 cot®«+3cotte)... 


®, = cos™! (tan 7, cote) 
I 

6 
F = (sin*® « — sin? y,)"? 


3 


7 , ‘ 
= 5 tan, cote + tan® y, cot® «+ oo ny, cote... 


. or B . ‘ eee R 
mands sin? n, cosec € — 8 sin’ n, cosec* « — = sin® n, cosec’€ ... . (20) 
and the forms of K(sin ,){u?} for o<p <5 are given by : 


K(sin n;){u}} =73 {®, + cos € cos? n, ®; — sin y, F} 


: 2» ‘ ft 
K(sin y4){ui} = ; rscos € (®, — , sin* n,) + ; r5 cos? nF 
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° 
w% 
"9 


K(sin 7){ui}= rib +cos*«)®, + 2 — sin‘ n,) cose D, + sin ny (1 —2sin? 4,)F 


‘ , 2 « —— , ‘ 2; 
K(sinm){ui} = 5 ricose(®,—, sin’ m,) + =r3(1 —sinty,)F— 2 BF 
b 


K(sin 7;){}} = 243 + cos? (sin?« +5))®,+ is — sin® n,) cose D; 


sin 
+ " rr 2°) (4 + sin? n, —8 sin! ee 


. 5, We mie 
K(sin 9,){41} = > cos €(®,— sin’ ,) + a (1 ~ sin’) F 


2r5 : 4r5 
— —(1+sin*® n,) F?-— — F°. I 
31 (1+ ™) 105 (21) 
The form of the integrals K(sin ,){u/{2} may be obtained by substituting in 
equation (9) cos’¢ sin ¢ for cos’¢. ‘The power of sin » in the first integral will 
remain unchanged. It has been found that the same four functions of equations 
(20) may be used to express these integrals as well, and the results for o <p <4 are 
quoted in equations (22). 
, eT ae sine sin — I _ 
K(sin 7,){u{ po} = Be sine D+ manent (sin?y, — 3) ®;— 3 cote sin? y, Fi 


. 


sin € COS € sine 


K(sin ))(Hif2} =72 {seco ®,+ 





(sin* », — 2 sin? », + 1) ®, 


cia 
— ~sin® n, cote F 
4 
ee Se ae - 
K(sin 9;){ujH2} = 5 eet om 91 (3 sin? y, — 5) ®, 
aft. >. 5 
—7r3 (; sin‘ n, cot « — Z sinecose) F 


sine cose, . r2 sine 
K(sin 7;){uj{H2} = 1, (eink +2) + 2 





(2 sin® n, — 3 sin'n, +1) ®, 


2 
a ‘ 
- F (sin? , cote ~ ,8inn sine cose) F 


ok 


r2sine sin? 7, 


K(sin {un} = 5273 resine D+ = 35 (5 sin?», —7)— ~sin* y, coteF 


2rs . ‘ . 2: “ 
+ et + sin? ,) F+ ie (22) 
= 


Finally if in equation (9) we replace sin? +*» by sin” *"» cos » we get the form of 
the integrals K(sin ,){u?43}. It is then immediately seen that the first integration 
yields (sin”+2y)/(p + 2), which taken between the limits of the integration vanishes 
identically. We may, therefore, omit from all further consideration the term 
Fu, in equation (17). 

The integrals K(sin 7) are, of course, obtained simply by substituting sin 7, 
for sin », in equations (20), (21), (22). Again we have to apply the operator 
K(sin n.) — K(sin n,) to two types of integrand which arise from the multiplication 
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of equations (18) and (4). ‘The term Fy, in equation ( 18) produces the following 





integrand ; 
ryLyhy > Cu.» kK “pt (23) 
7R® n=0 é ry : “J 
while the term Fy. produces the integrand 
nbPs & avy R \" , 
7R% ao . ial (24) 


Now Kopal has shown in (3), and it may be readily verified from the data given 
here, that 


: : Yy L + uD R n 
| K(sin n2) — K(sin m1) 4 aR 2 Cr (m1) i 


2 4 UD | » -Yn+l _, r+1 
ao they tne F Uo (at el _ fay” 
a\ R® nay U+1 | ro 


(to the order of accuracy adopted) so that it follows that the integrand (23) will be of 
the fifth order in 7/R, and consequently is negligible within the limits of our 
approximation. Similarly, as the reader may verify, we have 





P ha Fae ee : . 
(K (sin ya) — K(sin 7) }{miva}= —* 5 (sin" 1, — sin" 4H) sine (25) 


‘Taking into account the dependence of the coefficients of (24) on the ratio r,/R it is 
seen that equation (25) is of the order (7/R)* for all n. 

We are now in a position to write down the full expression for the integral in the 
numerator of equation (1). Let us denote its value in the full light zone by J, 
and in the penumbral zone by J, so that 


) ce Vaden) ¥ (4,9) do=1p +1, 


where 


2 3 
1p= Play (7) (1+cos«)*+ = (2) (4[(7 — €) cose + sine] — 3 sin*e) 


7 


I (r\* ‘ . 
a 2) (5 cost « — 12 cos ¢— 40086 +3)} 


( 2 3 
+ Fila\g (3) (1 + cose) sin « — aaa (3) (4[7 —«] sine — 3 sin*€ cose) 


157 
“ sine (1 “amt! 
32 R 





22 
(5 cos*e— 5 cose — 8) — vi sing} (26) 


and 
FL, sine rr? 

} Tom(qu,) Ri 30(1 ~My) + Omuy}. (27) 
The approximation in equation (27) arises from the neglect of one small term in the 
summation and involves an error of ~o-5 per cent. This is, in principle, the 
solution of the problem; there is, however, one further refinement to consider 
before applying the formula. 
5: Secondary reflection.—It has been pointed out that a reflection law which 
includes terms of the form rj73/R* must, to be consistent, include the effect of 
secondary reflection ; that is the effect of light which has been reflected from one 
star back to the other, and thence to the observer. . It is clear that this phenomenon 
effectively increases the luminosity of the illuminating star. Kopal has shown in 
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(3) that the effect of the increase is that L, in his formula for Lp(¢) (equation 10) 
should be replaced by L, + §7[L,/R?. It is obvious that, within our approximation, 
this is equivalent to adding an extra term in rir3/R* to the right hand side of equation 
(10). Since reflected light obeys a different law of limb darkening from that of the 
star’s proper light, one might expect the distribution of intensity over the surface 
of the illuminating star to be changed, with a consequent modification of 4 (¢, 7) 
in the penumbral zone. Such an effect, however, would be small since the res- 
triction we have made that both stars be spherical in fact amounts to an upper limit 
of r/R of ~ort or 0-2 so that reflected light is never more than a few per cent of the 
star’s proper light, and its effect on the distribution of intensity can be neglected. 
Nonetheless, this neglect represents a further danger in extrapolating our results to 
very close binary systems or to systems in which one component is so dark that 
radiation reflected from it becomes comparable in intensity to that of the star itself. 

6. The correction to the masses of binary components. By substituting the results 
of equations (10), (26) and (27) into equation (1), and remembering to write 
L, + 3r°L,/R, for L, we can obtain the explicit form of the correction Av(e) in terms 
of the enalienaiilen ‘on radii of the two components together with the rotational 
velocity of the reflecting component. ‘This formula, if written out explicitly, 
could, in principle, be used to rectify a velocity curve before a solution was made 
from it for the elements of the system. If, further, the notation of equation (19) 
is used then the phase law of the velocity variation is given by 


7 r2 (w,Mo + wo) 
AVe)= = TT) Litalg :) (1 +cose) 





: r2\ 16 31 3 —_— —1)(1 “OS 
+ ism %) oe ee 32\R (5 cose — 1) cos €) 
7m", t Lyrr, 8 
+ 4k Cos € + aL. Re 2 (1 + cose) }. (28) 
‘To obtain the correction Am to the mass of the reflecting component we require 
only to know po () = AK, 


We have at once that 
cose=cos%=0, sine=siny~=1 
and so Mo = Sint. 

The minus sign in front of the right hand side of equation (1) can be ignored 
since we are no longer interested in the direction of the radial velocity, but only in 
the absolute value of the semi-amplitude of the velocity curve. This is always 
increased by the correction AK, unless cosH<o and n,>1. ‘This exception 
corresponds to retrograde rotation, and for cos H = —1, my=1 we have AK,=o0 
However, retrograde rotation is an unknown phenomenon in binary systems, and 
the change of sign in F, would, in any case, ensure that the correction AK is always 
applied in the right sense. ‘I'he general formula, therefore, reduces to 


» {tfte\? 2 (r2\®_ 3 (72\* 1 rr ly 
Fia() “i a(#) 3 (3) + RRL, 


(La/Ly)+ { 2 (%2\*_ 1 ("2 + (4 : ay 
\37\R) ~ 8\R sm a onL,) Rt 


(29) 





AK,= 
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where F, is defined by equation (19) and L/L, by 
L, #"T3 
L,~ AT} 7” 
where 7’, , are the effective temperatures of the two stars. ‘The suffix 2 has been 
used to denote the reflecting star, so that AK for the other star is obtained simply 
by exchanging the suffixes 1 and 2 in equation (29) (except for F', which is adjusted 
by substituting in it , for m2). If rg, ; are expressed in kms and w in radians/sec. 
AK, , are obtained in km/sec in equation (29). ‘The AK’s must be divided by the 
factor sin i to obtain the true velocities of the system. If then Pis the period of the 
system in days, the true value of R in kms is given by 
Rsini = 13,750P(K, + K,+AK,+AK,). (31) 
The correction AR to be applied to any elements already deduced is 
ARsini= 13,750P (AK, + AK.) 
(K,, 2 are here supposed to be the observed semi-amplitudes of the velocity curves, 
and are, therefore, the true values multiplied by sin). 
If Ris then expressed in astronomical units and Pin years and m in solar masses, 
then, 
(m, +m, + Am) sin*i = R°/P? 
= 1-036 x 10-* P(K, + K, + AK, + AK,)* 
if Rand Pare to be expressed in the units of equation (31) but m still in solar masses, 
and 
Am sin’ 1 = 1-036 x 10°’ P[3(K, + K,)*(AK, + AK,) + 3 (K, + K,)(AK, + AK,)*] 
(32) 
since the term (AK, + AK,)? is of the order (r/R)*, provided L,/L, is not too small. 
Also 


m,+Am, K,+AK, 

m,+Am,  K,+AK,’ (33) 

Equations (32) and (33) give the corrections to be made to m, and m, since 

Am= Am, + Am,. 
In one-spectrum systems the corrections to the mass function is 
m} sin*i 

(m, +m)? 
7. The magnitude of the effect.—As has already been stated, the results of the 
previous section are only rigorously applicable to binary systems in which both 
components are sensibly spherical in shape. ‘This implies r/R <0-2 as otherwise 
tidal effects would begin to be significant. In such systems there is no constraint 
on either star to rotate with any particular speed or about any particular axis. 
Furthermore, even if we succeed in detecting rotation of the components of these 
systems we have no means of deriving either the speed or the inclination of the 
axis, but only a combination of the two. ‘The function F, is, therefore, indeter- 
minate. In closer systems the tidal bulge will tend to pull the axis of rotation 
perpendicular to the orbital plane (H =o) and also to make the period of rotation 
coincide with that of revolution (n,=1). In these systems, however, the assump- 
tion of sphericity breaks down. The formulae derived, therefore, are of only 
limited practical use until the effects of reflection from distorted stellar surfaces 
are considered—a task it is hoped to undertake soon. Nonetheless, it may be useful 





= 1036 x 10-7 P(3K3 AK, + 3(AK,)*K,). (34) 
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to employ equation (30) to estimate the order of magnitude of the maximum possible 
effect in a system in which the components rotate with the Keplerian angular 
velocity about an axis perpendicular to the orbital plane (F,;=r,wsini). The 
maximum effect is observed when the reflecting star is completely dark. Then the 
error, Am refers to the illuminating star only, for the illuminating star receives no 
light to reflect. Equation (30) reduces to 


oe #) 2 () 3 (3) 
ls SIN245 — TT —Ii — —_ —f — 
AK, = . (aa ts7\R) ~ 32\R) f 


= (%) - I(r. a 12 5+(m7—5)u, 1 \rir5 
37\R} 8\R 5m* 3, - 1) 


5s m, —_ 
K,= Rw sini. 
m, + ms, 


the determination of masses of close binary systems 531 





while 





In such a system, m, is probably several times larger than m, and, since we are 
interested only in orders of magnitude, we may ignore the factor m,/(m,+my,). If 
we also set u, =0-5, then 





AK, fof lo ro\2 r\3 
_— Be \2°5890 + 05470 7 ~01197(3) -0-0704(%) 


~0°2155 (2) ~0:3272 (Gey } 
J 


m,+Am, (Ky +AK,)! 
ao K3 


where m, is the observed mass of the primary component, whence 
Am, _ 34K, “ (EY + (=) 


~ me K, K, 


while 











which to our order of accuracy gives 
Aah + eet +0000 OY +120 BY +oeusl at 
." +1 7670 5 +2 6818 (2) + 1°7783 (3) +0 8354 (7) 
r2r2 


rr. rs 
— 06465 dl — 17345 7 (35) 


From this we see that for r,/R = 0-1 the maximum correction which may be made 
to the mass of the system is about one fifth of the observed mass of the heavier 
component. For r,/R=~o-3 the formula (35) certainly gives a value of the same 
order as the observed mass itself, as suggested by Ovenden, but it takes no account 
of the distortion which must most certainly exist insuchasystem. If the distortion 
of the reflecting star is small enough for the star to be regarded as an ellipsoid, then, 
although the precise corrections to equation (35) cannot yet be evaluated, they are 
unlikely to affect the order of magnitude of the result. However, many systems in 
which r,/R~o-3 are “semi-detached” systems in which the reflecting star is the 
contact component. As a result, gravity darkening will cause the distribution of 
the star’s proper light to be asymmetric in the opposite sense to the distribution of 
light it reflects, and so equation (35) would break down completely and give a gross 
over-estimation of Am,/m,. 
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By contrast we can also consider the case L, = L, which will give the minimum 
correction for any given value of r/R—this time to be applied to both stars. In this 
case we shall also take m, =m, so that 


K, .=}Rwsinz. 
AK\? AK\3 oo ae 
Again we take u,=0'5; (=) and x) are now negligible quantities to the 


order of accuracy adopted, so that 





Am,, 34412 
mi, Ky, 
and 
Am, : r,\* 
=~ "*s (2) +0°2547 (7) (36) 


For two stars of equal mass and luminosity, the radii may also be taken as equal, 
so that to obtain the total correction to the mass of the system equation (36) should 
be multiplied by 2. 

Substituting in equation (36) r,/R =0-1 we find the correction to the mass to be 
quite negligible, while for r,/R=0-3 it will be 2 per cent of the observed mass of 
one component, or a correction of about 4 per cent of the observed mass of one 
component to the whole system. It should be emphasized that equations (35) 
and (36) are only true if the component stars are rotating with the Keplerian 
angular velocity about an axis perpendicular to the orbital plane. Unless, how- 
ever, the angular velocity of rotation greatly exceeds that of revolution, both 
equations take a reasonable account of the effect of rotation. 

Summarizing, we may say that in certain cases corrections, due to reflection, 
to observed masses of the same order as the masses themselves may occasionally 
be expected, but that such corrections arise in circumstances where theory 
developed in this paper breaks down, so that their precise magnitude may not yet 
be estimated. On the other hand, in systems with r/R~o-1 and in which the two 
components do not differ greatly in luminosity no correction to the masses for the 
reflection effect need be applied. For the general run of systems in which 
0:1 r/R <0-2 and in which the components are reasonably similar in mass and 
luminosity corrections to the total mass of the system of the order of 10 per cent 
of the mass of the heavier component may be expected. This last conclusion 
confirms that reached by Kopal in (4) where the illuminating star was considered 
as a point source. 
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Summary 


Newly determined radial velocities are given for 25 [AU standard velocity 
stars and 343 other stars south of — 26. Many nearby stars and stars of high 
velocity are included in this list. Magnitudes and colours on the Johnson 
B, V system and spectral types on the MK system have been determined for all 
the stars in the main list. Diagrams showing the relation between absolute 
magnitude, spectral type and colour index are given for the stars with known 
trigonometrical parallaxes and from these are deduced photometric parallaxes 
for 142 stars of luminosity class V. Space coordinates and space motions 
have been computed for these stars. 








1. Introduction.—In 1951, as the result of an agreement between the Admiralty 
and the Radcliffe ‘Trustees, one third of the observing time with the 74-inch 
Radcliffe reflector was made available to the Cape Observatory. It was mutually 
agreed that a considerable portion of this time should be devoted to the determina- 
tion of the radial velocities of the later type stars while the Radcliffe observers were 
to concentrate on early type stars and special objects. ‘The first Cape programme, 
begun in 1951, included about 400 stars south of declination — 26° chosen either 
for their large parallaxes, or their large proper motions, or because they lie near one 
of the Kapteyn Selected Areas, Nos. 140-206, and have annual proper motions 
exceeding 0-1. ‘To this initial programme there were later added a number of 
visual double stars suggested by Dr W. H. van den Bos as likely to be worth observ- 
ing, a short list of interferometric double stars selected by Dr W. S. Finsen, 
and a list of 25 stars, which are, or have been, [AU standard velocity stars. ‘These 
latter were included as a result of the request by Dr J. A. Pearce, Chairman of the 
Sub-Commission on Radial Velocity Standards, that the equatorial and southern 
1AU standards be reobserved. As the observations proceeded, certain parts of 
the programme were worked out more rapidly than others and were filled in by 
selecting additional stars from the appropriate Selected Areas. 

The present list gives the results for those stars for which the radial velocity 
observations had been completed before 1956 June 1. Spectral types on the MK 
system have been estimated from the radial velocity plates and accurate magnitudes 
and colours for all the stars have been measured with the photoelectric photo- 
meters at the Cape Observatory. As many of the stars have well-determined 
trigonometrical parallaxes it has been possible to make a new calibration 
of the Hertzsprung—Russell main sequence for the later type stars. 

2. The radial velocity observations.—An account of the Radcliffe Cassegrain 
spectrograph and of its performance has been given by Feast, Thackeray and 
Wesselink (1). The three cameras that have been used for the Cape programme 
are those known as “ a”, “‘b” and “‘c” which give dispersions of 21, 29 and 49 A/mm 
at Hy. Neutral filters having absorptions of 1-3, 2-6 or 5-2 magnitudes have been 
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placed before the slit when observing bright stars in order to obtain an exposure 
time of at least ten minutes so that any errors produced by seeing variations should 
be properly smoothed out. 

The spectra have been measured with standard Hilger long screw machines, 
either with that belonging to the Radcliffe Observatory or with one or other of the 
two at the Cape Observatory. Most of those who have taken part in the observa- 
tions have measured some plates to gain experience but, because of the appreciable 
personal equations that have been detected, only the measures by Dr David S. 
Evans (DSE) and by Mr A. Menzies (AM) have been retained. ‘The selection 
of the lines to be measured and their adopted normal wave-lengths were taken from 
various publications of the Dominion Astrophysical Observatory ; from Con. 
D. A. O. Nos. 4 and 10 for the ‘‘a”’ and “‘b”’ camera plates, from Con. D. A. O. 
Nos. 11 and 12 for the ‘‘c”’ camera plates. Later, when Publ. D. A. O., 9, No. 3 
became available, the selection of lines given in it was used for reducing the spectra 
of F4-MB8 stars observed with the ‘‘b”’ camera. ‘The Victoria standards have 
proved themselves reliable and self-consistent in every way. 

All the data available at 1955 June 1 were used to derive the systematic 
differences between the two measurers and the corrections needed to reduce their 
measures’to the system of the Mt Wilson General Catalogue of Radial Velocities. 
At the same time eleven conveniently placed “‘reference’’ stars were chosen. 
Seven of these are LAU standards and are indicated by an “‘ R”’ number in Table I, 
while the remaining four are given at the end of ‘Table Il. ‘These stars are being 
observed frequently, usually one such star per full night of observation, and 
each of their spectra is being measured by both of the measurers in order to keep 
a continuous control on the necessary systematic corrections. Following a 
suggestion by Feast, Thackeray and Wesselink, the radial velocities of the brighter 
minor planets are being observed in an attempt to evaluate the correction required 
by the observed system of radial velocities to reduce it to an absolute one. 

The systematic corrections in km/s that have actually been applied to the 
measures are as follows :— 


7 “ ” re 
DSE +1°7 ek +0°7 +05 —-I'2 -o9g 
AM +1°7 +02 +07 —25 -—O5 


The first set of figures refers to plates taken before 1955 June 1 and the second 
to plates taken during the year between 1955 June 1 and 1956 June1. ‘The “a” 
dispersion was not used during this latter period. These corrections correspond 
to a linear displacement on the plate of a micron or less. 

The velocities of the 25 IAU standard stars given in Table | were determined 
in 1954 and 1955 from a minimum of six plates taken with either the “‘a”’ or the 
‘‘b’’ dispersion. Each plate was measured by both measurers and these measures 
were combined after the appropriate systematic corrections had been applied. 
Probable errors of the finally adopted velocities were derived both from the 
probable errors of the individual plates and from the interagreement of the 
separate measures. In every case the probable error of the adopted velocity 
derived by the first method is t+ o-2km/s. The probable error estimated from the 
interagreement of the individual plates is + 0-4km/s, or less, for 22 of the stars, 
the actual average being +0:25km/s. For HD 182572 it is + o-8km/s and for 
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HD 119971 and HD 223311 +05km/s, so that the velocities of these stars are 
probably variable. If, omitting these stars, the double star HD 51250 and the 
star HD 173009 for which there is an unduly large and unexplained difference, we 
compare the new velocities with those of a quality given in the Mt Wilson Catalogue 
we find that for 17 stars the average difference between the two sets of velocities 
is ookm/s. This is the result to be expected since the systematic corrections 
to the Cape measures were derived principally from a consideration of these stars. 
A more comprehensive comparison with the Mt Wilson Catalogue is given later. 

The measured radial velocities of the programme stars and their probable 
errors are given in the eighth column on Table I]. Each depends on a minimum 
of four plates and usually rather more. ‘The probable errors are the mean of 
those deduced from the interagreement of the plates and those deduced from the 
formal probable errors of the individual plates. ‘The cases in which the first 
estimate exceeds the second sufficiently to justify suspicion of variability in the 
velocity are indicated in the Notes. The last column of Table II gives the actual 
number of plates measured for each star. 

A fairly large number of stars in Table I] appears in the Mt Wilson Catalogue. 
Some were included in the first selection of stars to be observed in order to give a 
strong link between the new velocities and previous work. Using all the stars in 
Tables I and II for which a radial velocity of a or 6 quality is given in the 
Mt Wilson Catalogue, and arranging them in the usual groups according to 
spectral types, we get as the systematic differences between the observed velocities 
and the values given in the Catalogue the following : 


Sp. > : 

Gest B A F G K M 
a +1:2(1) -—1:7(3)  00(23) +0:1(31) +0°7 (8) 
b +28(1) +05(6) -—1:2(6) —1:6(9) —1-4(16) —1-0(4) 


a+b = +2:B(1) +06(7) —1-4(9) —0-4(32) —0-4 (47) +04 (12) 


the differences being in km/s and in the sense Mt Wilson Catalogue-Cape. 

In whatever way the data are handled, they produce a feeling of some con- 
fidence that any differences that may exist between the observed system of these 
new velocities and the system of the Mt Wilson Catalogue do not exceed the 
probable errors of the material under discussion. 

3. The photometric observations.—All the stars being observed for radial 
velocity have been included in a working list of about 4000 bright southern stars 
for which the magnitudes and colours are being measured with a photoelectric 
photometer attached to the 13-inch Cape Astrographic Refractor. The normal 
routine of measurement at the telescope has been blue sky, blue star, yellow star, 
yellow sky, yellow star, blue star, blue sky, the individual exposures being of 
approximately 30 seconds duration. The net deflexions (star-sky) are read off 
from the Brown Recording Potentiometer chart directly in magnitudes to the 
nearest o"-o1 by means of a perspex scale specially engraved for the purpose. The 
first, last and every third intermediate star observed is a standard. ‘The working 
list is so arranged that all stars between — 4° and — 64° declination are observed 
at 30° zenith distance, while the relatively few stars in the programme outside 
these limits have been observed near culmination and compared with E region 
standard stars at equal altitudes. Where necessary, differential extinction 
corrections have been applied, the assumed zenithal extinctions being 0™-35 for 
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the blue and o”-20 for the yellow, which correspond to average photometric 
conditions at the Cape. In actual practice, extinction corrections have to be 
applied to relatively few observations and very rarely exceed o"-o1. Each star 
receives a minimum of four observations on separate nights and the probable 
errors of the resulting magnitudes and colours as deduced from the internal 
agreement are +0"'-005 and +0"-004 respectively. 

The stars that are being used as standards for the “‘ 30°”’ portion of the photo- 
metric programme are those observed in the Cape “ Bright Star Programme”’, 
the results of which are given in Cape Mimeogram No. 1 (1953). ‘The Fabry 
magnitudes (BPg) and the photoelectric colours (Cpe) given there were approxi- 
mately transformed into the natural colour system (b,c) of the Astrographic 
observations by means of the relations 


b= BPg—0-18 Cpe 
b-y=c= I°Il Cpe — 9°35 
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after the following small zero point corrections, which resulted from a special 
series of observations, had been applied : 


+o"-o1 to the magnitudes of all stars between HR 1298 and HR 2451 
and between HR 4441 and HR 5068 

to the colours of all stars between HR 7063 and HR gog8 

—o™-o1 to the magnitudes of all stars between HR 3165 and HR 3940 
and between HR 7063 and HR 8204 

to the colours of all stars between HR 3165 and HR 4216 


An examination of the results of the first three years’ observations for this 
programme showed that the magnitudes and colours that had been thus derived 
for the standard stars were essentially reliable but that further small corrections, 
which could be partially explained as due to a “‘bandwidth”’ effect, were 
necessary. ‘These were as follows : 


c Ab c , Ac 
Mi m m th 
—o'ol —o°o! 
—0°50 —0'50 
0:00 0-00 
— 0°25 +050 
+00! +O°0I 
+0°65 
0-00 
+ 0°90 
—o-’ol 


The relations used for transforming the observed Astrographic visual magnitudes 
and colours into the Johnson V, B— V system are 
Y ay~Oy 
B-—V=0-92¢+0°52 
These were derived by a direct comparison of the magnitudes and colours of 


30 stars as observed with the Astrographic with those given by Johnson in Annales 
d’ Astrophysique (2) which, for several of the stars in question, differ appreciably 


37 
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from those given in the original paper of Johnson and Morgan (3). The average 
differences between the Johnson and Cape measures are 0-017 in magnitude and 
o™-o12 in colour, which are rather larger than might be expected from the quality 
of the data involved. 

The magnitudes and colours given in the fifth and sixth columns of Table I] 
are, with very few exceptions, derived from recent observations with the Astro- 
graphic Refractor. For a few stars brighter than 5"-o and south of — 64°, the 
magnitudes and colours have been deduced from the earlier observations given in 
Cape Mimeogram No. 1 according to the precepts given above. For the bright 
stars north of —64°, which have been used as standards for the normal “30°” 
programme, the magnitudes and colours given are those deduced from the actual 
Astrographic observations. 

The diaphragm of the Astrographic photometer has a diameter of 1’-5 so that 
it has not been possible to measure individual magnitudes and colours for a number 
of stars that are too close together to be effectively separated. Stars for which 
the magnitudes and colours fail for any reason to come up to the usual standard of 
accuracy are indicated in Table II by : and are referred to in the Notes. 

4. The spectral types.—The spectral types and luminosity classes given in the 
seventh column of Table II have been estimated from the radial velocity plates 
by Dr David S. Evans. They are intended to be on the MK system. The 
types were first estimated using the standards and criteria of the MKK atlas (4) 
and then the modifications set out by Johnson and Morgan (3) were applied. 
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Fic. 1.—Relation between colour and spectral type for Class V stars. Continuous line : Cape 
results; dotted line: Morgan and Fohnson. 


The general relation between spectral type and colour index is shown in Table III 
for those spectral classes for which there are at least three stars in Table II. The 
actual number of stars used for forming the mean colour index is given in the third 
column and the total range in colour index for these stars in the fourth. The 
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next three columns show similar data for Johnson and Morgan’s work taken from 
Johnson’s paper (2). The correspondence between the two sets of data is far 
from perfect, but can probably be considered as satisfactory if the limitations of 
the available data and the variations of colour index within the same nominal type 
are taken into account. 
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Fic. 2.—Hertzsprung-Russell diagram for stars of known parallax. The lower portion shows the 
breakdown of the data for Class V stars in terms of the probable errors of the absolute magnitudes. 


Fig. 1 shows the relation between colour index and spectral type for Class V 
stars, the only class for which the present data is at all extensive. The dotted line 
shows the intrinsic colours for class V stars as given by Johnson and Morgan in 
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Table 14 of their paper (3). The continuous line refers to the present Cape results 
which, for the later type stars, are based on many more stars. 

5. Photometric parallaxes for stars of Class V.—As a result of their initial 
selection, many of the stars in Table II have well determined trigonometrical 
parallaxes. Fig. 2 shows a plot of absolute magnitude against spectral type for 
the stars for which the probable error of the absolute magnitude is less than: 1"-5 
and which are not known to be multiple. The dashed and continuous line 
indicates the’main sequence defined by the stars of Class V. The lower half of 
the figure shows a breakdown of the data for these Class V stars according to the 
weight of the absolute magnitudes. Table IV gives the mean of the weighted 
values of the absolute magnitude for each spectral type for which three or more 
stars are available, together with comparable data from Keenan and Morgan (5), 
from Gliese (6) and from Allen(7). The general agreement between the various 
determinations of absolute magnitude is very satisfactory. 
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Fic. 3.—The full line represents the zero-age main sequence of Sandage (Ap. F., 125, 437, 1957); 


the dotted line the mean relation between My and B-V used in the derivation of photometric 
parallaxes. 


Fig. 3 shows a plot of the absolute magnitude against the B — V colours for the 
same Class V stars as were plotted in Fig. 2. The continuous line represents the 
main sequence for stars of zero age according to Sandage (8). A careful examina- 
tion of the data indicates that, with very few exceptions, the dispersion in the 
observed absolute magnitudes for any given colour can be completely attributed 
to errors in the observed trigonometrical parallaxes. If, then, we assume this 
dispersion to be purely the result of those errors, we can use the observed mean 
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curve for deriving the absolute magnitude for any Class V star for which the 
B-—V colour is known. 

The result of applying this process to the non-double Class V stars of Table II 
of types F3 and later is shown in Table V. The mean curve adopted is shown in 
Fig. 3 by a dotted line and is given implicitly in the second and third columns of 
Table V. The fourth column of Table V gives the “ photometric” parallax 
derived from the difference between the absolute magnitude given in the third 
column and the observed magnitude. The observed trigonometrical parallaxes 
given in the fifth column are taken from the Yale General Catalogue (g) or from 
Cape lists (10) which were not included in that catalogue. 

A comparison of the two sets of parallaxes reveals only five possible stars 
(R7, 94, 189, 208 and 254) for which there is any reason to suspect a difference 
that cannot be immediately attributed to the trigonometrical parallax observations, 
either because these have been observed at one observatory only (usually the Cape) 
or because the values determined at different observatories differ widely and 
bracket the value derived from the colour of the star. Stars R7, 94 and 254 
appear to be o"-g, o"-6 and o"-8 brighter than would be judged from their colour, 
while 189 and 208 are o"-7 and o"'-4 fainter. Star 189, which is No. 3669 in the 
Yale Catalogue, has been noted as a subluminous star by Eggen (11) while star 
208 is part of the system of 36 Oph from which it is separated by 3’. 

The space coordinates in parsecs and the space motions in km/s corresponding 
to the photometric parallaxes are given in the last six columns of Table V. They 
are relative to the Sun and are referred to axes defined by 
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l= 148°, b=0° 
l= 58°, b=0° 
b=go° 


‘The annual proper motions used for computing the space motions are given in the 
fifth and sixth columns of Table V. In general they are the weighted means of the 
values given in the Albany General Catalogue (12) and the values given in more 
recent catalogues or those derived in the course of the trigonometrical parallax 
measures. 

Stars 2, 28, 32, 68, 81, 200, 227, 278 and 323, which do not appear in Gliese’s 
recent catalogue (13), may be closer than 20 parsecs and may be worth the 
attention of trigonometrical parallax observers. 
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No. 5, 1957 Fundamental data for southern stars 543 
Taste II 
(1950) 
No. HD R.A. S. Dec V B-| Spec. Vel. (km/s) Pi. 
: ws 

1 105 | 00 03°3 | 42 02 | 7°52 -o'58 | Go V 3z;2+1°0 7 
2 834 10'r | 27 08 | 7°98 -o'95 | Ko V 6 l 4 
3 1187 13°6 | 31 43 | 5°68 -1°36 | Ks III 26°2 oO'4 6 
4 2151 ag°1 | 77 ga | aot +o61 | G2 IV 23°3 oa 6 
5 2885 29°3 | 63 14 A2 V 8:2 o's 5 
6 3443 34°38 | 25 03 | 5°58 0-70 | G5 V 16-o+0°8 4 
F 4 3737 37°4 4! 56 8-48 0°54 Fs \ i a o's - 
8 4229 39°6 | 85 58 | 6°81 1°28 | Ks III 12°O 0O°5 _ 
9 4309 42°3 | 74 32 | 7°58 | +047 | F8 V 25 I 4 
10 4293 42°6 | 42 57 | 5°94 0°28 7 Vv I 2 5 
11 5133 50°6 | 30 38 | 7°17 0°94 | K3 V I1r-oto'9g 6 
12 6269 | o1 009 | 29 48 | 6°29 0:94 | G5 IV 56:5 06 4 
13 01'3 | 60 22 Bo ee Fs V 12‘2 1°0 5 
14 0334 o1°3 | 60 22 —_ aia Fs V 13°70 1°O 5 
15 6623 03°5 | 72 00 |. 7°36 rio | K3 III 6-6 O'5 6 
16 6996 o7°3 | §© 52 1 7° 0°44 | Fs IV 12°6+0°6 7 
17 8224 18°5 | 57 36] 7:01 0°52 7 V 20°I O'S 4 
18 8391 20°2 | 43 §2 | 7°01 0:30 | Fo IV 4°2 1°4 5 
19 8638 22°4 | 28 05 | 8:30 0-69 | G3 V oy a Ft 4 
20 9379 28:9 | 59 55 | 7°96 0°57 | Go V - 2 FO 5 
21 9468 29°7 | 59 51 | 7°97 +o'45 | F5 V 4°rt 1°! 6 
22 9770 32°7 | 30 10 | 7°% -o'92 | K3 V 94°. 2° 7 
23 10360 ee aa a 4 Ko \ 20°1 O'4 5 
24 | 10361 37°9 | 5® 27 | 5°07 oe Ko V 22°0 O'4 4 
25 | 11695 51°7 | 46 33 | 4°4 ~158 | M4 III 570 04 4 
26 12438 59°0 | 30 14 | 5°35 o'88 | Gs Ill 26°4+03 5 
27 13445 | 02 08°4 | 51 04 | 6°12 0°81 | Ko \V 55°3 06 5 
28 13435 08:5 | 28 27 | 7°06 1°02 | K2 V 39°5 08 5 
29 14412 16°7 | 26 11 | 6°34 o'72 | Gs V + 79 06 5 
30 15339 24°9 | 46 13 | 7°14 111 | Ko IIl 16 0°5 4 
31 15590 27°2 | 42 18 | 7°97 066 | Gs IV 36°8+0°6 7 
32 16591 36°3 2 07 | 7°25 0:94 | Ko pd 48°5 o'8 8 
33 | 417155 418 | 46 39 | 9°04 | +105 | Ks V 347 10 | 4 
34 17865 489 | 44.17 | 8:15 055 | F8 V 2 ce ee oe 6 
35 17926 49°8 | 31 of | 6°39 +047 | F8 IV-\ 8-1 04 4 
36 18819 58°6 | 27 50 | 7°60 0°58 | Go = 1'6+0'8 5 
37 | 18907 59°5 | 28 17 | 5°88 | +0°79 | G5 IV 39°9 OF 4 
38 19743 | 03 06°6 | 61 54 | 7°06 -0'93 | G5 IV 6-5 o6 4 
39 20052 o0g'2 | 62 33 | 8°32 0°57 | Go V ; 26:3 1°3 5 
40 20010 o9'9 | 29 11 | 3°85 ost | F8 IV — 19°4 0%3 4 
41 20766 16°7 | 62 46 | 5°53 +063 | G2 V 12°3+0°5 5 
42 20807 17°1 | 62 42 | 5°26 058 | Gr V I1°3 0°4 4 
43 20855 17°7 | 58 10 | 7°43 tr-o5 | Ko IV — 95 O'S 5 
44 21749 26:1 | 63 40 | 8-08: | +1°16:| Ks Vv 57°3 0°9 4 
45 21899 28's | 41 32 | 6-11 | +0°47 F8 V 16°2 0°4 4 































































544 David S. Evans, A. Menzies and R. H. Stoy Vol 117 
TaBLe II (cont.) 
(1950) 
No. | HD R.A. S.Dec| | B-| Spec. Vel. (km/s) | PI. 
h m , 

46 | 22359 | 03 31°7 | 60 45 | 7°55 | +0-48 | F8 IV-V o-7t1°2 5 
47 22262 31°9 | 31 15 | 619 0-47 | F5 V 32°6 05 5 
48 | 23484 42°3 | 38 27 | 6-98 -o'89 | Ko V 31°9 «1°4 4 
49 | 23817 43°6 | 64 58 | 3°82 +116 | Ko IV 48:6 06 4 
50 24450 47°4 | 74 40 | 8:00 +037 | Fo V 54 08 5 
51 24512 48°0 | 74 24 | 3°31 + 1°61 | Mo Ill 14°7+ 04 4 
52 24331 48°9 | 42 43 | 8-59 >o-o1 | K2 V 22°4 O'S 4 
53 | 25587 | 04 00°7 | 27 37 | 7°39 | + 0°52 | F8 V 78 o°5 5 
54 | 25945 03°6 | 27 47 | 5°57 | +0°32 | Fo V 66-4 of | 5 
55 27274 14°7 | §3 26 | 7:62 +114 | Ks V 23°2 08 5 
56 | 28454 25°6 | 47 03 | 6-10 | +045 | F8 V 1764+ 1-1 4 
57 | 29875 39°O | 41 58 | 4°45 | + 0°33 | F2 V 17 04 | 6 
58 30003 39°5 | 59 02 | 6°53 +068 | Gs V IO'l O°5 4 
59 | 30361 43°3 | 47 30 | 8-33 o59 | Gi V 59°11 o8 | 4 
60 31081 45'8 | 76 24 | 7°73 +o-48 | F8 IV-\ o9 O04 5 
61 30694 46'1 | 44 21 | 8-12 +o-32 | F2 IV g:2+1°%3 4 
62 | 31093 49°7 | 34 59 | 5°85 | +0710 | Ar Vn 23. «2 4 
63 31975 53°9 | 72 29 | 6°26 +o'50 | F8 V 24°9 O°4 5 
64 | 33563 | 05 03-9 | 76 42 | 7°52 | +047 | Fs V 188 08 4 
65 33262 04°7 | 57 32 | 4°69 +o-51 | F8 V 2°0 O° 4 
66 | 34554 15‘1 | 31 20 | 7°46 0-45 | Fo V re her fe oe ak 4 
67 36137 25°9 | 46 08 | 7°70 o'40 | F3 V o8 1°3 4 
68 36435 27°1 | 60 27 | 6-95 +o-77 | Gs V 12°9 04 5 
69 | 36553 28°8 | 47 07 | 5:46 + 0°62 | G3 IV 16°5 0-4 5 
70 | 37655 36°5 | 43 00 | 7°43 o'59 | Go V 2s 3s 7 
71 37706 36°8 | 46 08 | 7-32 | +0°77 | Gs V 11°4+0°5 4 
72 | 38940 45°5 | 45 40 | 7-40 | +0-49 | Fo IV-V 35°55 O8 | 5 
73 | 39280 47°9 | 44 42 | 7°74 o'94 | G8 IV 29°5  0°7 4 
74 | 39655 59°3 | 44 02 | 855 | +035 | F2 V 293 08 | 5 
75 39962 52°3 | 42 14 | 7°96 | +039 | F2 V 8-2 08 4 
76 | 40292 53°7 | 52 39 | 5:28 | + 0-31 | Fo V 28-4+0°7 4 
77 | 41172 | 06 oo'1 | 27 25 | F-41 +044 | Fs IV-V 54 08 5 
78 | 41534 02°5 | 32 10 | 5°64 | —o-20 | B2 V +1051 1°8 5 
79 | 43834 117 | 74 44 | 5°07 | +071 | Gs V 348 o7 | 4 
80 | 48676 4°°9 | 42 31 | 7°94 | + 0°59 | Go IV-V 19°5 oO'9 sf 
50223 45°5 | 46 34 | 513 | +045 sy ¥ 20-9+ 0°8 4 

82 | 50503 49°7 | 47 16 | 7°32 | +1°16 | K2 III 31°3. 08 4 
83 50806 51°6 | 28 28 | 6-02 | +0-72 Gs IV 70°9 06 4 
84 | 51608 53°9 | 55 12 | 8:17 | +077 | G7 V 40°I O'7 5 
85 | 52395 57°8 | 29 38 | 7°79 | +0°53 | Go V 18:0 0-9 4 
52698 59° 25 53 | 6-72 | +o-91 | Ko V 12°6+0°6 4 

87 53349 | 07 00°3 | 58 52 | 6:02 | +0-29 | Fo V I0‘O 0-7 4 
88 | 56737 13°7 | 60 59 | 7°16 | +.0-41 | Fz V 8-7 06 | 4 
89 | 57095 16°1 | 46 54 | 6-69 | +0-99 | K2 V 59°90 08 3 
90 | 58895 23°3 | 58 24 | 659 | +0-70 | Gs IV + 197 o5 | 4 























































































No. §, 1957 


Tas_e II (cont.) 


Fundamental data for southern stars 





























| 

(1950) 
No. | HD R.A. 5S. Dec V B-V Spec. 

| h m ’ 
gI @7 27°79 | 43 12 Gs V 
92 62644 41°4 | 45 03 | 5°06 0°77 | Gs IV 
93 62850 41°7 | §9 11 | 7°18 +063 | G2 IV 
94 | 63077 43°7 | 34 04 | 5°36 | +0757 | Go V 
95 | 63685 45°7 | 61 19 | 7°42 074 | G5 V 
96 63744 46°9 | 46 57 | 4°71 1:07 | Ko III 
97 | 64379 50°4 | 34 35 | 5°03 044 | Fs V 
98 | 68978 | 08 11°6 | 31 35 | 6°68 o°61 | Gs IV-\ 
99 71243 19°9 | 76 46 | 4:06 0°38 | Fo IV 
100 | 74543 39°O | 73 54 | 6°82 104 | Ko IV 
101 74558 411 | 46 38 | 6-94 0:25 | A7 III 
102 74576 41°4 | 38 42 | 6°55 0-93 | Ki V 
103 74842 43°0 | 42 27 | 7:20 +074 | G5 V 
104 | 76483 53°4 | 27 29 | 4°87 | +013 | A3 V 
105 77140 57°2 | 47 02 | §°17 0:25 | Fo III 
106 81830 | og 24:2 | 61 44] 5°76 +o:16 | A2 Vn 
107 81783 24°3 | 47 33.1 7°71 1:28 | K3 IIl 
108 82241 27°6 | 44 19 | 6:96 +o-48 | F8 V 
109 82434 28:7 | 40 15 | 3°57 +o-35 | F2 IV 
110 84117 40°O | 23 41 | 4°95 +o:52 | Go V 
111 84850 44°4 | 58 34 | 6:22 0°44 | Fo IV-\V 
112 85512 49°1 | 43 15 | 7°65 118 | Ks V 
113 86006 52°6 | 45 30 | 8-16 +o71 | Gs IV 
114 87638 | 10 03°4 | 33 09 | 7°00 +o-30 | F3 IV 
115 87783 04'2 | 47 08 | 5:07 o'88 | Ko IV 
116 88201 07°3 | 32 36 | 7°46 0:54 | Go IV 
117 88746 11°0 | 47 14 | 8-13 +o-78 | G8 V 
118 88742 | 11-2 | 32 47 | 6°38 +o58 | Gr V 

| 119 go518 24°0 22 6°13 +213. | Ki III 
| 120 91881 33°7 | 26 25 | 6-28 0-45 | Fé V 

121 | 91981 34°2 | 47 35 | 7°31 0°56 | Go IV 
122 92449 37°3 [ 55 21 | 4°58 o'79 | G2 Il 
123 93497 44°6 | 49 09 | 2°69 -o-go | G5 III 
124 | 94444 51°2 | 44 09 | 8-12 049 | F8 IV-V 
125 | 95456 58:3 | 31 34 | 6°07 | +050 | Go V 
126 96557 | 11 04°8 | 32 19 | 6°58 0°34 | F2 IV 
127 | 96700 05°5 | 29 54 | 6°54 058 | G2 V 
128 | 97840 12°6 | 33 03 | 7°00 0°34 | Fs IV 
129 | 98718 18-7 | 54 13 | 3°89 | —o'17 | Bs Vn 
130 | 98818 19°4 | 60 57 | 7°35 1:13 | Ko IV 
131 99279 aa*¢ | Gr a2 | 7°22 127 “tie 
132 | 100623 92°t | 33 34.1 $°97 +080 | Ko V 
133 | 100733 32°8 | 47 06 | 5°75: | +1°66:| M3 III 
134 | 100773 32°9 | 60 37 | 6°57 | +0°35 | F2 IV 
135 | 101349 37°0 | 48 12 | 8-99 | +0°84 | Ko V 


Vel. (km/s) 
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546 David S. Evans, A. Menzies and R. H. Stoy Vol. 117 
TaABLe II (cont.) 
(1950) 
No. HD R.A. S. Dec V B-| Spec. Vel. (km/s) Pi. 
h m 8. 

136 | 101408 | 11 37°4 | 45 32 | 7°25 +104 | G8 IV 28'1+0°7 5 
137 | 101581 38°6 | 44 08 | 7°77 tro7 | Ks V 13°3 08 4 
138 | 102350 44°1 | 60 54] 4°10 ogo | G3 III 3°0 04 4 
139 | 102365 44°1 | 40 14] 4°91 0°64 | Gs V 14°7 0°6 4 
140 | 103746 54°2 | 46 48 | 6-27 +039 | F3 IV-V 76 16 5 
141 | 103932 55°5 | 27 25 | 6-97 116 | Ks V 40° 7+09 | 4 
142 | 103975 55°7 | 47 42 | 6°77 0-50 | Go V 21, 06 4 
143 | 109200 | 12 30°7 | 68 29 | 7-12 0°82 | Ko V o°3 «O90 5 
144 | 110287 38°6 | 45 52 | 5°84 1°52 | K3 Il 6-9 oO'4 6 
145 | 110304 38°38 | 48 41 | 2:17 -oror | Ao III I'g o79 4 
146 | 110477 40°2 | 60 52 | 7°77 0-42 | Fo IV Igt+i'l 6 
147 | 110829 2°7 | 60 42 | 4°71 ~1'06 | Kr III 52 06 5 
148 | 112164 2°2 | 43 §3 | 5°89 0°62 | G2 IV 32°0 O'5 4 
149 | 115577 | 13 15°6 | 28 04 | 6°81 o-98 | G8 IV I§2'2 0O°5 4 
150 | 116243 20°6 | 64 16 | 4°53 0°84 | Gs III-IV I2°9 04 5 
151 | 117440 28°1 | 39 09 | 3:90 1:20 | G8 III 1°3+0°3 5 
152 | 117558 28°8 | 27 51 | 6°47 toro | Ao III-V 4 2 5 
153 | 119070 39°0 | 48 13 |] g18 +o'77 | Gs V 570 oO” 4 
154 | 119629 42°4 | 48 33 | 6°77 to's: | F8 V — 21°7 £0 4 
155 | 120592 48°4 | 48 03 | 7-4 +-O°7 Gs V 6-9 I1°0 4 
156 | 120593 48°5 | 48 03 | 7°5 +or5 | F6 4°5+14 | 4 
157 | 121416 53°2 | 46 21 | 5:82 +115 | Ko IV — 10 o6 7 
158 | 121852 56'1 | 45 14 | 7°38 +o50 ¥ 38 1° 5 
159 | 123651 | 14 069 | 46 o2 | 8-21 +o°53 | Gr V 220 I'I 6 
160 | 123797 07°8 | 48 33 | 6°63 oor | Gs IV - 63:2 oO4 4 
161 | 125072 15°5 | 59 08 | 6°66 1°04 | K3 V -1t1'1 4 
162 | 125932 20'2 | 27 32 | 4°80 1°34 | Ks III 20°9 O'S 6 
163 | 125968 20°5 | 27 36 | 7°75 0°65 | Gs IV-V 25°9 O'7 5 
164 | 128266 34°I | 45 55 | 5°36 1°00 | Ko III - 159 O73 6 
165 | 128582 35°9 | 46 22 | 6°08 0-48 | F8 IV-V 12°99 O'5 4 
166 | 128931 37°5 | 29 05 | 7°82 +o'44 | F3.V 27°2+0°9 4 
167 | 128898 38°4 | 64 46 | 3:18 +026 | Fo III 65 o4 5 
168 | 129178 39°0 | 29 09 | 8-13 +o°47 | F6 V a6 5% 4 
169 | 129422 41°3 | 62 40 | 5-36 +o:29 | AZ Vn 5 3 6 
170 | 129747 42°5 | 45 40 | 8-49 +066 | G2 V §3°6 1°2 5 
131342 51°7 | 59 55 | 5:22 116 | Kr III 13°8+04 6 
172 | 132301 56°5 | 43 37 | 6°57 +046 | Fs V 21°99 08 4 
173 | 134060 | 15 06-7 | 61 14 | 6-28 +o61 | G3 IV 38°30 8 4 
174 | 134331 ©7°4 | 43 32] 7:o1 +o-61 | G5 V — 38 o7 5 
175 | 135379 13°6 | 58 37 | 4:07 | +o-10 | Az V g1 08 5 
136352 18-4 | 48 08 | 5-67 | +063 | G2 V — 7o:9+0°8 4 
177 | 138549 30°8 | 30 51 | 7°96 | +070 | Gs V 14°4 0°5 4 
178 | 139127 34°7 | 42 24 | 4°36 | +1-44 | Mo III — 64 o%4 6 
179 | 139664 37°38 | 44 30 | 4°64 | +0°39 Fs IV-V — 48 1°0 6 
180 | 139961 39°4 | 44 47 | 8-86 | +o0-10 | Ao V +145 3 5 

























































No. 5, 1957 Fundamental data for southern stars 






























TaBLe II (cont.) 




















(1950) 
No. HD R.A. 5S. Dec V B-V Spec. Vel. (km/s) Pi. 
h m > 
181 | 140690 | 15 43°3 | 43 05 | 8-10 +064 | G5 IV + 6o°5+1°5 4 
182 | 140901 44°2 | 37 46 | 6:03 +o-7o | G6 V — §9 o8 5 
183 | 141891 50°7 | 63 17 | 2°86 +o:28 | F2 IV + 26 o7 6 
184 | 142254 51°8 | 42 28 | 6-69 +o0-38 | Fo V — 1770 O99 4 
185 | 142529 53°5 | 48 o1 | 6:30 9°36 | F2 V + 183 oF 5 
186 | 143474 59°5 | 57 38 | 4°64 | +023 | As V - 8 +2 4 
187 | 144628 | 16 05°7 | 56 19 | 7-12 | +084 | K3 V + 37°9 10 6 
188 | 145158 07°8 | 45 12 | 6°65 +o-47 | F8 V — 5°70 o'9 4 
189 | 145417 09'8 | 57 25 | 7°54 to:80 | Ko V + 1°39 «I'2 5 
190 | 145544 10°8 | 63 34 | 3°85 ti-12 | G2 II 5°70 04 5 
gt | 146624 15°2 | 28 29 | 4°79 0:03 | A2 V 13;2+0-9 4 
192 | 146667 15°8 | 42 33 | 5°44 oro | A3 Vn 12 3 6 
193 | 146775 16-0 | 28 10 | 7°68 +059 | Go V 31'1 o6 5 
194 | 147722 aa-s {| 26 3 Xe : Go IV O72 O'5 5 
195 | 147723 21°5 | 29 : 5°40 to-s8 Go IV 29 06 5 
196 | 148587 28°5 | 63 44 | 7°36 0°57 | Go V 40°O+1°1 4 
197 | 149324 36-0 | 77 2 4°20 + 1:07 | Ko III 30°9_0'5 3 
198 | 150248 38-2 | 45 16 | 7-05 +063 | G3 V 65°4 oO'9 4 
199 | 150437 39°0 | 29 02 | 7°86 +068 | G2 V 12°5 08 4 
200 | 151337 45°3 | 47 38 | 7°38 | +091 | Ko V + 2974 o8 4 
201 | 152334 51°1 217 | 3°67 1°38 | Ks III 17°9t0'5 5 
202 | 153072 55°5 | 37 33 | 605 | +o-19 | A3 V - aS 7 
203 | 152923 55°6 | 59 15 | 7°12 +o-45 | Fo V 13°4 1:0 4 
204 | 153075 56°3 | 57 13 | 7:01 0°56 | Go V 98°7 o'9 5 
205 | 153950 | 17 00°9 | 43 14] 7°41 0°55 | G2 IV-V + 30°7 o'9 6 
206 | 154577 05°7 | 60 40 | 7°40 089 | Ko V 88+ 1-0 4 
207 | 155886 12°3 | 26 32 | 4°38 0°84 | Ko V 06 oO4 5 
208 | 156026 132 | 26 29 | 6°32 116 | Ks V + 38 oF | 
2 156274 15°3 | 46 35 | 5°48 080 | G8 V + 24°5 oO 5 
210 | 156751 18-5 | 58 25 | 6-77 025 | As V 6 I 4 
211 | 157750 23°3 | 32 56 | 8-02 | +066 | G2 V 18-6+ 1-2 5 
212 | 157919 242 | 29 49 | 4°30 | +039 | Fs IV r 36:9 04 | 4 
213 | 158311 28-1 | 62 12 | 7°50 -o:g0 | G8 IV - 24°6 06 4 
214 | 159809 35'2 | 45 44] 7°44 +102 | Kr IV + 68 1:0 4 
215 | 159868 35°4 | 43 07 | 7°23 to-71 | Gs V 23'8 06 6 
216 | 160691 40'2 | 51 49] 5‘10 to-7o | G5 V g'7+t0'5 4 
217 | 162396 49°3 | 41 59 | 6°17 to-52 | F8 V — 14°4 08 4 
218 | 165189 | 18 03:2 | 43 26] 4:92 +o:23 | As V 6 2 8 
219 | 165271 03°7 | 46 54 | 7°65 +o-65 | Gs IV — 88-1 o8 4 
220 | 165696 05°8 | 45 57 | 7°33 -o-49 | F8 V — 29°3 o8 4 
221 | 165753 06:1 | 44 57 | 7:02 +1r1 | Ko IV — o7+to0'8 5 
222 | 166348 08-7 | 43 27 | 8:37 | +1°31 | K7 V - 2°70 06 4 
223 | 167576 13°8 | 27 44 | 6°66 +127 | K3 Ill —- 104 I°O 6 
2 167665 14°2 | 28 18 | 6°35 +052 | Go V + 69 O'5 8 
225 | 169233 21°8 | 30 47 | 5°58 -1-16 | Ko HII-IV | — 18-7 0-4 6 
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David S. Evans, A. Menzies and R. H. Stoy Vol. 117 
TABLE II (cont.) 
(1950) 
No. HD R.A. S. Dec V B-I Spec. Vel. (km/s) Pi. 
h m 7 
226 | 164461 | 18 25-8 | 87 39 | 5-29 1°30 | K3 Il 33°3+0°6 5 
227 | 171627 33°9 | 28 33 | 6-77 | +097 | Ki V + 247 09 | 4 
228 | 173339 43°2 | 43 51 | 7°39 | +027 | A7 V ol ne 4 
229 | 174153 47°7 | 44 32 | 7°54 | +0°52 | Go V 74 09 | 4 
230 | 175329 54°2 | 60 16 | 5:15 +41°37 | Ki III-IV 179°6 o8 5 
231 | 175856 55°83 | 44 12 | 8:57 +o’52 | Fo V + a3te7 5 
232 | 176687 59°4 | 29 57 | 2°59 +o'ro | A2 III + 28:2 o8 8 
233 | 176578 59°5 | 47 07 | 6°85 +o'97 | Ko IV —-— 71 O08 4 
eee ; a ; F8 V 50°8 03 4 
235 | 177474 | 19 03°0 | 37 08 | 4:20 °'5° | bg vy — §3°2 03 4 
236 | 179140 10°5 | 58 05 | 7-21 0°63 | G2 V 30°3+1°0 4 
237 | 181428 18°5 | 29 42 | 7°09 0°56 | Go IV 36°4 o-9 4 
238 | 181544 18°9 | 29 37 | 7°08 +057 | Go IV 36-4 06 5 
239 | 181773 21°2 | 62 17 | 7°59 0°46 | Fs IV 32°90 «&I'! 8 
240 | 182156 21°4 | 30 54 | 7°63 +117 | Ko IV 14°8 0-9 4 
241 | 183216 26°5 | 30 54] 7-12 +o'59 | G2 V 42'8+0°6 5 
242 | 183312 27°O | 32 12 | 6°56 +039 | Fs IV-V 8-3. 0-7 4 
243 | 183877 29°6 | 28 07 | 7-14 0°68 | Gs IV 40°0 O08 6 
244 | 186651 43°9 | 43 28 | 7-11 +055 | Go V 93° 1rd 5 
245 | 186975 45°9 | 45 52 | 7°25 | +110 | Ko IV — 244 04 | 4 
246 | 188555 54'2 | 45 58 | 8-60 + 0°39 5 IV 43°94 t-2 5 
247 | 188228 54°9 | 73 03 | 3°93 oor | Ao V 4 ry 6 
248 | 188815 55°4 | 46 14 | 7°46 0-45 | Fo V 6-4 08 5 
249 | 189198 57°3 | 45 15 | 5°80 0°28 7 Ill 6 1°5 4 
250 | 189247 57°4 | 44 07 | 7°66 0-42 | Fs IV 22°¢ -5°0 5 
251 | 189567 | 20 00°6 67 27 | 6:08 0°63 | G2 V 11°7+0°6 4 
252 | 190333 02°8 | 43 21 | 9:20 +066 | G2 V 58:2 1°0 5 
253 | 189899 03°2 | 74 22 | 7-61 +o50 | F8 IV-V 2-9 +1°2 4 
254 | 190248 03°8 | 66 19 | 3°54 + O75 38 V 21'2 O-2 4 
255 | 190779 O5*! | 46 27 | 8-23: | + 0-42: Fs V 43°0 o'8 5 
256 | 191408 °7°9 | 36 14 | 5:30 0°87 | K3 V 128:9+0°3 6 
257 | 191584 O9'O | 42 56 | 6-21 +125 | Ke Ill 08 os 5 
258 | 191849 10'2 | 45 19 | 7:96 1:44 | Mo V 20°8 o8 5 
259 | 192310 12°2 | 27 11 | 5°72 +o-91 | Ko V 54°2 073 5 
| 260 | 194433 23°7 | 37 34 | 6°23 +o-98 | K2 IV-V 25°2 04 4 
261 | 194640 24°6 | 31 02 | 6°59 +0o'73 | Gs V o'ot+0°7 6 
262 | 196081 32°9 | 26 57 | 7-19 | +0-41 Fs IV — 26°3 08 7 
263 | 196067 35°8 | 75 32 G1 V 12°5 08 4 
264 | 196068 35°9 | 75 31 om | +ot Gs V — 11°23 «1S 5 
265 | 196051 35°9 | 76 22 | 6-00 +045 | Fs [Vn — 36 2 7 
266 | 196531 36°0 | 28 36 | 7-94 | 4 0°53 | F8 V — 32°6+08 4 
267 | 196917 38°3 | 31 47 | 5:76 | +1-55 | Mo III — 97°3 9°9 7 
268 | 197214 40°2 | 29 36 | 6-95 +067 | Gs V — 198 oF 4 
269 | 197484 42°2 | 43 08 | 8-86 +063 | G2 V + 211 o'9 4 
270 | 197416 42°6 | 60 28 | 8-15 +o-48 | F8 V — 36 o8 4 
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Fundamental data for southern stars 


TaBLe II (cont.) 








(1950) 

No. HD R.A. S. Dec 

h m , 
271 | 197900 | 20 44°9 | 44 23 
272 | 197937 45°I | 44 10 
273 | 198828 51°4 | 46 47 
274 | 199288 54°4 | 44 19 
275 | 200011 589 | 43 12 
276 | 200026 59°O | 43 12 
277 | 202560 | 21 14°3 | 39 04 
278 | 202457 14°8 | 61 33 
279 | 203448 20°I | 31 02 
280 | 203608 22°3 | 65 36 
281 | 203985 23°8 | 45 02 
282 | 205067 30°6 | 28 07 
283 | 205153 31°2 | 28 07 
284 | 205390 33°3 | 51 04 
285 | 206341 39°2. | 27 5g 
286 | 206395 39°8 | 43 43 
287 | 206868 2-7 | 28 s7 
288 | 206948 43°7 | 46 37 
289 | 207155 44°8 | 31 08 
290 | 207129 45°O | 47 32 
291 | 207964 51°5 | 62 07 
292 | 208215 52°8 | 47 10 
293 | 208812 56°8 | 43 43 
294 | 209100 59°6 57 00 
295 | 209742 | 22 03°6 | 45 38 
296 | 209952 Ost | 47 1 
297 | 210193 06°6 | 41 2 
298 | 210441 084 | 44 05 
299 | 211415 15°O | 53 52 
300 | 211998 20°4 | 72 30 
301 | 213042 26°4 | 30 16 
302 | 213135 27°0 | 27 22 
303 | 214094 33°8 | 43 44 
304 | 214308 3571 | 46 58 
305 | 214385 35°4 | 27 42 
306 | 214690 470 | 30 55 
3°07 | 214749 37°9 | 29 56 
308 | 214759 381 2 15 
309 | 214953 39°7 | 47 28 
310 | 216042 47°2 | 33 04 
31r | 216054 47°3 | 41 45 
312 | 216743 53°I | 42 49 
313 | 216803 53°6 | 31 50 
314 | 216956 54°9 | 29 53 
315 | 216989 55°2 | 45 26 
























V B-l Spec. 
6°46 +4118 | Kr IV 
5*10 +0735 | Fi V 
7°39 | +0°53 | F8 V 
6-50 | +058 | Go V 
6°63 +068 | G3 IV 
6°89 +o'97 | Ko IV 
6-68 +1°42 | Mo V 
6°58 +o-7o0 | G5 V 
7°86 +056 | Go IV 
4°23 +047 | F8 V 
7°45: | +0°94:| Ko V 
7°61 +065 | G2 V 
8-2 +054 | Go IV 
7°37 0°89 | K2 V 
7°68 +108 | Ko IV 
6:69 +053 | Go IV 
7°69 0°39 | F2 V 
7°56 118 | K3 III 
5°03 0°06 | Az IIIn 
5°60 0:59 | G2 V 
5°91 0o'40 | Fo IV 
6°55 o45 | F5 V 
8-17 0°47 | F8 IV-V 
4°73 1:06 5 V 
8-46 o'86 | K2 V 
“te fF o'14 | Bs V 
7°37 +o-65 | Gs V 
6°63 +1-00 | G8 IV 
5°40 +o'59 | G1 V 
5°29 0°64 | Go V 
7°68 tiro | Ks V 
5°99 | +033 | Fo V 
6°77 Lo-52 | F6 V 
7°70 +o-41 | Fs IV 
7°92 +o62 | G2 V 
5°89: | +1°31 | K3 III 
7°86 +114 | Ks V 
7°41 +o-80 | G8 V 
5°99 t+o56 | Gi V 
6°36 | +o-30 | F2 IV 
7°80 +o-73 | Gs V 
7°28 to-18 | Ar Vn 
6°47 +113 | Ks V 
1°18 to:12 | A3 V 
7°73 0°36 | Fo V 








Vel. (km/s) 
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TABLE II (cont.) 

(1950) = ioe ner 
No. HD R.A.  S. Dec V B-I Spec. Vel. (mk/s) Pl. 
h m i . ; " 
316 | 217595 | 22 59°6 | 45 34 | 7°22 | +043 | Fs V 20'5+0°3 | 7 
317,| 217766 | 23 00°7 | 43 21 | 7°78 +053 | F8 V 34°0 I'l 5 
318 | 217816 o1'r | 46 26 | 8°13 t+o50 | Fé IV-V t 5:0 O8 5 
319 | 217987 02°7 | 36 o9 | 7°36 | +1°50 | M2 ©: + 87 o8 5 
320 | 218227 o4't | 43 47 | 4:29 | +0°44 | Fo IV 4°7. 1°0 4 
321 | 219249 12°2 | §7 00 | 7°97 +069 | G5 V — 58+0°7 4 
322 | 219409 13'1 | 30 07 | 6°53 +109 | Ki III — 47°7 06 5 
323 | 219482 14°0 | 62 16 | 5°67 | +0°49 | F8 ¥ 070 0% 4 
324 | 219709 15°7 | 58 35 | 7°33 +0°69 G2 V + 20°38 o8 + 
325 | 220881 24°9 | 27 33 | 7°46 +o:28 | Fo IIl — 18°77 13 5 
326 | 221818 32°9 | 47 13 | 8°57 +o'78 | G8 V 5‘ot 1°0 5 
327 | 221839 33°0 | 27 46 | 6°68 +055 | Go IV 18 03 5 
328 | 222237 36-7 | 72 59 | 7°09 +1or | K3 ¥ 7o'1 o8 6 
329 | 222412 37°8 | 26 29 | 7°60 +043 | Fo IV — 4 1 7 
330 | 222433 38°0 | 32 21 | §°33 +097 | Ko III + 156 073 5 
: 331 | 222480 38°5 | 32 21 | 7°13 +066 | Gs IV 21‘5to9 5 
332 } 222508 38°38 | 41 51 | 7°82 | +046 | F7 V 10 2 5 
333 | 222803 41°4 | 45 22 | 611 | +0-99 | G8 IV 30°8 06 4 
334 | 223065 43°9 | 41 51 | 7°0 +or2 Az V — 33 2 9 
335 | 223633 48°8 | 42 39 | 7°57 | + 0°44 | Fs IV-V a ee 9 
336 | 224296 54°1 | 42 28 | 7:90 | +0-41 | Fs5 IV-V + O6u+1°1 6 
337 | 224360 54°7 | 46 23 | 7°72 | +044 | Fs V a 6 
338 | 224750 57°8 | 44 34 | 6:30 | +0°75 | G3 IV 370 O'S 5 
339 | 224889 59°1 | 77 20 | 4°79 +129 | K2 IIl 22°4 03 6 
R2 1581 | 00 17°5 | 65 10 | 4°23 | +0°56 | G2 V 9°5+0°4 4 
R4 | 20794 | 03 17°9 | 43 16 | 4:24 | +0-70 | G5 V 87°0 03 4 
R6 | roro2zr | 11 34°7 | 61 00 | 515 | +1°13 | Ki III 3°3. O'S 4 
R7 | 114837 | 13 111 | 58 50 | 4:92 | +0°46 | F8 V 64:8 1-0 5 
































Notes 


(We would like to express our indebtedness and thanks to the Union Astronomer, Dr 
W. S. Finsen, who has kindly provided up-to-date information from the Union Observatory 
card catalogue for many of the double stars. The magnitudes and dynamical parallaxes 
(pq) quoted in these Notes have been taken unchanged from several sources and no attempt 


has been made to render them homogeneous.) 
4 PBHyi. 


5 By Tuc. Bi. 4:9:5°9, 41°39; a=0"-48, pg=o"-024. Mean date of rad. vel. 
obs. 1954.22. (8,+f, Tuc. Bi. 4:48 : 4°52, 27°°1, 170°). 


6 ADS 520 Bi. 6-3 : 6-4, 25-0; a=0"'67, pg=o"-063. 


obs. 1953.98. Single lines only seen in spectra. 


10 ~—C Difficult spectrum, broad lines. 


11 The present plates do not support the designation as a sp. bi. in the Mount 


Wilson Cat. (No. 500). 


13, 14 Bi. 7°6 : 7-6, 5”-0, 128°. No certain motion. 


magnitude and colour. 





C: 11 mag., 100”, 41°. 


Mean date of rad. vel. 


Joint 
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16 
17 
18 
19 
22 


23, 24 


26 
37 
40 


41 
44 
46 
47 


64 


638 
69 
70 
71 
79 
So 
86 
89 
91 


g2 
96 
97 
99 
101 
102 





Bi. 7°9 : 8-2, 0-3, 314 (1945). Joint spectrum. 

Bi. 7°2 : 10, 2”-6, 229° (1942). py=o"-o19. 

Difficult spectrum, broad lines. 

Eggen gives mag. as 7°74. Knox-Shaw from southern Harv. vis. obs. gives 8-21. 

Triple 7°8 : 7-9: 11°2; ABa 0”°17,4°6 y ; AB, Ca=1"'58?, 144 y?. See O. J. 
Eggen (P.A.S.P., 64, 230, 1952) py=0o"’051. Mag. prob. var. A fourth star 
D of mag. 7 is such that AB, D is 140”, 20°. 

p Eri. Bi. 5°8: 5:9, 244 y, a=8"-36; py=o"'165. See J. L. Dessy (Pub. La 
Plata, 20, No. 3, 1949); also G. B. van Albada (Bosscha Cont. No. 5, 1956). 
‘The suggestion that these stars are peculiar does not seem well founded. The 
secondary and primary (Nos. 23 and 24 respectively) are Nos. 770 and 771 in the 
list of parallax stars given in Cape Annals, XV, 1952. The secondary is 
0-16 mag. fainter than the primary, giving, on the assumption of identical 
colours for the components, visual magnitudes of 5-91 and 5°75. The joint 
colour is consistent with the spectral types. The deduced photometric 
parallaxes are 0”-129 and 0”-120 for secondary and primary, while the observed 
Cape relative parallaxes are 0”:138 and 0”:129. 

is Phe. Mag. var. through o™-12. H. Spencer-Jones (Cape Ann., 10, Pt. 8, 1928) 
suspected var. vel. Present rad. vel obs. all close to mean date 1952°58. 

mw For. 

€ For. 

x For. ADS 2402 Bi. 4:0 : 6°5, 154°5 y, a=2"'70, py=0"'053. Companion 
may be var.; Spectra of bright star only. Joint mag. and col. 

C! Ret. 

¢* Ret. Common motion with 41. 

Mag. and col. measures impeded by neighbouring star of differing p.m. 

Bi. 8-2 : 8-3, 0”:5, 114° (1955). 

Bi. 6-7 : 7°3, 19°4 y, a=0"'24, p, —0"'024. Vel. refers to combined light, mean 
date 1953-92. 

Vel. possibly var. 

B Ret. Mount Wilson Cat. (No. 2114) designates as a sp. bi. Mean date of 
present obs. 1953°81; var. of vel. not confirmed. 

y Hyi. 

a Cae. Bi. 4°5 : 12, 6”-6, 121° (1933). 

Bi. 7°1 : 7°3, 2°°1, 91°; pq=0"'037. Parabolic orbit (Union Obs. Circ. 95). 
Spectra of combined light shows single lines. 

Interferometric bi., 0”:228, 40°°0 (1956°16). See W. S. Finsen (Obs., 72, 125, 
1952). Lines very broad, only one spectrum seen. 

Bi. 7:0 : 12, 0”°9, 95°. 

¢ Dor. 

H and K lines in emission, + 17 km/s. 

Triple 5-5 : 11°7 : 12°7, A, BC 25”°9, 257°; BC 0”°6, go’. 

Vel. probably var. 

Joint mag and col. with GC 7049. Bi. 7°4 : 9°5, 5”°1, 67° (1941), Pg =0"'037- 

x Men. 

Bi. 8-0 : 10°5, 2°-4, 186° (1951). 

H and K lines in emission, ~ 16 km/s. 

Bi. 7-1 : 7°9, 0”:16, 115° (1956). Hand K lines in emission, + 56 km/s. 

The faint companion of g Pup I, 22”,common motion. Rad. vel. obs. of o Pup I 
(D. S. Evans, WN ASSA, 16, 4, 1957) confirm, with only slight modifications, 
the work of R. E. Wilson (Lick Obs. Bull., 9, 117, 1918). The stars are too 
close for the photometry of the faint companion. 

Vel. possibly var. 

Q Pup. 

Bi. 5-1 : 8:1, 3°°3, 275 (1954), Pg =O” 021. 

a Cha. 

Difficult spectrum, broad lines. 

H and K lines in emission, ~ 11 km/s. 

5 Pyx. 
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TABLE II (cont.) 
(1950) 
No. HD R.A. S.Dec| V B-V Spec. Vel. (mk/s) Pl. 
h m f 
316 | 217595 | 22 59°6 | 45 34 | 7°22 | +0°43 | Fs V + 2o5to3 | 7 
317. | 217766 | 23 00°7 | 43 21 | 7°78 +053 | F8 V 34°0 I'l 5 
318 | 217816 ort | 46 26 | 8-13 to:50 | Fé IV-V + 5:0 O8 5 
319 | 217987 02'7 | 36 09 | 7°36 | +1°50 | M2 V + 87 o8 5 
320 | 218227 04°! | 43 47 | 4°29 +044 | Fé IV 4°7 10 4 
321 | 219249 12°2 | 57 00 | 7°97 | +0°69 | G5 V — §8+0°7 4 
322 | 219409 13°1 | 30 07 | 6°53 +109 | Kr III —~ 47°7. 06 5 
323 | 219482 14°0 | 62 16 | 5:67 +o-49 | F8 V oO 0% 4 
324 | 219709 15°7 | 58 35 | 7°33 | +0°69 | G2 V + 208 o8 | 7 
325 | 220881 24°9 | 27 33 | 7°46 to:28 | Fo Ill — 18:7 173 5 
326 | 221818 32°9 | 47 13 | 8°57 +o'78 | G8 V - §'0+1°0 5 
327 | 221839 33°0 | 27 46 | 6°68 | +0°55 | Go IV + 38 of 5 
328 | 222237 36°7 | 72 59 | 7°09 +ror | K3 V + Fo1 o8 6 
329 | 222412 37°8 | 26 29 | 7°60 | +043 | Fo IV — 4 I 7 
330 | 222433 38°0 | 32 21 | §°33 +097 | Ko III + 1§6 0% 5 
331 | 222480 38°5 | 32 21 | 7°13 | +0°66 | Gs IV + 21'§5to'9 5 
332 }| 222508 38°8 | 41 51 | 7°82 | +046 | F7 V + 10 2 5 
333 | 222803 41°4 | 45 22 | 6-11 | +0-99 | G8 IV — 308 o6 4 
334 | 223065 43°9 | 41 51 | 770 | +02 | Az V ae 9 
335 | 223633 48°8 | 42 39 | 7°57 | + 0°44 | Fs IV-V t+ 22 1 9 
336 | 224296 54°1 | 42 28 | 7-90 | +041 | Fs IV-V + O61+1°1 6 
337 | 224360 54°7 | 46 23 | 7°72 | +044 | Fs V ~~ 2 6 
338 | 224750 57°8 | 44 34 | 6:30 | +075 | G3 IV + 3° oO'5 5 
339 | 224889 59°! | 77 20 | 4°79 +1°:29 | Ke IIl + 22°4 O73 6 
R2 1581 | 00 17°5 | 65 10 | 4:23 | +0°56 | G2 V + g§t04 4 
R4 | 20794 | 03 17°9 | 43 16 | 4:24 | +0-70 | G5 V * 87'0 03 4 
R6 | 1roro2zr | 11 34°7 | 61 00 | 5°15 | +1°13 | Kr III +t 943 oO 4 
R7 | 114837 | 13 11'1 | 58 50 | 4:92 | +0°46 | F8 V — 64°83 1'o 5 



































Notes 


(We would like to express our indebtedness and thanks to the Union Astronomer, Dr 
W. S. Finsen, who has kindly provided up-to-date information from the Union Observatory 
card catalogue for many of the double stars. The magnitudes and dynamical parallaxes 
(pq) quoted in these Notes have been taken unchanged from several sources and no attempt 
has been made to render them homogeneous.) 
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10 
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B Hyi. 

B, Tuc. Bi. 4°9:5°9, 41°39; a=07-48, pyg=o"-024. Mean date of rad. vel. 
obs. 1954.22. (8,+ 8, Tuc. Bi. 4°48 : 4°52, 27”*1, 170°). 

ADS 520 Bi. 6-3 : 6-4, 25-0 y; a=0"-67, py=0"-063. Mean date of rad. vel. 
obs. 1953.98. Single lines only seen in spectra. 

Difficult spectrum, broad lines. 

The present plates do not support the designation as a sp. bi. in the Mount 
Wilson Cat. (No. 500). 


13, 14 Bi. 7°6 : 7°6, 5”-0, 128°. No certain motion. C: 11 mag., 100", 41°. Joint 


magnitude and colour. 
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Fundamental data for southern stars 


Bi. 7°9 : 8:2, 0°°3, 314 (1945). Joint spectrum. 

Bi. 7-2 : 10, 2”°6, 229° (1942). py=o”-org. 

Difficult spectrum, broad lines. 

Eggen gives mag. as 7°74. Knox-Shaw from southern Harv. vis. obs. gives 8-21. 

Triple 7°8 : 7-9: 11:2; ABa=o"'17, 4-6 y ; AB, Ca=1"-58?, 144 y?. See O.J. 
Eggen (P.A.S.P., 64, 230, 1952) py=0":051. Mag. prob. var. A fourth star 
D of mag. 7 is such that AB, D is 140”, 20°. 

p En. Bi. 5°8:5°9, 244, a=8"-36; py=o"'165. See J. L. Dessy (Pub. La 
Plata, 20, No. 3, 1949); also G. B. van Albada (Bosscha Cont. No. 5, 1956). 
‘The suggestion that these stars are peculiar does not seem well founded. The 
secondary and primary (Nos. 23 and 24 respectively) are Nos. 770 and 771 in the 
list of parallax stars given in Cape Annals, XV, 1952. The secondary is 
0°16 mag. fainter than the primary, giving, on the assumption of identical 
colours for the components, visual magnitudes of 5-91 and 5°75. ‘The joint 
colour is consistent with the spectral types. The deduced photometric 
parallaxes are 0-129 and 0”-129 for secondary and primary, while the observed 
Cape relative parallaxes are 0”-138 and 0”:129. 

is Phe. Mag. var. through o™-12. H. Spencer-Jones (Cape Ann., 10, Pt. 8, 1928) 
suspected var. vel. Present rad. vel obs. all close to mean date 1952°58. 

aw For. 

€ For. 

« For. ADS 2402 Bi. 4:0 : 6°5, 154°5 y, a=2"°70, py=0"'053. Companion 
may be var.; Spectra of bright star only. Joint mag. and col. 

C! Ret. 

¢* Ret. Common motion with 41. 

Mag. and col. measures impeded by neighbouring star of differing p.m. 

Bi. 8-2 : 8-3, 0”°5, 114° (1955). 

Bi. 6-7 : 7°3, 19°4 y, @=0"'24, py —0"'024. Vel. refers to combined light, mean 
date 1953°92. 

Vel. possibly var. 

B Ret. Mount Wilson Cat. (No. 2114) designates as a sp. bi. Mean date of 
present obs. 1953°81; var. of vel. not confirmed. 

y Hyi. 

a Cae. Bi. 4°5 : 12, 6-6, 121° (1933). 

Bi. 771: 7°3, 2°°1, 91°; pq=0"'037. Parabolic orbit (Union Obs. Circ. 95). 
Spectra of combined light shows single lines. 

Interferometric bi., 0”:228, 40°-0 (1956°16). See W. S. Finsen (Obs., 72, 125, 
1952). Lines very broad, only one spectrum seen. 

Bi. 7°0 : 12, 0”°9, 95°. 

¢ Dor. 

H and K lines in emission, + 17 km/s. 

Triple 5:5 : 11°7 : 12°7, A, BC 25-9, 257°; BC 0”°6, go’. 

Vel. probably var. 

Joint mag and col. with GC 7049. Bi. 7°4 : 9°5, 5”°1, 67° (1941), Pg =0"'037.- 

x Men. 

Bi. 8-0 : 10°5, 2°°4, 186° (1951). 

H and K lines in emission, — 16 km/s. 

Bi. 7°1 : 7°9, 0”16, 115° (1956). Hand K lines in emission, +56 km/s. 

The faint companion of g Pup I, 22”, common motion. Rad. vel. obs. of o Pup I 
(D. S. Evans, MNASSA, 16, 4, 1957) confirm, with only slight modifications, 
the work of R. E. Wilson (Lick Obs. Bull., 9, 117, 1918). The stars are too 
close for the photometry of the faint companion. 

Vel. possibly var. 

Q Pup. 

Bi. 5-1 : 8-1, 3°°3, 275° (1954), Pg=0”'02I. 

a Cha. 

Difficult spectrum, broad lines. 

H and K lines in emission, + 11 km/s. 

5 Pyx. 
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Bi. 5:2: 11, 24”°1, 338° (1929). 

Bi. 6-0 : 10, 8-7, 328° (1933). Difficult spectrum, lines broad and faint. 

This star has a common p.m. with HD 82207 from which it is distant 108"-4 in 

a. 97°. 

ys Vel Bi. 3°9:4°8, 34°1¥, @=0"'92; py=0"'061. Companion reported as 
variable 4™-5-5™-1. Lines broad. Only one spectrum recorded but its 
appearance is suspected of variation with lines narrower in 1953. Cape 
photometric obs. of massed light give no sign of var. 

Bi. 5°4 : 7°1, 0”°5, 330° (1956). 

Bi. 8-3 : 10°8, 5*1, 128° (1947). 

Bi. 6°5 : 7°7, 211 y, a@==1"°0, pg=0"'019. Mean date of vel. obs. 1954°70. 

Bi. 7°6 : 11, 1”°9, 156° (1935). 

x Vel. A multiple star. AB 4-4: 6-6, 517-6, 105° (1938); AD 4-4: 14, 187-4, 
86° (1880); AE 4°4 : 13, 33°71, 225° (1880); BC 6-6: 11, 21”-0, 174° (1934). 

ut Vel. Bi. 3:0 : 6°5, 0”°65, 90° (1943). This pair has closed in and the companion 
has been invisible since 1945. Period unknown; pg=0”'022. Vel. of brighter 
star only, possibly var. 

m Cen. Bi. 4:7: 5°5, 42°8 y, a=07'275, py=O"oll. 

Bi. 7°6 : 8-6, 265 y, a=4"°01, py=0"'105. 

Companion, 15™, 17”, which, according to a note in the Yale ‘‘ General Catalogue 
of Parallaxes ’’, has the same proper motion. 

Bi. 5°8 : 13, 6-7, 163° (1949). Possibly an optical double. Mag. and col. 
probably var. 

H and K lines in emission, -; 16 km/s. 

Bi. 4°1 : 13, 25”, 320° (1929). 

Difficult spectrum, broad lines. 

Mount Wilson Cat. (No. 7208) designation as a sp. bi. seems doubtful. Eggen 
gives mag. as 7°23. Knox-Shaw from southern Harv. vis. obs. gives 6°94. 

y Cen. Bi. 3:1 : 3°1, 84°5 y, a=0"'930, py=0"'024. All obs. refer to combined 
light. Mean date of vel. obs. 1952-49. 

«Cru. Bi. 4-7: 9°5, 26’:9, 22° (1942). 

m Cen. 

d Cen. Bi. 4-6: 4:8, 62-6 y, a=0"°155, pg=0"'003. Obs. refer to combined 
light. 

Difficult spectrum. 

Bi. 7°4.: 7°5, 30°°5, 131° (1932). A wide pair making measurement of mag. and 
col. difficult. 

51 Hy:. 

a Lup. Tri. AB 6-1 : 6-1, 0”-142, 41° (1955-6); AB, C 5-4: 7°9, 19”°0, 24° (1933)- 
Mean date of vel. obs. is 1954:10. Lines suspected of duplicity in 1953 Feb. 

a Cir. Bi. 3-2 : 8-6, 16”, 232° (1943). 

Difficult spectrum, broad lines. 

Forms a bi. with HD 134330 7°0:7°5, 50”°6, 22° (1902). Common p.m. 

B Cir. 

Vo Lup. 

w Lup. Bi. 4-3 : 11, 11”-8, 28° (1933). Present obs. do not support suggestion 
of var. vel. by H. Spencer Jones (Cape Ann., 10, Pt. 8, 1928). 

g Lup. Difficult spectrum, broad lines. 

Difficult spectrum with diffuse lines. This star is underluminous. See J. 
Jackson (MNASSA, 8, 29, 1949). 

Bi. 8-3 : 8-8, 0”-15, 324° (1935). 

Bi. 6-0 : 13°5, 15”, 133° (1953). 

BTrA. A 14™ star which s.p. 8 Tri. by 155” possibly shares its p.m. 

Nor. Tri. AB 5-6:5°7, 26-4, a=0"-356, py=0"-028; AB, C 46:8, 
10°°8,245°. Difficult spectrum, broad lines. Vel. refers to mean date 1953°12- 

8 TrA. Bi. 4:0: 11, 25”, 105° (1892). 

d Sco. 

A Nor. Bi. 6-0 : 6-9, 0”-4, 120° (1956). 

Bi. 5-94 : 6°59, 5-7, 353° (1935), Py=0"'036. Joint mag. and col. 
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8 Aps. 

f Sco. 

Bs. 6-9 > 7-1; 31-2 ¥, a=0"-166, pq -O-'orr. Spectrum difficult with broad lines 
lines but no indication of var. vel. or double lines. Mean date of vel. obs. 
1953°73- 

Bi. 7-1 : 11, 26°°1, 12° (1918). 

36 Oph. Bi. 5:29 : 5-33, 4°°4, 166 (1955). ‘The mag. and col. refer to the joint 
light but the rad. vel. and spectral type to the north component only. Hand K 
lines in emission, — 3 km/s. 

H and K lines in emission, --2 km/s. Vel. possibly var. This star and 207 have 
a common motion. 

41 Ara. Multiple star. AB 5:5 : 8-7, 550 y, a=8"°8, py=o"'113. AC 5-5: 13, 
41°°8, 279° (1900); AD 5°5 : 8-5, 47”, 30° (1900). Mag. and col. refer to joint 
light, rad. vel. and spectrum to bright star only. 

Bi. 6-8 : 9-6, 9”-4, 283° (1917). 

Vel. possibly var. 


45 Oph. 

Bi. 7-4: 12, 1”°8, 246° (1948). 

ym Ara. 

Bi. 5°8:5°8, 214 ¥, a=1":208, Pp, —-0"'022. Difficult spectrum, broad lines. 


Vel. refers to the combined light, and is possibly, var. Mean date of rad. vel. 
obs. 1952°74. 

H and K lines in emission, — 3 km/s. 

x Oct. Spectra all rather weak. 

Extremely difficult spectrum, lines very broad. 

w Pav. 

C Ser. Bi. 3-4: 3°6, 20°8 vy, a=0"'520, py -0”-040. All data for combined light. 
Vel. not considered var. 

y CrA. Bi. 4°95 : 5°05, 119°3¥, @=2":07, py -0”'064. Vel. of components 
observed separately, mean date 1953:9. 

Bi. 7°3 : 9°8, 075, 195° (1942). 

Bi. 7°7 : 10°8, 2”-9, 53° (1937). 

€ Pav. 

Difficult spectrum, broad lines. 

Bi. 9°4.: 11, 1°°7, 299° (1934). 

5 Pav. 

With HD 190803 forms Bi. 8-2 : 9:2, 63”:1, 102°. This is probably an optical 
pair only as the p.m.s differ. Mag. and col. not sufficiently observed to deserve 
standard weight. 

Bi. §°3 : 11°5, 7°‘°I, 123° (1949). 

Mount Wilson Cat. (No. 12602) designation as sp. bi. not supported. 

Bi. 6-5 : 8-2, 0°25, 245° (1952), Pp, 0"'024. lll data refer to joint light. Mean 
date of vel. obs. 1953°69. 


263, 264 Bi. 7°08 : 7°60, 17"°2, 17° (1940). Joint mag. and col. 


“#, Oct. Difficult spectrum with broad lines. 
t Mic. Bi. 5:1 : 15, 4°°3, 271 (1932). Difficult spectrum, broad lines. 
Bi. 7°9 : 8:4, 0”°4, 138° (1952). All data for combined light. 


76 Bi. 6°51 : 6:94, 57°°4, 73° (1930). Common motion. 


H and K lines in emission, 23 km/s. 

y Pav. 

Mag. and col. not sufficiently observed to deserve standard weight. 

@PsA. Tri. AB 5-7: 5°8, 0°-179, 35° (1954); AB, C 5:0: 11, 35°°5, 339 
(1918). Difficult spectrum, broad lines. 

Bi. 6°6 : 9°7, 55”°0, 356 (1931). 

Bi. 6-5 : 6°7, 25°8 y, @=0"'237, Pg 0”-019. All data for combined light. Diffi- 
cult spectrum, broad lines. 

e Ind. H and K lines in emission, — 42 km/s. 

xy Gru. Bi. 1-7: 12, 28"*5, 149° (1947). Mag. suspected of variability. 
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Bi. 5°4 : 10°0, 3°°2, 29° (1948). 

vInd. Bi. 6-1 :6°2, <0 +115 (1953). Not seen as Bi. since 1933. Duplicity 
suspect. 

Mag. possibly variable. 

H and K lines in emission, — 5 km/s. 

Bi. 6-0 : 10°0, 7”°8, 129° (1943). 

Bi. 7°0:7°2, 27°2, @=0"'21, py -o”'015. All data refer to combined light. 
Mean date of vel. obs. 1954-07. 

Difficult spectrum, broad lines. 

This star has approximately the same parallax and proper motion as 314, but the 
rad. vel. differ by 10 km/s. 

a PsA., Fomalhaut. Standard star in MKK Atlas. 

Bi. 8-1 : 9°7, 3”°1, 107° (1935). 

H and K lines in emission, -- 6 km/s. 

8Gru. Multiple star. AB 4°5:7-0, 1°°3, 68° (1955); AC 4°5:8:5, 161”, 
293° (1904); CD 8-5 : 12, 95", 335 - pq (from AB)=0":027. All data refer 
to combined light of AB. Mean date of vel. obs. 1953:12. 

Bi. 7-3 : 9°5, 34°°3, 322° (1955). 





Broad lines. 
Bi. 6-9 : 8-7, 146 v, a~0"'656, pg = 0”'025. 
u Scl. 


Bi. 8-2 : 8-7, 0”-8, 192° (1956). All data refer to combined light. 

SX Phe This is the var. star with the shortest known period (80 min.). Maxima 
7™-1o-6™-54: Minima 7™-41-7™-53: Amplitude o™-8-o™-3. Absolute 
mag. too faint to be classed as RR Lyrae on present ideas. See J. Jackson 
(MNASSA, 8, 29, 1949), O. J. Eggen (P.A.S.P., 64, 31 and 305, 1952), Th. 
Walraven (BAN, 12, 57, 1953). 

Bi. 6°8 : 7-3, 0”-36, 291° (1954). 

8 Oct. 


82 Eri. 
Bi. 4°9 : 10°3, 2”°7, 341 (1928). 








Relation between spectral type and colour index 
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Fundamental data for southern stars 


















ee ETS ee 


























| 
Cape Johnson and Morgan 
‘Type 
C.. No. | Range C.I. No. | Range 
Fo Ill ~0'26 3 0°03 
Ko III 1‘05 5 O'14 + 1°00 10 o'10 
Kr Ill I°tl 5 o'10 I‘12 2 0°03 
K2 Ill + 1°23 3 0°13 I‘17 5 0°08 
K3 III 1°23 5 o°21 1°30 4 O'ls 
Ks Ill 1°34 4 0°10 1°52 2 o'o! 
Mo III 1°53 3 0°17 
F2 1V + 0°32 6 0°07 0°35 I 
Fs 1V o'4! 9 oO'12 
F6 IV 0°42 4 0:06 0°44 I 
Go IV -O°55 9 0°05 o'61 2 0°06 
G21V 0°62 3 0°02 
G3 IV 0°66 4 O14 
Gs IV o'75 13 0°30 
G8 IV 0°97 |' 6 O14 0°86 I 
Ko IV 1‘O7 12 0°29 0°92 I 
Fs IV-V o'4! 5 0°05 0°40 I 
F6 IV-V 0°48 3 0°06 
F8 IV-V ~ 0°48 i 0°03 0°54 2 ool 
A3 V O13 5 0°09 0°09 3 0°03 
As V 0°24 . 0°02 O'l4 4 0°03 
A7 V 0°28 3 0°02 o'18 3 ool 
Fo V 0°34 7 0:09 0°33 2 0°02 
F2V 0°36 5 0:06 ~0°38 3 0°03 
F3 V 0°42 3 0°04 + 0°38 I 
Fs V 0°45 rl 0°05 0°44 3 0'03 
-F6 V 0°47 7 0°07 0°46 4 0°04 
F7 V 0°49 3 0°06 "50 7 0°04 
F8 V 0"50 21 oo fe) 0°53 3 0°04 
~GoV 0°56 18 o'14 0:60 7 0°05 
“GiV 0°57 6 0:06 0°61 I 
G2 V +0°63 16 0°13 0°64 2 o'ol 
Gs V Lo-71 21 0°16 0°68 4 0°05 
G8 V 0:78 5 0°05 Lo-74 2 0°03 
Ko V 088 15 O'15 o'82 2 0°07 
K2V 0°93 5 o'16 0°88 2 0°02 
K3 V 0°94 6 0°20 1‘O1 | I 
Ks V I‘12 II 0°13 118 2 0°03 
| 
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TABLE IV 
Relation between spectral type and absolute visual magnitude 
M, 
Type M, No. 
Gliese Allen 
IV 3°7 4 31 
Vv 4°4 5 
Vv 4°7 6 4°5 4°6 
Vv 4°9 5 
Vv 48 8 4°7 
Vv 5°2 II 5"0 5°3 
Vv 5°4 4 5°5 
Vv 6:2 7 5°8 6-2 
Vv 6°3 4 6-3 
V 6°8 5 6°5 
Vv 7% 10 7s 7°6 
Note 


Following Johnson and Morgan (3), the M and K magnitudes have been adjusted by 
—om-y, 

































TABLE V 


Fundamental data for southern stars 












Hs 


























+0°44 

Type Fs V 
| 21 | +0°45 
64 | +0°47 
81 | +0°45 
97 | 10°44 
172 | +0°46 
| 255 | +0°42: 
| 292 | +0745 
| 316 | +0°43 
337 | +0°44 

Type F6 V 
66 | +0°45 
168 | +0°47 


203 | +0°45 


231 +o'5§2 
248 | +0°45 
303 | +0°52 
Type F7 V 
17 | 70°52 
158 | +o"50 
Type F8 V 
R7 | +0°46 
7 | +0°54 
9 | +0°47 
34 | +0°55 
45 | 0°47 
53 | +o'52 
56 | +0°45 
63 | +-o'50 
65 | +o'51 
108 | +-0°48 
154 | +051 
188 | +047 
217 | +0°52 
220 | +0°49 
266 | +0°53 
270 | +0°48 
280 | +0°47 
317 | +0°53 
Fo°49 




















4°08 


4°24 
4°08 
4°00 
4°16 
3°82 
4°08 
3°91 


4°08 
4°24 
4°08 
4°58 
4°08 
4°58 


4°58 
4°47 


4°16 
4°68 
4°24 
4°73 
4°24 
4°58 
4°08 
4°47 
4°53 
4°32 
4°53 
4°24 
4°58 
4°40 
4°63 
4°32 
4°24 
4°63 














O'O17 
0'022 
0062 
0062 
0°033 
0'013 
0°032 
0°022 
o'o18 


0'021 
O'O17 
0°025 
o'016 
0'021 








O'O4I 
0°'067 
0'026 


0°036 | 0'033 


0033 
0'026 


0°070 
0'017 
0°021 
0°021 
0°042 
0°027 
0°039 
0°044 
0°093 
0030 
0'036 
0°033 
0048 
0°026 
0°022 
0'O17 
0*100 
0°023 








0°047 


0'058 


0'026 


0'078 


oO'llit 








| 





0°172 
-0°027 
0'030 
—0'200 
~O°194 
0061 
-0°127 
—o'o18 
0109 


O°OI5 
+0038 
+0021 
— 0°34 
+o°o81 
+0°245 


~0'°165 
—0'057 


—0'264 
—0°099 
+0098 
—0O'OoI7 
—0°023 
+-0°095 
+0°048 
—0°057 
—0'042 
—o0:o89 
—0°147 
—0°093 
+ 0°134 
+-0°O31 
+0:088 
+o'110 
+0°'096 
—0°003 











+0'020 
—OoO'l21 
+0-386 
+0°241 
—0°133 

-0°OQI 
—0'063 
—0'008 
+-O°O14 


+ 07164 
—O'103 
—o-108 
+0°022 
+-0°020 


+0078 
—0'067 


—o-169 
—o-096 
—0°046 
—o°280 
—o'187 
+0083 
—0°270 
+ 0°273 
+O°122 
+-0°009 
-+0°023 
—0°'070 
—o'lg2 
—o'l3l 
—0'007 
—0°'075 
--0°807 
—O'124 
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TaBLe V_ (cont.) 
No. | B—V | My Tob Ttr My | My a | y | z u | v | uw 
Type Go V P . s 
I 0°58 | 4°83 | 0-029 + 0°087 0'084 8/— 5 33 6 — 19 
20 | + 0°57 | 4°80 | 0°023 +O°1Ig o-o10 9 22 37 18 15 
36 | +0°58 | 4°83 | 0-028 -0°005 O'110 13 12 31 | — 15 Il 
39 | +0°57 | 4°80 | 0'020 ~0°034 | —O'114 5 33 37 | — 290] — 20 
7° 0°59 | 4°86 | 0-031 | 0'055 | +0°094| + 0°273 11 | —26 | —16| + 50 23 
85 0°53 | 4°63 | 0023 + O°O10 o'120| +21 -37 3i-— 26 
94 | +0°57 | 4°80 | 0°077 | 0°057 | —0°293 1663 5 12 I + 130 68 
IIO | +-0°§2 | 4°58 | 0°084 | 0°073 0°392] + 0°241 3 11 5 32 26 
125 | + 0°50] 4°47 | 0048 0083 | ~ 0-099 2 19 9 12 3 
142 | + 0°50] 4°47 | 0°035 -0'128 0°092 II 25 + 13 8 
193 | 10°59 | 4°86 | 0°027 | 0035 0°054 0°270 35 = 10 25 2 
196 0°57 | 4°80 | 0°031 | 0'030 0°448 O'1gl 26 18 7 2 82 
204 0°56 | 4°77 | 0:036 0'166 0°292 2 13 5 64 85 
224 | +0°52] 4°58 | 0°044 0°133 0-162 23 2 3 5 9 
229 0°52 | 4°58 | 0:026 0°O17 0109 36 5 13 58 29 
244 0°55 | 4°73 | 0°033 Oro! o-o81 27 2 15 4 17 
274 | + 0°58 | 4°83 | 0°046 | 0-048 —0°532 0°975 16 I 15 2 103 
Type G1 V 
42 | + 0°58 | 4°83 | 0082 | 0-081 1°324 0664 | I 8 9 7O 46 | 
59 | + 0°59 | 4°86 | 0-020 | 0'023 0°265 | —0°047 II 37 | —32 15 85 | 
118 | +0°58 | 4°83 | 0-049 0°049 | —0°364 0056 I 19 7 33 44 
159 | 1+ 9°53 | 4°63 | o-o19 0°090 | —o-112 38 34 14 33 ut 
299 | + 0°59 | 4°86 | 0°078 | 0077 0°416 0662 3 10 27 40 | 
309 | +.0°56 | 4:77 | 0-057 | 0-060 | —o-or2 | - 0°324 9 2 15 13, | — 28] 
Type G2 V 
R2 | +0°56 | 4:77 | 0-128 0°134 1-688 1°163 3 4 6 66 ° 
4! 0°63 | 4°92 | 0°075 | o-105 | + 1°330| ~ 0669 I 9 10 77 50 
127 | ~ 0°58 | 4°83 | 0-045 | 0-042 | —0-517 | —0-133 | — 3 19 iI 39 34 
170 0°66 | 5:00 | 0-020 0°183 0070 39 | —30 10 7° 4 
176 | + 0°63 | 4:92 | 0-071 0'059 | —1°625 | —o-271 | —12 8 2 118 43 
199 | +0°68 | 5-10 | 0-028 + 0°034 | —o'107 35 6 6] — 13 12 
211 | +0°66 | 5:00 | 0-025 ~O'O10 O'125 40 4 1} + 20] — 19 
241 | + 0°59 | 4°86 | 0°035 + 0008 0°265 26 4 II 39 39 
251 | + 0°63 | 4°92 | 0-059 0'047 0°837 0681 | —12 8 9 69 29 
269 | +.0°63 | 4:92 | 0-016 —0-'061 | —0-200 48 3 | —40 - 23 61 
282 | +0°65 | 4:97 | 0-030 +0°208 | —0-073 | —ar 8 | —25]4 41 19 
290 | +0°59 | 4°86 | 0-071 | 0:068 0°142 | —0'297 9 2 II Ir | — 20 
305 | +0°62 | 4-91 | 0-025 | 0-013 | + 0°459 | +-0'009 17 7 36 | + 55 | — 10 
324 | +0°69 | 5-15 | 0°037 0°035 | +0°216| —o-150] —12 10 22 rr | — 36 
Type G3 V 
19 | +0°69 | 5-15 | 0-023 | 0068 | +0-309 —0°277 | 5 4 -43 a7 | - ‘9 
198 | +0°63 | 4-92 0°037 *0°0§5 | —o'100 | — 25 10 o| — 59 | — 28 
| 
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TABLE V_ (cont.) 


Fundamental data for southern stars 
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| 
No. | B-V | My | =m, | mH | M, | M, | - | , r ‘ . ‘a 
| 
Type G5 V 5 > 
| R4 | +0°70 | 5°20 | 0-155 | 0-156 3°047 0°744 I 4 5 33 95 21 
| 29 | + 0°72 | 5°30 | 0-062 | 0-091 0:230 0°452 5 3 15 14 35 12 
| 68 0°77 | 5°64 | 0°055 | 0-040 0°147 0°093 I 15 it) 8 4 18 
| 71 0°77 | 5°64 | 0°046 | 0°027 o-118 0°471 6 18 11 43 14 | —23 
| 79 0°71 | §°25 | O°109 | O-115 O'121 0:218 2 8 4 18 29 12 
95 0°74 | §°43 | 0°040 | 0:007 0:064 0°275 I 2 7 31 31 ° 
| 103 | +0°74| §°43 | 0°044 | 0°042 0002 0286 4 22 I 22 9 -20 
139 | 0°64 | 4°94 | O-102 | o-092 1°538 | © 0°393 3 9 4 65 39 6 
| 153 | 0°77] 5°64 | 0-020 0:176 0164 32 37 12 29 41 ~27 
| 174 | 0°61 | 4°90 | 0°038 0'022 O'110 22 14 5 4 3 14 
| 177 0°70 | 5°20 | 0°028 O'rdss 0°020 32 12 12 I 28 21 
| 215 0°71 | 5°25 | 0-040 0°22 o'1s! 24 6 3 28 26 14 
| 216 0°70 | 5°20 | 0'105 | 0071 0'021 0°197 9 3 | 2 | 13 4 I 
| 261 | 0°73 | 5°36 | 0°057 | 0051 0°022 0°532 14 3 10 4 43 7 
| 268 0°67 | 5:05 | 0-042 0°027 0°200 18 5 14 | 11 | 2 10 
278 0°70 | 5°20 | 0°053 | 0°038 0°483 0°429 12 6 13 54 | 29 6 
| 297 0°65 | 4°97 | 0°026 0°177 0:062 | — 22 Oo} —32 20 17 24 
| 311 0°73 | 5°36 | 0°032 | 0°035 0°369 0209 | 14 I | 25 63 16 9 
| 321 | +-0°69 | 5°15 | 0027 0°251 0'090 | 17 12 31 44 I 17 
Type G6 V 
182 | + 0°70 | 5°20 | 0°068 | 0-069 0°428 | —o-212 | 13 | 6 | : 19 | 27 6 
Type G7 V 
84 | | 0°77 | 5°64 | 0031 | 0-065 | 07035 | — 0°178 | 31 i -s8 20 | 35 27 
Type G8 V 
209 | --o-8o | 5°85 | o-119 | 07125 | 1010} 07190 | 8 | 3 | I | 29 | 21 | —31 
isa | sore lesol?™7| — | -ragel rags —21-2/-2|/ 7 a7|- 81-5 
254 | 0°75 15°50 a 0'170 | I al I rad + | 3 | 3 47 | 12 ~13 
| 308 | 0°80 | 5°85 | 0049 0050 | 0°363 | 0°034 9 |} 2 18 39 | 8 rs 
| 326 | ~0°78 | 5-71 | 0-027 | 0-041 | 0°196 | — 07300 | ~— 13 | 7| 34 | 16 59 11 
Type Ko V 
2 | + 0°95 | 6°63 | 0°054 | O-014 0°287] +O'11l 2 I 18 26 3 -— 9 
27 o'81 | 5-92 | o-og1 | 0°083 2°142 0662 fe) 5 10 101 76 | —25 
32 0°94 | 6°59 | 0°074 0°128 0°045 2 6 12 14 — 25 “= $7. 
45 | +0°89 | 6°35 | 0°075 | 0060 0°213 0'286 4 7 II 31 16 7 
86 o'gI | 6°45 | 0-088 | 0:063 0 198 0°052 6 9 2 6 14 8 
132 080 | 5°85 | 0-095 | 0°099 0°679 0818 3 9 5 52 20 15 
135 | +0°84 | 6-10 | 0-026 0°169 0°O17 14 35 8 —. ~1I 
143 | 0°82 | 5°98 | 0-059 | 0°070 | - 0°537 | —0°306 8 15 2 27 ~ i -= 
0046 ae ake = 18 12 2 3 | —142 | —42 
189 | ~o'80} 5°85} * wena $72 1°397 13 9 : 58 105 | —31 
200 0°91 | 6°45 | 0065 0°104 | --0°033 14 6 I 25 it fae 
206 0°89 | 6°35 | 0:062 | 0-076 0081 0*600 14 8 4+] - 33 29 18 
259 | ~o-91 | 6-45 | 0-140 | o-116 1'247 | —0°213 6 2 4 68 | — 12 9 
281 0°94:| 6°59 | 0:067 | 0°063 0°249 O°IgI | 10 1 I! 6 12 20 
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TasLce V_ (cont.) | 

4 ae et % | Si. | 

No. | B—V | My | 7 | Ttr | Mh, My % y z u | v w | | 

| | 

Type Ki V - - es J 
102 | + 0°93 | 6°55 | o'100 | o-og1 0°312 0°347 | 24 884 Fs 24 Z | O | 
227 | +0°97 | 6°71 | 0°097 | 0°008 0°130 | —10 | 1}/-— 2 4 3 8 | 

} 

Type K2 V | @ | 
28 | + 1°02 | 6°88 | 0-092 | | -; 0-094 | —o-o16 | 3 | 2 10 | 13 12 | 36| @ 
52 | +0-91 | 6°45 | 0°037 | 0°023 | 0°195 | 0°630 | 7 -. 20 87 8 7 F | 

284 | +0°89 | 6°35 | 0°069 | 0069 | +.0°424 | —0'208 9 | 2 II 5 22 | 37 | 

295 | +0°86 | 6-21 | 0°035 | 0-049 0°377 0°258 | rz | 2| —23 52 - 40 | 4 | | 

| 

Type K3 V | 
II | +0°94 | 6°59 | 0°077 | 0-103 | + 0°615 | + 0-040 | o| —1| —13 33 | 18 12 | | 
161 | +1°04 | 6°92 | 0-113 0:008 | 0°455 | —0°816 | 6 | 6 ° 2 9|-—26| | | 
187 | +0°84 | 6-10 | 0-062 | 0'071 0°130 | +0°326 14 | — 9 = 3 13 22 | 

256 | +0°87 | 6°26 | 0-154 | 0°174 | ~ 0-442 1°555 | — 6 I 3 +118 53 | +49 | 

328 | + 1-01 | 6°85 | 0-089 | o-102 | +0115 0°737 | 5 | 6 $i — 31 71 | —24 | | 

Type Ks V | { 

33 | +1°05 | 6-94 | 0°038 | 0°013 | +0°067 | — 0-516 | 3 13 23 42 60 | —11 | 
44 | +:1°16 | 7°09 | 0°063 | 0-060 0°363 O°241 | 2 I! I! . 8 63 18 | 
55 | + 1°14] 7°03 | 0°076 | 0-069 0°763 | +0401 2 9 9 36 ‘ +46 | 

112 | + 1°18 | 7-15 | 0°079 | 0087 | + .0°446 | —0-433 o| —13 2 37 9] -— 6] 

137 | + 1°07 | 6-98 | 0-069 | 0°075 | — 0°660} + 0-241 5 13 4 41 27 8 

141 | +1716 | 7-09 | 0-106 | 0-096 1‘075 0627 3 9 5 18 7° 5 | 

208 | + 1°16 | 7-09 | 0°143 | 0-172 | —0°479 | —1°143 7 ° I I 40 4 

294 1°06 | 6-96 | 0°279 | 0°285 | +3:941 -2°569 2 I 3 79 4° 6 

301 | + 1°10 | 7°00 | 0°073 | 0070 0°226 0809 sj 2 ta} - 7 56 | —16 ¥ 

307 | + 1°14 | 7°03 | 0°068 | 0-065 | + 0°375 | —0:027 6 2 13 21 8 13 

313 | + 1°13 | 7°02 | 0°129 | 0123 | +0°321 | —o-'171 3 ee ance 2 8 | —18 

Type K7 V 

222 | +1°31 | 7°74 | 0-075 | 0-073 | +o-122 | —0-438| —13 | — 3 | 3 | 10 | 20 | —18 

Type Mo V 

258 | +-1°44 | 8-94 | 0°157 | 0-158 | -+0°778 | —o-152 5{|—1] 4 | 38|+ 2 3 

277 | +1°42 | 8-74 | 0-258 | 0-255 | —3°283 | — 1-146 3 | ° Li - = | ae a8 

Type M2 V 

319 | +1°50 | 9°54 | 0°273 | 0-273 | + 7 : et — 4 | o| — 3| +106 | 14 | —S$S 
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In addition to those given in ‘lable V, trigonometrical parallaxes from (9) 
and (10) are available for the following stars : 



































No. 7 No. 7 No. | 7 
| — _ | 
4 0°153 117 o'0b1 205 O'O14 
5 30 122 15 207 183 
6 7° 124 53 212 15 
22 62 131 85 230 II 
23 148 138 6 232 20 

24 145 6 234 
25 6 147 29 235 + 
26 10 148 38 243 24 
37 32 150 4 247 10 
40 70 151 7 249 12 
49 42 157 15 252 15 | 
51 I 162 32 260 35 | 
54 37 163 33 267 7 | 
| 57 38 165 15 272 43 | 
58 59 167 49 275 12 
76 29 169 6 276 6 | 
| 33 35 171 24 279 12 | 
8o 69 178 10 289 2 
92 50° 179 53 296 SI 
96 I 180 5 300 30 
98 19 181 9 306 9 
99 46 183 78 314 144 | 
104 29 186 15 320 12 | 
105 | Igo 22 330 21 
109 59 19! 25 333 12 
: <a 23 197 27 334 2 
11s O’oll 201 0°021 339 ool! 








The Royal Observatory, 
Cape of Good Hope: 
1957 July. 
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THE FORMATION OF POPULATION I STARS 
PART I. GRAVITATIONAL CONTRACTION 


W. H. McCrea 
(Received 1957 May 3) 


Summary 

The present paper is preliminary to the second part in which it will be 
suggested that the gravitational contraction that is primarily significant for 
the formation of Population I stars is that of the material destined to form 
a cluster of such stars. The relevant type of gravitational instability is that 
recently studied by Ebert and by Bonnor, which concerns material under 
external pressure. A more elementary treatment than theirs is here developed 
from a general form of the virial theorem. ‘This shows that the essential 
features of their results hold good under somewhat relaxed conditions which 
may be more like those encountered in applications. But also it shows that if 
we employ the exact results of Ebert and Bonnor we may have to regard them 
as applying to somewhat larger masses (up to about fifty per cent larger) than 
those to which they explicitly refer, on account of incomplete realization of 
the theoretical conditions. The present treatment appears to afford some 
additional physical insight into the results obtained by Ebert and Bonnor and 
their relation to other criteria for gravitational instability. This is aided 
by its extension to cylindrical and plane-stratified distributions, which had 
been shown by Ebert to provide instructive comparisons. Among more 
detailed results we exhibit the behaviour for varying temperature of isothermal 
gravitating material under given external pressure and we note some 
numerical results required for subsequent use. 





1. Introduction.—A theory of the formation in interstellar matter of 
condensations that might lead to the production of fresh stars has been discussed 
in an important paper by Ebert (1955). In Part II of the present work it will be 
suggested that a condensation of the sort immediately considered by Ebert 
comprises the material of a cluster, rather than of a star, but that the stars of the 
cluster may be produced by a repetition of the same basic process. Bonnor (1956) 
has also studied the process independently, though with regard to a different 
astronomical application. ‘The object of this paper is to review the process and 
to obtain some numerical results for subsequent use. 

Ebert and Bonnor investigate the gravitational instability of a mass of gas 
under an applied pressure. The type of instability discovered by them is likely 
always to be more significant than the types previously discussed by other authors. 
The reasons for this are stated more particularly by Bonnor. Certainly, it appears 
to be the type that is relevant to the applications here in view. However, the 
results given by the authors are those for an unstable state that is approached 
through a sequence of states of strict equilibrium. Before using the results for 
the proposed applications, we need to be satisfied that they remain valid, to a 
sufficiently good approximation, under less idealized conditions. ‘The somewhat 
elementary discussion about to be given is presented in confirmation of this. 
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Also it appears to give useful physical insight into the phenomenon as well as 
a better understanding of the employment of the virial theorem in this connection. 

2. A general form of the virial theorem.—Consider a body of gas that can be 
regarded as permanently composed of the same set of molecules or “ particles”. 
Let there be a Newtonian frame such that there is no mass-motion of the gas 
relative tc this frame at the epoch of interest, say fy. Let a typical particle of 
mass m have position vector r in the frame at time f, and let P be the force acting 
upon it. Its equation of motion is 


mr = P 
and we have 
1° d . , ws 
Fall) =25, (0) = 2h +2r.F, 


Hence, summing over all particles of the set, 
1d SN mr2 mre SY 
5 gp hme) = Umer +=r.P. (2.1) 
26 


We have <mr*?=/, the moment of inertia about the origin. At ¢, we have 
“mr? =7, the total thermal kinetic energy, since there is then no mass-motion. 
The quantity =r. P is the virial. ‘The summation in the virial is by definition 
over all particles of the set. But this is equivalent to summation over all forces 
acting in the system such that, if P be now regarded as a typical force, r is the 
position vector of its point of application. 

We assume that the only forces are those acting in collisions of the particles 
with each other or with particles of any other material present, together with the 
gravitational attraction upon the particles exerted by the other particles of the 
system or by other material present. 

The forces in collisions between particles of the system make no contribution 
to the virial because they consist of equal and opposite pairs having the same 
point of application. If the forces in collisions with other material may be treated 
as producing a pressure p at a boundary S, then their contribution is f, pr.dS, 
where dS has the direction of the inward normal. If pdv is the mass of the 
particles in a volume-element dz at r, and if g is the total gravitational force per 
unit mass at r, then the rest of the virial is [, pr.gdv, where V is the volume 
enclosed by S. Thus (2.1) becomes 

1d7J 

2 df? 
This is a form of the virial theorem somewhat more general than that usually 
quoted (Jeans, 1925; Chandrasekhar, 1939). 

Equation (2.2) applies to a body of gas released at time /, from rest in a given 
configuration and under given external pressure. So at t, we have 


=2F + | pr.dS+ | pr.gdzv (t=ty). (2.2) 
/S ~v 


dI/dt=o (t =t). (2.3) 
If the given state is an equilibrium state we must also have 
d?]/dt?=0 (t=ty), (2.4) 


and so the vanishing of the right-hand side of (2.2) gives a necessary condition 
for equilibrium. The condition is not, of course, sufficient. But if (2.3), (2-4) 
are satisfied, we shall call the state a state of quasi-equilibrium in the present 
context. Otherwise, if the given initial conditions are such that (2.2) yields 
d@I/d?®>o, <o, we take this as a criterion that the gas as a whole begins to expand 
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or contract, respectively. This is the criterion generally used in such work. But 
it has to be noted that (2.2) will not apply further to the motion itself, since once 
the expansion or contraction is under way this will constitute mass-motion which 
is not taken into account in (2.2). However, the equation may show that there is 
no equilibrium state to which the motion could lead, in which case we are able to 
conclude that the expansion must proceed indefinitely or that the contraction 
leads to collapse. 
If the material has uniform temperature 7 then 


27 =3MRBT/p, (2.5) 


where M is the total mass of the material concerned and yp is its mean molecular 
weight (where » =1 for atomic hydrogen). 
If p has the same value at all points of S then 


| pr-dS=p | Fd = = | divede = ~ 3pV. (2.6) 


If g arises solely from the self-gravitation of the material concerned then 
(Chandrasekhar, 1939) 


| pr-gdo=0 (2.7) 


where {2 ( <o) is the potential energy compared with a state of complete dispersal. 


If all these conditions are satisfied, (2.2) becomes 


1d? : 
5 gp =3MBT|u-3pV+Q — (t=b). (2.8) 


For equilibrium or quasi-equilibrium this gives 
pV = MRT /p + 4Q. (2.9) 


In the case of equilibrium, if the last term is negligible, this reproduces Boyle’s 
law. If Q is not negligible, (2.9) expresses the appropriate modification of 
Boyle’s law. But we need to know more about the equilibrium state in order 
actually to evaluate Q. 

3. Isolated body of isothermal gas: temperature given.—We here confine 
attention to material at a given uniform temperature 7’, so that (2.5) holds good 
with the same value of 7 for any different states of the same material. 

We consider a body of gas in equilibrium under no forces other than its self- 
gravitation and a uniform pressure p over its boundary. This pressure is assumed 
to be exerted by other gases of the same or a different species. Then equation (2.8) 
is applicable. But also, because of the way in which the pressure is applied, the 
equilibrium state must possess spherical symmetry (cf. Lyttleton, 1953). There- 
fore, on account of its physical dimensions, Q is here of the form 

= — AyM?/R, (3-1) 
where y is the gravitational constant, R is the radius of the configuration, and 
A is a positive dimensionless number depending upon the density distribution. 

The quantity A is the same for homologous states of bodies of different mass. 
It is not necessarily the same for different states of compression of the same mass. 
But on general grounds we expect it to vary only slowly with R for equilibrium 
states and we shall provisionally neglect any dependence of A upon R. We shall 
verify later that this is in fact a good approximation. 
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Using (3.1), we find from (2.9) 
— 3 MRT p) = AyM? 


—4nR! gn R®’ (3.2) 
This is the modified form of Boyle’s law for the present case. 
It is convenient to define quantities R,, p,, p; having the dimensions of length, 
pressure, density, by 
__yM _ (BT \n)' (RT [us 
= (BT)? ge 3-3) 


Then (3.2) becomes 
R,\3 R,\4 
n3(R) ~4(R)- a 


We are considering the case of given M, 7, i.e. given R,, p,. If Ris sufficiently 
great, the second term on the right in (3.4) is small compared with the first. If 
then the pressure is slowly increased, so that the system may be regarded as 
traversing equilibrium states, R will decrease. It is seen from (3.4) how the 
proportionate effect of gravitation in assisting the compression increases as R 
decreases. In fact, there is a maximum value of the pressure for which any equilibrium 
state is possible. This behaviour is illustrated by Fig. 1, curve UU. 

The maximum occurs at R= R,, where, from (3.4), 








R.,/R, = 44/9, (3-5) 
and the critical pressure is given by 
PerlPi = i(44/9)~. (3-6) 
The corresponding value of the mean density p is given by 
Per/P1 = 3(44/9) *. (3-7) 


We may notice that this is four times the density that would be produced by the 
same pressure without the assistance of gravitation, in other words the density 
that a small specimen of the gas would have in the same surroundings. For the 
unmodified Boyle’s law in the present units is simply p/p, = p/p. 

If the pressure applied at the boundary is increased beyond the critical value p,, 
the system must collapse under its own gravitation. ‘his is evident from the fact 
that if p> p,,, equation (3.4) cannot be satisfied for any value of R. But it is also 
confirmed by going back to (2.8) and observing that if p>p,, then d*I/dt?<o 
for all R, i.e. the system cannot remain at rest in any state but must always contract. 
This conclusion does not depend upon the assumption that 4 is approximately 
constant; it follows from the fact that A has a positive lower bound. 

We see from (3.4) that for p less than the critical value there are equilibrium 
states having R<R,,; these are shown by the dotted part of curve UU in Fig. 1. 
But these states are unstable with respect to gravitational collapse. For, if the 
radius is slightly decreased below the equilibrium value, the presssure remaining 
constant, (2.8) gives d?J/dt?<o so that the contraction proceeds further, and 
there is no other equilibrium state in which it can be halted. If we retain the 
approximation of treating A as constant, even such unstable states exist only in 
Rpin<R<R,,, where Rmin is the value of R for which p=o in (3.4), i.e. 


Rwin/ R, = 44. (3.3) 
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For any smaller value of X the material would collapse under its own gravitation 
without the assistance of any external pressure, as can again be confirmed from 
(2.8). In fact, we have here the familiar condition for collapse that is inferred from 
the virial theorem in the form 

27 +Q<o (p=o). (3-9) 
But this more usual application of the virial theorem cannot be so significant as 
the preceding one. For the system cannot reach the point of collapse from Rmin 
without first being brought to the point of collapse from R,,. The essential 
physical difference is even more obvious from the cases considered in Section ro. 

Bonnor (1955, Pp- 353, 357) mentions briefly some consequences of the 
usual virial theorem, but he is not concerned to pursue them. Part of our present 
purpose is to show what sort of conclusions may legitimately be drawn from a 
suitable form of the theorem. Broadly speaking, we find it to give useful results 
so long as these are viewed as a whole. But it may be misleading to draw isolated 
inferences from it such as taking (3.9) as a generally significant condition for 
gravitational instability. 

It is to be emphasized that the discussion in this section concerns material 
that is maintained at a given temperature. The question as to whether this is the 
case required for applications will arise in Part II. 

4. Isolated body of isothermal gas: pressure given.—The properties of the 
critical state naturally do not depend upon its being approached through states 
of the same temperature. By way of illustration, we may consider instead states 
for which the applied pressure is given. 

Equation (3.2) may be rewritten as 

RT 47pR® AyM 
7 3M + 3R F (4.1) 





We write 
Ry =(4r)-M4yM4p-14 M22, RT| = (477) "4344 M2, 


pPo= (4a) Y4-**p4 1-28, (4.2) 
Then (4.1) becomes 
T 1/R\8 R 
7,-3(z) +48: (+3) 
The right-hand side of (4.3) has a minimum at R= R.,. where 
R.,/R_=(A/3)"" (4-4) 
for which the value of T is T., where . 
T >| T z= $(A/3)*" (4-5) 
leading to 
PerlP2 = 3(A/3)**. (4-6) 


It is easily verified that (3.5), (4.4); (3-6), (4.5); (3-7), (4.6) agree in giving the 
same critical state. 

The immediate interpretation here is that, given M, p, there is a minimum 
temperature below which no equilibrium state is possible. 

The general interpretation of the work is that (3.6) or (4.5), which may be 


written 
7 (RT.,/p)* 
Per = ap Soe , (4-7) 
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gives the locus of critical states in the 7, p-diagram so that for given M there 
is no equilibrium state for any pair of values of 7, p corresponding to a point 
above this locus (see Fig. 4). 

5. Gravttational collapse. Returning to the behaviour of material at a given 
temperature 7’, we have seen that for any value of this temperature, and for any 
value of the pressure p, the material will certainly collapse if it starts with a small 
enough radius. Now a point we have sought to make is that, so long as the 
material behaves in accordance with the modified Boyle’s law, it cannot reach 
such a small radius without first being brought to a critical point on (4.7) from 
which it would collapse anyhow. Nevertheless the question arises as to whether 
the material must inevitably obey the law—must it climb to the peak of the 
curve UU (Fig. 1) or can it in any way slip through? The answer appears to be 
that the gas must obey the law so long as there is no singularity in its density 
and that a singularity can arise only after the system as a whole has attained a 
critical state. 

However, this conclusion applies directly only to the case we have been 
considering of an isolated body of gas. If there are gravitating bodies immersed 
in it these can act as condensation nuclei (in the manner of Bondi’s symmetric 
accretion (Bondi, 1952; McCrea, 1956) although the steady-approximation may 
not be adequate). ‘Thus there appears to be no permanent equilibrium state 
in this case. But in the proposed applications of the theory the time-scale for 
the effect of condensation nuclei to become significant is almost certainly long 
compared with the times of other changes. So the effect may be ignored in these 
applications. 

6. Exact theory (Ebert).— Before making further use of the results of Section 3, 
we recall Ebert’s exact solution of the same problem. This is, of course, obtained 
by solving the equation of hydrostatic equilibrium under the appropriate boundary 
conditions. ‘The equation reduces to Emden’s equation for an isothermal dis- 
tribution and the problem is that of adapting Emden’s solution so as to exhibit 
the equilibrium of a given mass of material under a given applied pressure. 

Let {(&) be the solution of the Emden equation 


2 

wat agte= (6.1) 
such that {(0)=0, ¢’(0)=0. Then the properties of a mass M under pressure p 
are expressed in terms of £ as parameter and the quantities R,, p,, p, of (3-3) 
as follows: 


Applied pressure p=pye(&C')*, (6.2) 
Radius R=R,(-€/)"', (6.3) 
Central density Peet =Pill'C }*, (6.4) 
Mean density P = 3p1( — €2')*, (6.5) 
Surface density p=p,e'(€2Z’)?, (6.6) 
Potential energy Q = M(RT/p) —3-e€C—), (6.7) 
Virial coefficient A=3(-&)'-et'. (6.8) 


Expressions (6.2), (6.3) were given by Ebert. The others are readily 
obtainable from his work: in particular the expression (6.7) for Q follows most 
easily by using (2.9), and then (6.8) follows by using (3.1). Itisto be emphasized 
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that to any particular value of € there corresponds a particular configuration that 
can exist under the applied pressure given by (6.2) and that the remaining 
quantities characterize the configuration as a whole; the formulae do not express 
the variation of any quantity through the configuration. It may be noted that 
small values of £ give small values of ¢ and —¢’ and consequently small values 
of p and large values of R. 





1 1 J T 
400F\G 4 
p 
hp, 

B 
100F- = 


40r 
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0.2 0-5 1.0 


Fic. 1.—Relation between radius R and external pressure p for a given mass of gas at given 
temperature : EE according to solution of Emden equation (exact theory of Ebert and Bonnor); 
UU for uniform density according to the approximate theory of the present paper; BB according 
to Boyle’s law neglecting gravitation. The critical points on EE and UU are the maxima on 
these curves; the broken parts represent states that could not be attained in practice; the point G 
on BB corresponds to the radius from which gravitational contraction would occur according to a 
common interpretation of the virial theorem. (Scales are logarithmic.) 


The behaviour resulting from these formulae is obtained by using Emden’s 
(1907) numerical solution for ((€). In the significant range of configurations 
it is found to be qualitatively like that predicted by the elementary theory in 
Section 3 (and we shall see later that it is also quantitatively very similar). This 
is shown by the curves drawn by Ebert and by Bonnor and the authors’ discussion 
of these curves, Part of the p, R curve is reproduced in Fig. 1; a few numerical 
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values are given in Table I, chiefly to show the approach to the critical state given 
by the last line. 


TABLe [ 


Isothermal sphere, temperature T, under external pressure p 





| 











a | : | 
| 2 | I fixed _ fixed faaid | A 
| | Pip, R R, | 1 T; R R, 
| 2 | 2°04 +~=«+1°06 | 0°837. 1°27 1°41 0°616 
4 | 13°0 0°506 | 0°527 0°'961 2°70 0°666 
6 | as o'418 | 0489 0855 5°02 | 0°720 
65 | 17°6 o-410 0°488 0-840 5°78 | 0°733 | 
Ls | | | 


This exact theory again shows the existence of a critical state corresponding 
to that given by (3.5), (3.6). For given M, T, the critical pressure is that given 
by the smallest value of € (0) for which dp/df=o, that is from (6.2) 

eG (20" + 02+ 42-10) =0. 
Using (6.1), this gives 


¢'2—2e' =o. (6.9) 
From Emden’s numerical solution of (6.1) we find that (6.9) yields approximately 
£=6°'5, e* = 0-070, ' = —0°375. (6.10) 
Hence we obtain from (6.2)-(6.8) the critical values, approximately, 
Re, = 0°41 R,, (6.11) 
Per = 17° OP, (6.12) 
Per = 43°4P1» (6.13) 
(Peent/Per = 5°8 (Pcent! P)er = 14°35 (6.14) 
A,, = 0°733- (6.15) 


The values (6.11), (6.12) agree with those given by Ebert. 

From (6.12) the relation between critical values of p,, for given M, corre- 
sponding to (4.7) is 
(BT) 

Using (6.15) in (4.7) would give a coefficient 21:7 in place of 17-6; this is 
some indication of the error introduced by treating A as constant, as distinct 
from using an approximation to A,, itself. As another illustration (3.5) would 
give 0°33 in place of the o-41 in (6.11). 

We also note that the elimination of M between (6.11), (6.12) gives 


Per =17 


2 RT/p)° * 
Ror "Per = 2-96 = yry (6.17) 
Finally, for given M, p and using the definitions (4.2), we find 
R=R,(eé ¢’-*)"4 (6.18) 
T= T,[e(2 CP). (6.19) 


‘These together give in parametric form the relation corresponding to (4.3). ‘I'he 
39 
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relation is shown in graphical form in Fig. 2. As betore, these results lead to the 
same critical states as those got from given M, T. 

7. Approximate theory : the parameter A.—The theory given in sections 3-5 
would be exact if we were to use the value of A given by the exact calculations of 
equilibrium states. But in that case there would be no object in having the former 
theory at all; the results of section 6 would supply everything we need. We are 
going, however, to make wider use of the former theory. 
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Fig. 2.—Relation between radius R and temperature T for a given mass of isothermal gas 
under given external pressure : EE according to solution of Emden’s equation ; UU for uniform 
density according to the approximate theory of the present paper ; BB according to Boyle’s lau 
neglecting gravitation. The critical points on EE and UU are the minima on these curves ; the 
broken parts represent states that could not be attained in practice. (Scales are logarithmic.) 


Asa preliminary, we see from Table I that A does vary only slowly in the signi- 
ficant range of states preceding the critical state. Also we have seen from a couple 
of illustrations that treating A as not varying at all would not introduce a consider- 
ableerror. ‘Thus we verify that the approximate theory with A treated as constant 
would have given a reliable first approximation in circumstances in which we do in 
fact know the exact results. 

Now in order for collapse to occur as calculated according to the exact theory at 
some particular values of p,,, T,.,, just before collapse the system must be in the 
corresponding equilibrium state. In particular it must have developed the degree 
of concentration of its mass towards the centre indicated by equations (6.14). 
But we wish to apply the theory to a cloud of interstellar matter acted upon by 
external pressure arising in ways we shall describe. The requirements which we 
shall envisage occur in what must be regarded as rather rapidly changing conditions. 


—— 








NC nS «nbn ana nada el 














aes « 





SS ALE AE Sonate 





are Rasim arent 





No. §, 1957 The formation of Population I stars $71 


So it would be an over-idealization to suppose that, if gravitational collapse is 
produced, it is preceded by the establishment of a strict equilibrium state. Rather 
than to suppose the cloud to possess the degree of central condensation required for 
this, it would be more realistic to ignore the possibility of systematic density- 
gradients and to treat the density as being approximately uniform. 

We are interested here only in a cloud that does eventually collapse. In the 
first place this means, as we shall see, that the cloud must be roughly spherical. 

In the second place, we should not apparently be giving the process over- 
favourable treatment if we suppose the compression, prior to collapse, to proceed 
through a sequence of quasi-equilibrium states. For we thereby ignore the 
mass-motion associated with the compression and this mass-motion is in the sense 
that would actually assist ultimate collapse. 

For these various reasons we shall now write down results corresponding to 
those of sections 3-5, taking the system to be merely a uniform sphere of gas satisfying 
(2.9) regarded as the equation of quasi-equilibrium. For the reasons stated, these 
results may be more significant for applications than those of the exact theory. 
But we shall see that, considering all the approximations involved in the 
applications, the differences between the two sets of results are not great. Thus 
we may also conclude that the exact theory is not misleading for the required 
applications in spite of its high degree of idealization. 

8. Uniform sphere.—The coefficient A for a uniform sphere is 3/5. In this 
case therefore the relations (3.5)-(3-7), (4-4)-(4.7) become approximately (the mean 
density being here the uniform density) 


R.,=0°27R, (=0°40R,*) (8.1) 
Per=39°5P: (=17'6p,*) (8.2) 
Per=158p; (=70p,*) (8.3) 
R,,=0°67 R, (=0'82 R,*) (8.4) 
Top =0'40T, (=0-4973*) (8.5) 
Per = 10 Ps ( =8-2 p,*) (8.6) 
a T iy 
er= 39°53 oar ( =1 ata) (8-7) 


The curve UU in Fig. 1 is drawn for the case of a uniform sphere. It is seen 
that the curve closely follows that for true equilibrium states up to nearly the critical 
point for the latter. But its own critical point occurs for a somewhat higher 
pressure and smaller density. This is because the compressibility of the true 
equilibrium configurations is assisted by the concentration of material towards 
the centre. However, the difference is not such as to be of much importance 
in the type of application we wish to make. 

In this case the elimination of M between (8.1), (8.2) gives 

RT/p)? 
Rotten RTP a 
‘The coefficient here is nearly the same as in (6.17). Consequently if we take a 
particular value of p,, here it will cause a mass M to collapse from a radius R,,, 
and in the exact theory it will cause a mass M* to collapse from practically the 
same radius. We find that a good approximation is 


M*=3M. (8.9) 
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This is shown by the fact that the values in brackets in (8.1)-(8.7) are got by replacing 
M by M* given by (8.9) in the definitions of R,, etc. so giving R,*, etc. It is 
seen that the coefficients of R,*, p,* in (8.1), (8.2) are almost identical with those of 
R,, p; in (6.11), (6.12). 

It is therefore of interest to compare the p, R graphs for a mass 3M according to 
the exact theory of section 6 and a mass M according to the theory of the present 
section. These are shown in Fig. 3 where it has seemed more natural not to use 
logarithmic scales as in Figs. 1,2. We see that in the region of interest, that is as 
the critical point is approached by compression, the curves are very similar. ‘The 
inference to be drawn from this is that in applications we may use the exact theory 
while estimating that the departure of actual conditions from the conditions of the 
theory could require the masses concerned to be increased by up to about fifty per cent. 
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Fic. 3.—Comparison of p-R relations for mass 3M according to solution of Emden’s equation 
(EE) and mass M for uniform density according to approximate theory (UU). 


g. Critical values.—Returning to Ebert’s results, if for given T, p we take the 
value M,,, say, of M determined by (6.16), then at temperature JT any mass 
exceeding M.,, will be caused to collapse by external pressure p. 

If we substitute for M,, in terms of 7, p in (6.11) we obtain the critical radius 
R,, for this mass at the point of collapse. | We can assert that no quantity of the 
material at temperature T under external pressure p can exist in equilibrium so as 
to have radius exceeding R,,. 

If we substitute for M,, in terms of 7, p in (6.13) we obtain the critical mean 
density f,, for this mass at the point of collapse. We can assert that no quantity 
of the material at temperature TJ under external pressure p can exist in equilibrium 
so as to have a mean density exceeding p,,._ We find, in fact, that 


Per = 2°46p (RT /p)" (9-1) 
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so that the greatest departure from Boyle’s law that can be produced with the 

assistance of gravitation before actual gravitational collapse gives a density nearly 

2} times the density in the absence of appreciable gravitation. (The elementary 

theory gave from (3.6), (3.7) the factor 4.) 
Following Ebert, we use for interstellar matter the value 
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B=1'5. (9.2) 

(The present writer previously used the value 1-4 (McCrea 1953, 1956).) ‘Then 
the numerical results may be conveniently expressed as 

M., =8-91 T?(p/k)? solar masses (9.3) 

R,, =2:85 T(p/k)!? _parsecs (9.4) 


fg = 2:46 T-\(p/h) 
where k is Boltzmann’s constant and 7,, is the particle-density corresponding to 
Per= Some illustration of the values given by (9.3) is shown graphically in Fig. 4 
mainly for the purposes of Part II. This figure serves also as the graph of (6.16), 
and qualitatively of (4.7). 


particles cm-* (9.5) 
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Fic. 4.—Relation between critical temperature and critical pressure for given mass M according 
to solution of Emden’s equation (M=1, M= 100, M=500 solar masses). 


For comparison, the criterion (3.8) or (3.9) applied to uniform material of 7 

particles cm-* and using (9.2) gives a critical mass and radius 
M,,=41'1 T??n-"? solar masses (9.6) 
R,, = 64 T'?n-'? _ parsecs. (9.7) 

10. Other distributions.—Ebert obtained instructive results by considering 
distributions other than the spherically symmetric one. We shall treat these by 
the elementary theory of section 2 for comparison with his work. 

Cylindrical distribution. Consider a cylindrical distribution of gas in equili- 
brium under no forces other than its own gravitation and a uniform pressure p 
over its boundary, this pressure being exerted by further gas of the same or a 
different species. All properties are supposed not to vary in the direction of the 
generators of the cylinder, which is therefore regarded as of effectively infinite 
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length. For equilibrium there must be cylindrical symmetry. Let s be the radius 
of the boundary and m the mass of material per unit length of the cylinder. 

We apply the virial theorem to the material in a length / of the cylinder between 
right sections. Equation (2.5) holds good with M replaced by m/l. The first 
equality in (2.6) no longer holds good because the pressure is not uniform over the 
plane sections. But we may take it as defining a mean pressure p’, say, over the 
whole boundary of the region, so that the final member of (2.6) becomes ~— 378/p’. 
Equation (2.7) is not applicable because the field g arises partly from the rest of the 
cylinder outside the region considered. But from dimensional considerations we 
have 


| pr. gdv= — Bym?! (10.1) 
» 


where B is a positive dimensionless number depending upon the density distri- 
bution. As in the case of the coefficient in (3.1), we provisionally neglect any 
dependence of B upon s. 

The virial theorem (2.2) then becomes for equilibrium or quasi-equilibrrum 


0 = 3ml (RT /n) — 3781p’ — Bym*l 


p= (*- - \Bym) = (10.2) 


ms” 


giving 


If we ignore gravitation the second term on the right of (10.2) is absent and the 
pressure over the whole boundary is p; we should have again the standard form 
of Boyle’s law. ‘Taking account of gravitation, (10.2) gives the modified form of 
Boyle’s law for the present case. 

We see from (10.2) that no equilibrium state is possible at any radius if 

m> “en , (10.3) 
Otherwise, whatever the value of p’, there is always an equilibrium state for which 
the radius s is given by (10.2). This follows from the boundedness of B without 
assuming its constancy. Thus in the case of a cylindrical distribution the material 
can never be brought to the point of gravitational collapse by compression at any given 
temperature. 

Like (3.8), the condition (10.3) is equivalent to the familiar condition (3.9). 
But we now see how the fact that (3.9) happens to give the only condition for 
gravitational collapse in one case, the present one, is no indication that it is the 
significant condition in other cases, e.g. the spherically symmetric case. A similar 
remark applies to the use of the criterion for gravitational instability obtained 
by Jeans. 

Ebert’s treatment. Ebert gave a concise rigorous proof that 

pos? 
corresponding to (10.2). His work shows that in this case the various states of 
compression are homologous so that actually B is constant and p’ is a constant 
multiple of p. 

However, Ebert appears not to have remarked upon the existence of an upper 
bound for m for equilibrium, that is upon the property associated with (10.3). 
This property arises because a cylindrical distribution of gas at given temperature 
in equilibrium under its own gravitation alone, and with no singularity in the 
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density along the axis, is of unique finite mass per unit length. The mass that can 
be in equilibrium under an applied pressure must be less than this amount. If we 
attempt to add mass in excess of this amount, the system will collapse. 

This accounts for the difference from the spherical case. For there the dis- 
tribution in equilibrium under its own gravitation and with no central singularity 
is of infinite mass. Consequently, any finite mass of sufficiently large radius can 
be held in equilibrium by an applied pressure. However, in that case, the various 
states of compression are not homologous and so, as we have seen, gravitational 
collapse will occur when the radius is made sufficiently small. 

The writer is not aware of any published proof of the stated property of the 
cylindrical distribution. In Ebert’s notation, where € is proportional to the radial 
distance from the axis and » is proportional to the density at é, the appropriate 
“Emden” equation is 


1d (Edy 
rp(*g)- eat, (10.4) 
Also the mass within distance € is proportional to 
d 
_ “a . (10.5) 


Making the substitutions 


(10.4) becomes 
d*y/dr* = e* (10.6) 


and the expression (10.5) becomes 
2 + dy/dr. (10.7) 


If we suppose that y exists as + -> 00, then it can be seen from (10.6) that the 
only consistent possibilities for x, dx/dr are 


x-> 0, dx/dr-—> constant (7-> 2%). 


Hence either the distribution does not tend to infinity or the expression (10.7) tends 
toa finite limit. This is sufficient to establish the required property. 

Possibility of break-up. We wish to make some estimate of the stability of the 
distribution against longitudinal disturbances. _ It will be adequate here to use the 
approximation of treating the undisturbed density p as uniform because a radial 
variation of p will not have a direct effect on the problem. 

We have then 

m = 7rS"p (10.8) 


Ata point distant g from the axis of the cylinder g is directed towards the axis and 
\g|=2ymg)s? (q<3) (10.9) 


It is easily seen that the integral in (10.1) is independent of the choice of origin of r. 
If we take this to be a fixed point on the axis, the projection of r upon the direction 
of gis — q, SO 

r.g= —2ymq"/s* 
Thus (10.1) becomes 


- | pr. g dv = 2ymp(1/s*) | amg? dg = wympls* = ym*1, 
ov 0 
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using (10.8), so that in this approximation 


B=1. (10.10) 
in this case (10.3) may be written, again using (10.8), 
> 3 (87) (10.11) 
“me yp 


|For comparison, (3.8) leads in the case of uniform density to the analogous 
condition 
R®> 15 (BTu) (10.12) 
47 =P 
‘Thus, as might be expected, (10.11), (10.12) give roughly the same critical distance. | 
Now the Jeans criterion (Bonnor, 1956, equation (3.7)) gives the condition for 
gravitational instability for isothermal disturbances of wavelength A 


aie (BT /u) 
: yp 


If we apply this to disturbances propagated along the cylindrical distribution, it 
indicates that any tendency for the cylinder to break up by gravitational action 
would produce portions having length of the order of A given by (10.13). But 
(10.13) is derived from the theory of propagation of plane waves in an infinite 
medium. It can be a significant approximation in the present case only if the 
lateral extent is at any rate comparable with A. However, (10.11) shows that the 
medium could not then be in equilibrium but must in fact be collapsing laterally. 

On the other hand, if the radius is considerably less than the lower bound of s 
given by (10.11) and, consequently, than the lower bound of A given by (10.13), 
then the gravitational efficiency of longitudinal disturbances would be reduced. 
Thus the significant lower bound of A would be even greater than that given by 
(10.13). Hence any gravitational break-up of the cylinder in this case could 
result only in portions that are themselves still of length large compared with their 
radius. 

Combining all the results of this section we conclude that a long uniform fila- 
ment of intersteller material, provided it can withstand any lateral pressure at all, 
can withstand any amount of such pressure and that, moreover, it is gravitationally 
stable against break-up into short portions. Had the latter part of the conclusion 
not been reached, the first part by itself would not possess much significance. 

Plane-stratified distribution. Consider a plane stratified layer of gas in equili- 
brium under no forces other than its own gravitation and a uniform pressure over 
its two parallel plane boundaries. All properties are supposed not to vary in 
any direction parallel to these planes. For equilibrium, the distribution must be 
symmetrical about a median plane. Let 22 be the thickness of the layer and let o 
be the mass per unit area. We apply the virial theorem to the material inside a 
volume /x/x 2z of the layer. As in the previous case, we may define a mean 
pressure p’ over the boundary such that the pressure term in (2.2) is again 
— 3p’ x volume, that is — 6p’/?z in thiscase. Analogously to (10.1) we have 


(10.13) 


‘ pr. gdv= —Cyo%l*s, (10.14) 


where C is dimensionless. Then corresponding to (10.2) we obtain 
0 = 30/°( RT.) — 6F*2p' — Cyol*z 
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giving 
p= — — 4 Cyo’. 


as (10.15) 

Ignoring the gravitational term in (10.15) we again get the elementary form ot 
Boyle’s law. But we also seg in this case, unlike that of a spherical distribution, 
how even allowing for gravitation the material behaves more closely in accordance 
with this law the greater the compression, i.e. the smaller the value of z. This is 
because for an infinite plane-stratified distribution the gravitational attraction 
at distance q, say, from the median plane depends only upon the amount of material 
between +g, and not upon the value of q itself. ‘Therefore gravitational effects 
are not increased by compression, whereas the thermal pressure is increased. 

Interpreting (10.15) as in previous cases, this equation suffices to show that 
there must be a value of z yielding an equilibrium state for any applied pressure. 
Then the material can never be brought to the point of collapse by compression, 
as we should also expect from the argument just given. 

In this case there is no possibility at all of gravitational collapse of the sort we 
have seen to be possible for a spherical distribution. This again illustrates the 
untrustworthiness of attempting to argue from the behaviour of gravitational effects 
in one dimension to that in three dimensions. The possibility of genuine gravita- 
tional collapse, as distinct from some sort of contraction from an unrealizable 
initial state, seems to have been first made clear in Ebert’s discussion of the three- 
dimensional case. 

An exact solution for the present plane-stratified case was also given by Ebert. 
It confirms the inferences we have drawn from the virial treatment. Actually the 
latter may not be pushed much further than we have taken it because it is nota good 
approximation to ignore the variation of the coefficient C over a large range of =. 

11. General interpretation.—Interstellar matter is often observed in the form of 
elongated distributions and probably also in “‘sheets”’. It is sometimes remarked 
that these give the appearance of being under compression by other parts of the 
material. We now see from Ebert’s work and the present discussion that gravita- 
tional collapse cannot occur directly from these forms. In any case, the distri- 
butions that are most commonly recognized as filamentary or sheet-like are not 
thought to have any immediate connection with star formation. Also their 
behaviour is not likely to be much influenced by self-gravitation. The only 
interest of the present work in relation to such distributions is that itconfirms their 
gravitational stability even if their self-gravitation is significant. 

Nevertheless, the bodies of gas concerned in star-formation probably arise 
from temporary distributions that have either an elongated form, as the “ elephant 
trunks” considered by Ebert, or a flattened form, as the compressed regions 
considered by Oort (1955). These undoubtedly break up quitt rapidly. But 
their break-up must be owing to irregularities in the density-distribution and in the 
external pressure to which it is subjected, or to differential motions. Our 
discussion, based on Ebert’s work, shows that gravitational effects cannot play a 
significant part at this stage. It shows that it is only after such break-up has 
produced a body of gas of roughly the same diameter in all directions that gravita- 
tional collapse may ensue. Whether it does so or not can then be determined 
by Ebert’s and Bonnor’s theory for spherical distributions. The present work is 
claimed to show that the use of their results is not vitiated by a lack of exact 


realization of the conditions under which they are obtained. 
40 
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Since this was written, a paper by Mestel and Spitzer (1956) on star formation 
has been published. While their paper is concerned mainly with magnetic 
effects which are not dealt with here, it does make extensive use of the virial 
theorem and conditions for gravitational contraction derived therefrom. ‘Those 
conditions are open to the criticisms mentioned in section 3 above, although any 
consequent amendment would probably not substantially affect the authors’ 
conclusions. 


A preliminary account of the present work was given during the spring of 1956 
inlectures in the Berkeley Astronomical Department of the University of California 
where the writer was privileged to hold a visiting professorship and also a 
Fulbright travel grant. ' 


Royal Holloway College, 
Englefield Green, Surrey: 
1957 April 20. 
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RADIAL VELOCITIES OF SOUTHERN B STARS DETERMINED 
AT THE RADCLIFFE OBSERVATORY* 


(Paper II) 


M. W. Feast, A. D. Thackeray and A. }. Wesselink 


Summary 

in the second stage of the Radcliffe programme on radial velocities of 
southern O and B type stars, velocities of 130 stars selected for distance 
(chiefly mg-..M> 11-0) are presented together with 59 stars with previously 
observed velocities. All the stars have been classified on the MK system. 

The fainter stars demanded the use of faster emulsions compared with 
Paper I, and an investigation of residuals of the lines He 4471, 4026 has 
shown that for these emulsions significant deviations from LAU wave-lengths 
occur which vary with line-width, in accordance with Petrie’s findings. 
The Radcliffe procedure as regards wave-lengths of 4471, 4026 has been 
revised accordingly. 

The stars of this paper with newly determined velocities are at an average 
distance more than twice as great as that of the stars in Paper I. 


*The full text of this paper is published in Memoirs of the R.A.S., 68, Part 1, 1957. 
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